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The Generalized Plane Strain of an Elasto-Plastic Material 


By Davib R. BLAND, Manchester, Great Britain?) 


Introduction 


The majority of solutions to problems in plasticity obtained to date have 
been for a rigid-plastic non-hardening material, generally in plane strain or in 
axial symmetry for which general methods of solution are known [5, 10] 2), For 
an elasto-plastic material there have been about a score of particular solutions 
into some of which workhardening has been incorporated and into a few of 
which the effect of temperature gradients [1]. In each of these particular solu- 
tions there was such a high degree of symmetry that the dependent variables 
were functions of one space variable only. It is the object of this paper to put 
forward two methods of solution for certain problems of elasto-plastic materials 
in generalized plane strain, i.e. the simplest case of dependence on two space 
variables; the first method is for non-hardening materials, the second for work- 
hardening. 

It has been found by experiment that most metals satisfy the Mises yield 
criterion and its associated flow rule—the principal exception is the initial 
yielding of mild steel— and not the Tresca criterion [11, 7]. However by a suit- 
able choice of constants the absolute error in the yield stress predicted by the 
Tresca criterion can always be made less than 8%. It has already been found 


[10] that for rigid-plastic materials in axial symmetry a general method of | 


solution is feasible using the Tresca criterion but not using the Mises criterion. 


The same situation will arise in this paper. The error in a particular problem’ 


‘can be made considerably less than 8% if the stress vector in the plastic region 
covers only a small area of the yield surface in stress space but it may exceed 
8% if the choice of yield criterion appreciably effects the position of the elasto- 
plastic boundaries. 

Temperature effects are generally included in the classical theory of elasti- 
city by an additional term in the generalized HooKe’s law. In fact it is more 
correct to say that the classical theory is a particular case of the more general 
thermo-elastic theory obtained by assuming either isothermal or adiabatic con- 
ditions. A concise treatment of this topic has recently been given by LessEN [6]. 
The specifically plastic equations of an elasto-plastic material, namely the yield 


criterion and the associated flow rule, are unaffected by temperature provided ~ 


1) Department of Mathematics, University of Manchester. 
2) Numbers in brackets refer to References, page 132. 
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that the variation of temperature is not large enough to affect the value of 
the yield stress. It is assumed in this paper that the values of all physical con- 
stants that occur are effectively independent of temperature. Temperature 
effects are therefore included in these theories of elasto-plastic materials only 
by the presence of the temperature term in the generalized HOoKE’s law. 

In the treatment of the non-hardening material, it is assumed that the rota- 
tion of an element is small. This assumption is implicit in most treatments of 
the theory of plasticity. If ,;; is the rotation tensor for a material element, if 
D/Dt denotes the time derivative following a particle of a tensor referred to axes 
fixed in space and if d/dt denotes the time derivative following a particle of the 
tensor referred to axes fixed in the material element, then for a second order 
tensor, see e.g. NOLL [9], 

de; ; De;; 


ear De + €nj Opi + Cin Mp; - 


It is the derivative de,,/dt of the strain tensor that has the greater physical 
significance because it is de® |dt, not De? /Dt, that appears in the theory of the 
plastic potential. The axes. used in particular problems are axes fixed in space’ 
and therefore the two rotation terms in the above equation should be included 
at some stage of the analysis. They are generally omitted, i.e. it is implicitly 
assumed that they are small. = . 

The theory for the non-hardening elasto-plastic material is similar to the 
well-known theory for the rigid-plastic material. In fact the latter is a special 
case of the former obtained by letting the elastic moduli tend to infinity. 

The treatment in the workhardening theory is limited to linear workharden- 
ing materials, i.e. materials for which the yield stress is a linear function of the 
equivalent plastic strain, and to small strains because the linear strain-displace- 
ment equation is used. The fractional error introduced by the latter assumption 
is of the order of the strain and, since the strain is small, the approximation of. 
linear workhardening is not unrealistic for actual materials. 

For a given total permanent deformation of a body, workhardening has the 
effect of distributing the deformation more evenly through the body, i.e. one 
tends to have a large area of small permanent deformation rather than a small 
area of large permanent deformation. This is illustrated by a tension test of a 
cylindrical specimen (Figure 1). First yielding will occur at some section, say 
AB, determined: by circumstances fortuitous ona macroscopic scale, and this 
yielding will produce a small reduction in area at this cross-section. If the 
material is non-hardening, the specimen is now weakest at this section and all 
further yielding takes place there. However if the material is workhardening, 
the increase in yield stress due to permanent strain more than balances the 
decrease in cross-sectional area and the yielding then occurs at another section. 
The over-all effect is that the cross-section of the cylinder decreases uniformly 


I 


. 
: 
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along its length until a strain is reached where the increase in workhardening 
no longer balances the decrease in area. 


Figure 1 


A tension test specimen. 


The approximation is made in the workhardening theory that the rotation 
of the principal axes of stress in any material element during plastic deformation 
can be neglected. The stress-plastic strain increment equations can then be 
integrated to give stress-total plastic strain equations. Without this assumption 
the subsequent analysis would not be possible. The limitation to small strains 
tends to reduce the error introduced by this assumption because in the limiting 
case of vanishingly small strains there is no approximation —the total strain 
has become an incremental strain. Workhardening also tends to reduce the 
error as will be described in the next paragraph. 

Consider a workhardening body in plane strain, subject to external forces of 
which the directions and the ratios of the magnitudes are maintained approximate- 
ly constant, and which deforms in such a way that the mean rotation of the body 
is zero, i.e. over any cross-section lying in a plane of plane strain, /w dS = 0. 

Then if any part of the body rotates it will produce strain in the surrounding 
parts. These strains if more than about 0-3°% will be plastic and workhardening 
will occur. Consequently greater stress will be required for the further rotation 
of the part of the body first considered. The effect of workhardening is to re- 
duce differential rotation between one part of a body and another and therefore 
if the mean rotation is zero, to reduce the absolute magnitude of the rotation 
in all parts of the body. It does not follow that the principal axes of stress will 


not rotate during deformation if the external forces satisfy the above conditions, 


i 


ras 


but the rotation under these conditions will be reduced. Consequently the rota- 
tion of the principal axes of stress in any material element is reduced if the 
material is workhardening and if the external forces are maintained as stated. 


The Yield Criterion 


If the z-axis is normal to the planes of plane strain of a plastic rigid material 
de’, = 0, and therefore the equation for de?, derived from the plastic potential 
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can be used to eliminate o,, from the yield criterion. Then for a non-hardening 
plastic-rigid material, the yield criterion and the equations of equilibrium pro- 
vide three equations for the three components of stress, 0, 7, Ory ando,,, which 
can be solved when the boundaries of the plastic region and suitable boundary 
conditions are known. If de?, had not been known in advance it would have 
been necessary to solve the stress and velocity equations simultaneously. The 
difference is so crucial that it is now considered under what condition does 
de?, = 0 in the generalized plane strain*) of an elasto-plastic material when the 
material satisfies the yield criterion either of von MIsEs or of TRESCA. 
The yield criterion of von MIsEs when o,, = 6,, = 0 is 


/ = (Fx rae, Oyy)” ats (Cyy a Gua) at; (Oe <o en)" a 0 Oy =6R?. (1) 


By the theory of the plastic potential [3] 


of 1 
de, = 55 dh=4 & —~5 Gre oyy)| da. (2) 
If de?, = 0, 
1 
O22 = ee (Cre ate yy) , (3) 
é,, = C and e?,=0 since de?, = 0 throughout the deformation. Hence ¢, = 


Eg, — ef, =C. 

The generalized HOOKE’s law in the presence of temperature gradients is 
1+y 

of er et ~ Onn Oi, + « T Oj; (4) 

and therefore 

EC=E &, =62—¥ Cag t+ Gy,) +EaTl. 
Whence 

Oz2 = ¥ (xe + Oyy) + E(C —aT). (5) 


Equations (3) and (5) are consistent only if 


1 
oy (6) 


a 
and 
ie ene 
ive Consta= ae (7) 


Equations (6) and (7) are necessary conditions for de’, = 0 for a von Mises 
elasto-plastic material in generalized plane strain. They are also sufficient as 
can be seen by substitution from equations (6) and (7) into equations (5) and (2). 


3) Generalized plane strain is defined by u = u(x, y), v= v(x, y), w= Cz. In plane strain C = 0 


ee 


OS fe 
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The yield criterion of TREsca is 
f= (Cpe = 6y,)2 + 4.03, =4 (8) 


whenever o, , is the intermediate principal stress. From the theory of the plastic 
potential and equation (8), 4 

det, ="! ano. 

oO, & 

Therefore a sufficient condition for de®. = 0 for a Tresca material is that o,, 
is the intermediate principal stress. 
If o,,4, and o,,;, are the algebraically greater and lesser principal stresses 
in the (x, y)-plane, then from equation (5) the sufficient condition is 


Onax > ¥ (Gx = Oy uv) a E (C rau T) = Omin » 
Or, using Cmax — Onin = 2k and Omax oc Omin = Ox x Oyy, 
1 il 
k + 2 (Cra + Fyy) > ¥ (Fx~ + Oyy) 35 E(C= aL) > hs (Fax + Fyy) » 


x. 6 

(5 -») Cee + vy) —E(C-aT)| <b. (9) 
Note that, when equations (6) and (7) are satisfied, the inequality (9) is always 
satisfied and equation (1) reduces to equation (8). It follows that the class of 
problems for a von Mises material in generalized plane strain for which de?,=0 
is included as a special case of the more general class of problems for which 
de’. = 0 for a Tresca material. For the rest of this paper it is assumed that the 
yield criterion is equation (8) and that the inequality (9) is satisfied at all points 
in the plastic region. 


Non-Hardening Materials 


The method of solution of a set of partial differential equations generally 
depends upon their nature, i.e. whether the equations have or have not real 
characteristics. If the characteristics are real and distinct (hyperbolic equations) 
it is convenient to refer the equations to the characteristics and to solve parti- 
cular problems by integration along them. If the characteristics are imaginary 
(elliptic equations), weak discontinuities of the dependent variables do not 
occur (except where the given variables are discontinuous) and it is immaterial 


~ in which direction the equations are integrated. It is well known that the equa- 


tions for the plane strain of a plastic-rigid material are hyperbolic and that the 


- two sets of characteristics of stress and velocity are orthogonal and coincident. 
In this section the nature of the equations for a non-hardening elasto-plastic 
- material is investigated, and in the next section the work-hardening case. 
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For a non-hardening material, the yield criterion equation (8) and the 
equations of equilibrium 


00 y x Cc, nae 06 ny Wop os. 1 
Fixit Roy en ane (10) 
are three equations containing only three unknowns o,,, 6,, and o,, because 
k = const. These equations are identical to those occurring in plastic rigid 
theory. It is well known ([5], chapter 6) that the equations (8) and (10) are 
hyperbolic with orthogonal characteristics. By substituting 


Orzg=—P—ksin2y, oy,=—ptksin2p, Gyy=Rceos2p (11) 


(p + 2/4 is the anticlockwise angular rotation of the direction of the maximum 
principal stress from the x-axis), the equations referred to the characteristics 
can be expressed in the form 


p+ Z2Zky = const »on a =tang, ana-line, 
(12) 


dy bs : 
p—2kgw=const on = tan(p +3), a f-line . 


Next consider the velocity field. The plastic strain rates are given by the 
theory of the plastic potential and by equations (8) and (11) as 


oo 5 =—kisin2g, 
ip = ov — See) * _ pi sindy, (13) 


éo =? oF yA=kicos2o. 


The elastic strain rates are obtained from the generalized HookE’s law equation | 


(4), after substitution for o,, from equation (5) and after differentiation with 
respect to time (the derivative following a particle) as 


pean aaccoasdcwwadgs re Read 7 
= (1+) (— (1-20) p— 2h cos2 g-g — 72 © + +EaT), 

SN a ee (14) 
ST Ra tet tae) ahh yi ot 
( ”) (= ( v) p+ 2Rcos2g-@ Tape ae 


Eesy = (1+) G,,=—2 (14+) ksin2o-@. 


a 
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The generalized plane strain constant C is a function of time only. 
The strain-rate velocity equations are 


Sires 20 SE RE ei a RE Ula 12) 
ne dx’ Eyy = dy anc 4 Egy = dx dy , 


(15) 
where U and V are the components of velocity parallel to the x and y-axes 
respectively. Then 

aus + a = Coxe + Eyy = Exe + Soy + ey + eb, - 
On substitution from equations (13) and (14) 


0U OV eae : 
a tt ee Pee aL TSO) (16) 


where the modulus of rigidity 


From equations (13) 


cos2 p (Gee — 4g) + sinZ io (Eng — Eos) = 9. 


Substituting from equations (14) and (15) 


0U 


Loe k : : : 
; on c 
gd. faere b+ GE cos2Zgy-ptyve (l+) a ) 


cos2 yp ( 


+ sin2 p [ (= + a +2 sin 2 p-9| = 0 


or 


0U ee OV 0U 1—2» ; Re. 
cos2p G+ > sin29 (5. + cai 26 cos2p- Pt EP 


| phe! | ig 
+ cos2y (vC — (1+)aT)=0. 


| Equations (16) and (18) are two equations for U and V. The characteristics of 


the equations are 


dy dy _ als 
Gy a tang and aa tan(p + ). 


‘The characteristics of stress and velocity coincide. If u and v are the compo- 
nents of velocity and ds, and ds, the increments of length along the a-lines 
_and f-lines respectively, then the equations holding along the characteristics 
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can be shown to be 


du —vde + (“Fer p+ egerrc—(ts)aT) ds,=0 onvane-line; | 


dv + udp + (ze? t p+rC—(L+%) aT) ds,=0 on a #-line. 


Equations (19) are the analogues of the Geiringer equations for the plastic- 
rigid material. When T = C = const and y = 1/2 (an incompressible material) 
the equations reduce to those derived by Hit ([5], p. 241). When G + oo and 
T = C = const, equations (19) reduce to the Geiringer equations. Alternatively 
one can say that the terms inside the brackets in equations (19) are the additional 
terms introduced into the Geiringer equations by including the elastic terms. 

In general it is necessary to follow the course of the deformation in solving 
equations (12) and (19), ie. it is necessary to work out a series of solutions at 
successive times. In practice this is prohibitively time consuming. As in plastic- 
rigid theory, however, there are certain classes of problems in which only one 
solution is required for a finite deformation. One class is steady flow for which 
the time derivative at a fixed point is zero. In this case the time derivatives 
in equation (19) become 
f+ 60 Tee 

OsB 


OSe 


: 0 0 ; 0 
b p Op P 


Re et jane g=u-—t+v (20) 


OSx 


The values of the partial derivatives of # and m are known from equations (13) 
and those of 7 from the solution of 4?7 = 0. Equations (19) provide two 
equations for u and v. 

In practice the elasto-plastic equations (19) may only be required when the 
elastic and plastic strains are of comparable order of magnitude. An important 
example is the cold rolling of metals for passes of less than 3% reduction for 
which no satisfactory theory at present exists. It should be noted that it is 
necessary to solve the elastic equations in the elastic regions and to obtain 
continuity of stress and displacement on the elastic plastic boundary. This is 
a more difficult procedure than that for plastic-rigid materials, but it seems 
reasonable to suppose that the lack of uniqueness, that occurs for plastic-rigid 
materials, will not arise for the elasto-plastic and that the elasto-plastic solution 


for steady flow problems will also give the residual stress distribution in the 
deformed material. 


Workhardening Material 
Nature of the Equations 


As stated in the introduction, the analysis is limited to small strains and 
to linear workhardening and the approximation is made that the rotation of 


‘ 
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the principal axes of stress in any material element, during the permanent de- 
formation of this element, can be neglected. 

The three components of stress in-the (x, y)-plane still satisfy the yield cri- 
terion, equation (8), and the equations of equilibrium, equations (10). But k 
is no longer constant and therefore it is no longer possible to solve equations 
(8) and (10) to completely determine the components of stress before consider- 
ing the components of strain. This is in contra-distinction with the situation for 
non-hardening materials, both rigid-plastic and elasto-plastic. p and gy will 
again be introduced by equations (11). This introduction is equivalent to satis- 
fying equation (8). 

The limitation to small strain means that the strain-displacement relations 
are 
Ou ov 1 / dv 
= Fe ene Aa a (5° i a) = Exy- (21) 
In this section « and v denote the components of displacement parallel to the 
x and y-axes. 


The yield stress of an element of a linear workhardening material satisfies 
k=k, (1+ 728), (22) 


where f, is the yield stress of the annealled material, 7 is a constant and é? 
the equivalent plastic strain undergone by the element since annealling. The 
author has shown [2] that equivalent plastic strain and plastic work done per 
unit volume, W,, are equivalent measures of workhardening for a material 
satisfying TREsCa’s criterion, provided that no yielding takes place at a singular 
point on the yield surface. This condition is satisfied in this paper since a, , is 
assumed to be the intermediate principal stress. It is more convenient to use 


_ é? than W, as the measure of workhardening in this analysis. 


“h* 


Consider axes OX, OY fixed in a material element. The plastic strain in- 
crements at any stage of the deformation, for a material having the yield 
criterion and plastic potential given by equation (8), are 


—ksin2@d/, 


~ 
x 
es 

| 
= 
> 

Q 
Le 
Lal 
dS 

I 


del,» = (Oyy — Sx x) Gp = Psin2 Pdi, 
de’. = Oxy A=kcos2Pdh, dee, =0, 
where equations (11) have been used with ® + 2/4 as the anticlockwise angular 


rotation of the direction of the maximum principal stress from the positive 
_ X-axis. The approximation that the rotation of the principal axes of stress in 


ASS 


Ss 
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any material element can be neglected means that @ is effectively constant 
and the plastic strain increment equations can be integrated. If the directions 
of the axes, OX, OY in each material element are chosen so that they coincide 
with the Ox, Oy-axes at the and of the permanent deformation, then 


eb,=—Aksin2y, e&,=Aksin2o, e,=Akcos2y, &,=0. (23) 


The elastic part of the components of strain are given by the generalized 
Hooke’s law, equations (4). Substitution for o,, from equation (5) and use of 
equation (17) give 


eo, = 5 a Por sin2g+(1+v)aT—yvC, 

a : 

y= 5G Pt ae sn2e rT i py aT oC (24) 
ee =. c0s2 


Substituting from equations (23) and (24) into equations (21), 


Ou ; 1—2 

a = ban = bea + the ——ge (A+, Gg) ksin2p+ (L+»)aT—vC, 
ov iy ‘ 1 — 

dy buy = Syy + yy = —— pan “pt(A+a, cg) hsin2@ + (1 +v)aT—yvC, 
1 / ov Ou 

Plaats | = oan = oy thoy = (A+ 3G) Reos2y. 


The equivalent plastic strain, ¢’, is defined by 


= Ves | (act, def.)ue 


the integral being taken from the annealled to the current condition. When, 
as in the present instance, the plastic strain components increase in proportion, 


then 
EP eh ie (ef, é?) 1/2 . 


Substituting from equations (23), 
Ee? = —_ )k. 
Ve. 
Substituting in equations (22), 


bs V3 (k — Ry) 
: “V2 Re k # 


Oe 
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whence 


1 1 /3 V3 
eae ie Ae oe 
2G 7 2G 27 ko 2nk- 


Substituting back into equations (25) for 1/2G +4, 


Ou a 1 
ee cue = —— a— b- es +35, -) ksin2.p + VS sin 


+(1+r)«T—vC, 


oO . 1-2» 1 y3 3 
Oy = uy Se b+ (ag + Hp) bsin2y — V3 sinag 


+(1+r«T—rC, 


1 /odv Ou 1 V 
a ae + 3) fay (ee eed 5) cos2p— V3 cos2g. 


The equations of equilibrium in terms of #, & and @ are 


ty) 0 : 0 

— Se — gp (bsin2g) + = (eos 29) =0, 
Op 0 0 : 

mee tk C082 p) (ksin2 y) =0. 


| (26) 


(27) 


Equations (26) and (27) are five equations for u, v, p, k and gy. The nature of 
the equations, i.e. whether the characteristics are real or imaginary, will now 


be investigated. 
From the first two of equations (26) by addition 


mi G Ou Ov 2(1+ )G 2G 
rP in (& =o Le ia as amas Erk 
and by subtraction 
Se Fy , 
yhsin2g = — 4 4 V3 snag, 
where 
Take 
Be ye 


The last of equations (29) can be written 


0 3 
ypkcos2p= - + = ae = cos2 9. 


(28). 


(29) 


(30) 


(31) 


Equations (28), (29) and (31) are used to eliminate p, k sin 2 and & cos 2p 
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from equation (27). This gives 


yG O7u Ou 2(1+ %)yGa_ oT 
1—2y» a Ox OY 1—2y» OX 
. 2y3 Op Aa ee ec Op 
1 724 Fs cos2 py -- Fit Ne a0 == (i, 
(32) 
yG O2u 02v 2(1+ v) yG a oT 
“Ss ror t H)- ee, ay 
DY Si OP 2/3 Op Bor 
bed weno ee Trt cos2 Oy On 
Multiply equation (29) by cos 2m, equation (31) by sin 2m and subtract: 
Ov Ou Ov Ou 
— —- = ——  —— 33 
oy ” = tan2 (52 + aa (33) 
Differentiating equation (33) with respect to x, 
02v Ou 02 OP \ a tae Ou\ Op 
“ae ere Ae + ae aS 29 (s+ oa On 
Using equation (31), 
0? O*u 07u 02u V3 0p 
Ox Oy Ox? tan29 (Gea + Ox wl =2(ya—W) sec2 p > 
similarly (34) 


02u 02u y3 
Op Saray — tan2@ (sry te oat =2(ya— 1 Fe) seca $e oy” 


y can now be eliminated from equation (32) using equations (33) and (34). 
However, since no further differentiation is required for the present purpose 
of determining the nature of the characteristics at any point in the field, first 
choose any one point,.orientate the axes so that gy = 0 at this point, and then 
substitute from equations (34) into equation (32): 


pS pices 2(1+)yGa oT 
Lgpeey \0x? Peal 1-2» ‘Ox 


3 Ox Oy Ox? | 
¥G ( Oru +3 2(1+»)yGa oT 
1—2» \dx dy * dy? 1—2y * “Oy : 
k +t (0% O7u 
+ V2 (Ae _ 4 ( a BS 
ae V3 oy? Ox oy 04 


a 
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1.e. 
O*u Oru 0®u 2(1 +) yG 
7 x ie — a ogc 
B ox2 Oy? d OX Om wl = 2 Ox? | 
ome, Ot dy  -2(14+r)yGa oT ei) 
Ox dy Ox? +8 = A EN Ber) Oy? 
where 


x yG nyk —1 1G b ‘s 


dy 


B (az) +04 6 - a (2) +p =0 


or 


Sf a- Bt a EA. 


Real characteristics can only exist if at least one of the two values of (dy/dx)? 
is real and positive. The condition for this to be the case is that 


1+ p?— 93 q2— (1+ Bye 


> eee Tinks > OF 
Substitution from equations (36) and (30) shows 
Me tt BYE 
2B 


is always negative. Consequently there are no real characteristics at the point 
considered. But the point chosen was an arbitrary point in the plastic region. | 
It follows that the characteristics are imaginary at all points in the plastic 
region. Hence the equations for an elasto-plastic workhardening material are 
elliptic everywhere —it is well known that the characteristics are imaginary in 
the elastic region. Since the temperature satisfies LAPLACE’s equation V?7'=0, 
there are no discontinuities in 07 /0x and 07/0y in the interior of the body. 
Since the right-hand sides of equations (35) are continuous, the displacements, 
strains and stresses and all their derivatives are continuous in the plastic region. 


The Stress Function . 


In generalized plane strain the stress function yp is introduced by means of 


the equations 
Ory On _ Oy 
3 Onze = “Oy? ? Ony = — Ox oy and Oyy — Ox ° (37) 


aA ue Ne er 
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The conditions of compatability required to introduce these second order par- 
tial derivatives are identical to the equilibrium equations and are therefore 
satisfied. In an elastic material the stress function y satisfies 


Vy =0 (38) 


and this equation has proved suitable for the solution of many particular prob- 
lems. 

In an elasto-plastic non-hardening material, y will satisfy an equation which 
reduces to equation (38) when the plastic strain is zero. The equation will there- 
fore contain fourth or higher derivatives of y. But due to the hyperbolic nature 
of the equations for a non-hardening material in the plastic region, disconti- 
nuities can occur in the derivatives of stress, i.e. in the third order derivatives 
of y. It follows that the equation for y is not suited to the solution of particular 
problems since its fourth order derivates do not exist along certain lines. It is 
necessary to refer the equations to these lines, i.e. to the characteristics. 

In an elasto-plastic work-hardening material the derivatives of all orders 
of stress, and therefore of stress function, are continuous and, by analogy with 
elasticity, it is reasonable to suppose that solutions of the equation for the 


stress function will provide a convenient method of solving some particular 


problems. This equation will now be found. 
The strain compatibility equation is 


tr, Ae Oe 
4x Oy ~~ dye + ona 39) 
Since, from equations (11) and (37), 
oh ae a ifs il O*y F O2y Oy 
p= —{V%, ksin2 p 5 (= ml and eM Deal tae (40) 


substitution from equations (26) for ¢,,, ¢,, and ¢,, into equation (39) yields 


sae 1 V3 02 02 
Vy + (7 ut ae es sak Le Sicaames pr 2 ed 


V3 (/ & o?). 02 3 
+o | (ae _ iat) Sin2p + 2 ara cos2p} =O 
or 
1—~p V3 02 02 
( 2G + LiF. rey + VA f( 2 jaz) 8in2p + 259 7 cos2g}—0. (41) 


From the last two of equations (40), 


ory ory Oy 
2 tan2 57 OF een ye O. (42) 
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Equations (41) and (42) are two equations to determine the stress function yp. 
p can be eliminated between them but the resulting equation is too complicated 
to be of any use. Note that as 7 > co equation (41) reduces to the elastic 
equation V4y = 0. An elasto-plastic material with an infinite rate of work- 
hardening is indistinguishible from an elastic material with an infinite yield 
stress. 


Some Analytic Solutions 


In general equations (41) and (42) have to be solved numerically. However 
in three simple degenerate cases, A, B and C below, an analytic solution can 
be found and the corresponding practical problem solved. In a fourth case, 
D below, the stresses are functions of one variable only but the equation for p 
is non-linear. The complete solution for this problem has not yet been found. 

Case A. p = p(x): From equation (42) m = 2/4, and then, from equation (46), 
p satisfies d*yp/dx* = 0. 

Whence 


A B 
ioe 31 ay ee Oe D.. 


Only the first two terms in the expression for y give rise to non-zero stress com- 
ponents. They are 
Greg 5 Ug, HO Ope Ae +S. 


This is the stress system in the plastic regions of a bar bent in plane strain 
(Figure 2). Note that solution A, like solutions B and C, is a degenerate solution 


mirage) Saree 
Elastic y , 


Plastic 


Figure 2 
Plastic and elastic regions in pure bending. 


in that the equation for the stress function has become linear. Consequently 


_ the yield criterion is no longer identically satisfied and it is necessary to choose 


the constants so that the criterion is satisfied for initial yielding on the elasto- 


plastic boundary. The constants are different in the two plastic regions. This 


involves no contradiction because the two regions are not connected. 


| Vy +K (52 a5 — Fae et) 28 


Case B. yp = y(r): In plane polar coordinates equations (41) and (42) become 
1 0? 3 0 0? 


as ieetee 4 (43) 
ele. DOO on 5) 00829] Lay 
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and 
eal Oy LO OG VE Ory. ay Oy 1 Ogre 
2, tan2 P (; OF 00 7 v 4) t Ore y OV v2 062 0 , (44) 
where - 
1 3G 
i fle 45 
K alt oe) (45) 
and 
A= — aos (46) 


Putting y = (7) in equation (44), p = 2/4. Then from equation (43), w satisfies 
a 1 d \2 
Cr ee 


Whence 
y=A+ Blnv+C74DrInz, 
and the components of stress are 


Oy C= Sete Din cele aKa 


Cpe 99h Pp aor. 
0? B 
69 = pr = or + 2Dinr4+ 3 D420, 
Dye te WeOw aes 
Grp = ae (apap) = Os 


The constants B, C and D are not independent. The condition that the dis- 
placement should be single-valued imposes a homogeneous linear relation be- 


Figure 3 


Pressures and temperatures on a thick cylinder. 
tween them. This stress field arises in the expansion of a thick-walled right 


circular cylinder by internal pressure—the tube may be subject to tempera- 
ture gradients. The solution to this problem has been given by the author [1]. 


| 
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Case C.y = A. Equations (43) and (44) are satisfied with p = 0. The stress 


system is y 


er 


(47) 


This stress field arises in a thick-walled right circular cylinder subject to shear- 
ing stress on its surface. As far as the author is aware, this problem has not 
previously been solved explicitly and so the complete solution will now be 
given. 


\ 
N 


S PlastiGee 
eS 


Figure 4 
Shearing stresses on a thick cylinder. 


The stress field in the elastic region is also given by equation (47). If the 
torque per unit width is S, then 


an 
|r 2nr ~40=S or A 


0 


Whence 
Ora 4 ee” 


The radius c of the elasto-plastic interface is given by 079 = kp or 


Vs 


DM 


va = 


The current yield stress at a point, radial coordinate 7, in the plastic region is 


| { S c* ho 
4A k= LC MeO eat v2." 


vf 
Bc 


The equivalent plastic strain at such a point is, by equation (22); 
nS 
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In particular at the internal surface 7 = @, it is 

> eriee 

(6?) = (= 7 1) 
The parameter A is given by 


Wer kee Tae y3 (1 “), BETES, Co 


From equation (23), 


3 (ee 
sere tee en tae = ee 
ef = 25, = On, 6 = Fy (= 1), Cans 
From the generalized HooKE’s law 

Z 1 i Ge 


Chee seeee iad, 05 
&, = &e=9, {6- 26 
Therefore 


—- =— , 


k c2 3 (ce 
Le ae (4 =i}, ax<r<c 


lay GF 
2G. ee 


E,, = &9, =0, e+ = 


CAG e 


The strain-displacement equations in plane polar coordinates are 


Ou uU il Ov ee ieee Ov v 
Or yt 208 ail pie O60 tare Aerie OS 
Therefore 
hn VE) eye V3 he 

[- (ae psoas = ; tiny + 5 (cal eee 

u=0, v= aS See 
Rig C? Ro 
a Se = 
| NOTA ep, AI ae: 


where the constants of integration associated with the rigid body motion are 
chosen so that v is continuous and w= v = 0 at ry =c. The normal derivative of 
the transverse displacement is discontinuous across the elasto-plastic boundary. 


Figure 5 
Wedge with one edge under pressure. 


Case D. p=r* (6): In this case the stresses are functions of 0 only. This is 


the type of stress field one would expect in the plastic regions of a semi-infinite 
f 1 
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wedge of finite angle subject to uniform stress on its faces. It has been shown 
[8, 4] that when the wedge is subject to uniform increasing normal stress on 
one face, the plasticity commences simultaneously at both faces. As the stress 
increases further, one expects to have an elastic central region bounded by two 
plastic regions as shown. The components of stress in both regions are given by 


Se = KX +24, O9=24, O4=—-7'. 
In the elastic region y satisfies /4y = 0 or 
%=A+BO+4+ Ccos2 6+ Dsin2 6. 
In the plastic region y and @ satisfy 


0 " as 0 1 - O 
a5 (y"+4y7+K aa, sn2o+Zkcos2g=C, C=const (48) 


and 


W 


tan 2 YP = ae é (49) 


The stresses and displacements are continuous on the elasto-plastic boundary. 
The solution of equations (48) and (49) has not yet been found. 


Conclusions 


In this paper the author has considered elasto-plastic materials in gener- 


alized plane strain and has shown that the nature of the equations is different 


in the two cases of non-hardening and linear work-hardening. In the former 
case the characteristics and the equations of stress and velocity referred to them 


_ were found, and in the latter case the equation satisfied by the stress function. 


The procedure for solution of particular problems has yet to be formulated. 
The author thinks that the next stage in building up the workhardening 

theory is to look for analytic solutions of the type of case D above, e.g. try if 

a solution of the form y = y(xy) exists and then solve the non-linear ordinary 


differential equation for y, if necessary numerically. Dependent on the choice 
of boundary conditions, this solution will give a possible region with its stress 
and displacement fields. Next attempt to fix elastic regions adjacent to the 


plastic regions so that the final body and external force distribution bear some 


resemblance to a practical problem. In this way it may be possible to construct 


a 


f 


a few analytic solutions which will serve as useful guides to the first approxi- 


mation to be chosen for the solution of a given problem —in particular for the 


_ position of the elasto-plastic boundary. 
~ — Since the equations for y in both the elastic and plastic regions are fourth 


_ order elliptic, two conditions are required at all points on the boundary of the 


a's alee \ let 
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body and four relations between the variables in the elastic and plastic regions 
on the elasto-plastic interface. If it is assumed that the plastic region gradually 
expands as the external forces increase, then it is reasonable to suppose that 
the stresses on the elastic side of the interface are about to satisfy the yield 
criterion. Hence there are five conditions to be satisfied on the interface, 
namely, the continuity of the three components of stress and of the two com- 
ponents of displacement. Since only four of these conditions are needed for a 
solution in both regions for a given boundary, the author presumes, in absence 
of evidence to the contrary, that the fifth condition will only be satisfied if the 
interface has been correctly chosen. In practice in a numerical solution, it is 
unlikely that the interface will be correctly chosen initially. The magnitude and 
sign of the error in the fifth condition at each point on the interface should give 
an indication of what the next choice of interface should be. 

In problems solved by either elasto-plastic theory, it will often be possible 
to calculate the residual stress distribution when the external forces are re- 
moved. To do this, the stress distribution, obtained by applying to the body 
forces equal in magnitude but opposite in sign to the original forces and calcu- 
lated on the assumption that the stress-strain relation is elastic throughout, is 
added to the original stress distribution. The sum represents the residual stress 
distribution provided that the yield criterion is nowhere violated. 
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Zusammenfassung 


Unter der Voraussetzung, dass der hydrostatische Druck nicht gross ist, tritt 
die Hauptspannung normal zur Ebene des verallgemeinerten ebenen Verzerrungs- 
zustandes in der Fliessbedingung von Tresca nicht auf, Man kann zeigen, dass die 
Spannungs- und Geschwindigkeitsgleichungen fiir ein elastisch-ideal-plastisches 
Material hyperbolisch sind, wobei die beiden Charakteristikenscharen paarweise 
koinzident sind. Die auf die Charakteristiken bezogenen Gleichungen unterscheiden 
sich von denjenigen des starr-plastischen Materials nur durch die Gegenwart von 
vier Zusatztermen in jeder der Geiringer-Gleichungen. 
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Fiir ein Material mit Verfestigung werden nur kleine Spannungen beriicksichtigt, 
und es wird angenommen, dass die Gleichungen zwischen den Spannungen und den 
Gesamtdehnungen verwendet werden kinnen. Unter diesen Voraussetzungen wer- 
den die Gleichungen fiir die Spannungen und Verschiebungen elliptisch. Da alle 
Ableitungen stetig sind, existiert im plastischen Bereich eine Airysche Spannungs- 
funktion samt speziellen Ableitungen aller Ordnungen, und es kann gezeigt werden, 
dass sie einer nichtlinearen partiellen Differentialgleichung vierter Ordnung geniigt. 
Fiir einige einfache Probleme, bei denen die abhiangigen Veranderlichen nur Funk- 
tionen einer der beiden unabhangigen Variablen sind, werden die Lésungen an- 
gegeben. 

Der Einfluss von Temperaturschwankungen kann sowohl bei idealplastischen 
Materialien wie auch bei solchen mit Verfestigung mitberiicksichtigt werdeu. 


(Received: August 11, 1958.) 


Eine Methode zur Berechnung quasistationarer 
Temperaturzustande bei gesteuerten Warmeleitvorgangen 
in festen Korpern 


Von HAns JACKEL, Freiberg, Sachsen, Deutschland +) 


1. Vorbemerkungen 


In einem der letzten Hefte dieser Zeitschrift behandelte V. VoDICKA ein 
spezielles, gesteuertes, quasistationares Warmeleitproblem [1]?). Der hierbei 
verwendete Ansatz scheint im Falle beliebiger zeit- und ortsabhangiger Steuer- 
funktionen zur Lésung wenig geeignet. Deshalb soll nachfolgend ein Losungs- 
verfahren fiir ganz allgemeine Aufgaben dieser Art entwickelt werden. Die 
dazu notigen Gedankengange bauen auf einer Arbeit des Verfassers auf [2], in 

welcher mittels gewisser Einflussfunktionen das allgemeinste gesteuerte Warme- 

leitproblem gelést wurde. Aus den dort angefiihrten Ergebnissen lasst sich 
-ndmlich unter Vernachlassigung aller Gréssen, die mit wachsender Zeit gegen 
Null gehen, ein Theorem zur Bestimmung beliebiger quasistationarer Tempe- 

raturzustande ableiten. Die auf Grund der physikalischen Gegebenheiten an die 
- Steuerfunktionen zu stellenden Bedingungen ergeben sich hierbei ganz zwangs- 
-léufig. Bemerkt sei noch, dass aus [2] auch der Lésungsweg des von T.-S. CHow 
n dieser Zeitschrift erérterten Warmeleitproblems [3] entnommen werden 
kann. Als Beispiel wird schliesslich die von V. VopIcKA behandelte Aufgabe 
-gelést und die Aquivalenz beider Ergebnisse nachgewiesen. 


i ii i r kademie. 
u titut fiir angewandte Mathematik der Berga Ke" 
2 fe Ziffern in kes Klammern verweisen auf das Literaturverzeichnis, Seite 142 


PSAP SOY 
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2. Aufgabenstellung 


Wir betrachten einen beziiglich der Warmeleitung homogenen und iso- 
tropen festen Kérper K, dessen Oberflache / aus 7 analytischen Teilflachen 
F,, %,..., & besteht, die sich stetig aneinanderschliessen. Mit P bezeichnen 
wir die Punkte des dreidimensionalen Raumes, da wir uns bei den folgenden 
Untersuchungen auf kein spezielles Koordinatensystem beziehen wollen. Die 
Anfangstemperaturverteilung des Kérpers K sei durch T(P, 0) = /(P) gegeben. 
Die Temperaturen der Grenzschichten der die 7 Teilflachen umgebenden Medien 
sollen von der Zeit ¢ und vom Ort P abhangen. Diese sogenannten Steuer- 
funktionen bezeichnen wir mit ®,(P, t) (k =1, 2,...,7). 

Gefragt wird nun, unter welchen Bedingungen die Temperaturverteilung 
T(P,t) des Kérpers K im Laufe der Zeit einem quasistationaren Zustand 
T(P, t) zustrebt und wie dieser im Falle seiner Existenz berechnet werden 
kann. Dabei verstehen wir unter dem quasistationdren Temperaturzustand 
T(P,t) die Temperaturverteilung, die sich fiir grosse Werte der Zeit ¢ im Kérper 
K einstellt. 

Bemerkt sei noch, dass Probleme, bei denen im Inneren des K6rpers K 
zusatzlich Warme entwickelt wird, durch eine einfache Transformation auf die 
eben formulierte Aufgabe zuriickgefiihrt werden koénnen [2]. Wir gehen deshalb 
hier nicht weiter darauf ein. 


3. Lésung des Problems 


Um derartige Aufgaben zu lésen, kann man im wesentlichen zwei Wege 
gehen. 

Der erste basiert auf der physikalisch plausiblen Tatsache, dass die Anfangs- 
temperatur /(P) keinen Einfluss auf den quasistationaren Zustand hat und die 
Steuerfunktionen ®,(P, t) jedenfalls fiir alle in Betracht kommenden Werte 
der Zeit ¢ beschrankt sein miissen, wenn eine solche Temperaturverteilung 
T(P, t) iberhaupt existieren soll. Man wird so, nach der Fourierschen Theorie 
der Warmeleitung, zu dem Randwertproblem 

oT 


ape @AT =D{T(P. m= o, 


(1) 


Ar ee = 
(2 seach ie LEE, t)}n, =D(P, i) (hel, 2a) 
, 
gefihrt. Dabei ist wie iiblich A der Laplacesche Operator, die Aussere Norma- 
lenrichtung, aA die Temperaturleitzahl, A die Warmeleitzahl und «, die zur 
Teilflache F, gehérige Warmeiibergangszahl. Das Problem (1) versucht man 


- 


— 
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nun in jedem speziellen Fall durch einen direkten Ansatz fiir T(P, t) so zu 
lésen, dass alle Bedingungen erfiillt sind (vgl. zum Beispiel {1, 4]). 

Fiir die folgenden, ganz allgemeinen Betrachtungen wird es jedoch zweck- 
massig sein, den zweiten Weg zu wihlen. Diesem liegt die Lésung 7(P, t) des 
Randwertproblems 


D{T(P,H}=0, T(P,0) =f(P), 


2 
L{T(P, tke, =Dx(P, 2), Bie Qyy) Gase, Bos, ay #) - 
zugrunde, welche fiir alle Werte der Zeit 0 < ¢ < co den Warmeleitvorgang 
richtig beschreibt. Die quasistationiare Temperaturverteilung 7T(P, t) ergibt 
sich dann im Falle ihrer Existenz, wenn wir in dem Ausdruck fiir T(P, #) alle die 
Gréssen vernachlassigen, die mit wachsender Zeit ¢ beliebig klein werden. Auf 
diese Weise erhalt man neben der allgemeinen Lésungsmethode fiir die eingangs 
formulierte Aufgabe zwangslaufig auch die an die ®,(P, ¢) zu stellenden Be- 
dingungen. 
Die Lésung des Randwertproblems (2) ist bekannt [2]. Sie lautet 


T(P, t) = T,(P, ?) 
‘| ie : (3) 
2 G(P, 0, 2) +| Es GE(P, 7,1 —7)| de — HIP, 2) 


a 
0 


mit 
D{T(P,t)}=0, TWP, 0) =f(P), | i. 
L{T(P,t)}7=0, #>0, | 
L{H,(P, 1) hp, = —6;,P;,(P, 7) (5) 
f=1,2,...,7; k=1,2,0..,7; Oj,—1 fort =k, Ope 0 U0 tee), 
D{G#(P,t,1)}=0, GE(P, 1, 0) = H,{P, 2), & 


L{GEP, 1, f},=0, t>0 (F=1,2,...,7). 


In (5) und (6) ist t als Parameter aufzufassen. Ausserdem sei erwahnt, dass die 
G*(P, t,t —t) als Einflussfunktionen gedeutet werden konnen, wodurch eine 
physikalisch sinnvolle Interpretation des Ausdruckes (3) méglich istelD, 2h 
Bezeichnen wir schliesslich mit w,,(P) die aus (4) bzw. (6) resultierenden ortho- 


a * . . . 2 . was - kann 
normierten Eigenfunktionen, mit 67 die entsprechenden Eigenwerte, so 
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fiir (3) geschrieben werden 


es: 
w=1 


k=1 |p=l1 
(7) 
oe dB\ (x) 820-*\y_ HIP, 2) 
Ee eg? P —@a le —T me ; : 
+] 2% 5 u,,(P) e Tt w( 
Hierin ist 
A,= [ H(P) u,(P) do, B® (z) = [ HAP, t) u,(P) dv. (8) 
K K : 


Um die gleichmassige Konvergenz der in (7) auftretenden Reihen sicherzu- 
stellen, verlangen wir von f(P) und @,(P, t), dass die zweiten Ableitungen nach 
den Ortskoordinaten in K existieren und stetig sind [6]. Dariiber hinaus miissen 
die ©,(P, t) einmal stetig nach ¢ (0 <¢ < oo) differenzierbar sein. Diese letzte 
Voraussetzung garantiert ndmlich, wie man anhand der allgemeinen Lésungs- 
form des 3. Randwertproblems der Potentialtheorie [7] leicht verifiziert, die 
stetige Differenzierbarkeit der H,(P, ¢) nach ¢. Es gilt demzufolge die wichtige 
Relation 


Bhat atest AP) dv) <M (R=1,2....9 (9) 


fir 0 St St) < ov. Damit ist die Giiltigkeit von Gleichung (7) fir 0 <t< co) 
gesichert, vorausgesetzt natiirlich, dass sie fiir alle diese Werte der Zeit auch 
physikalisch sinnvoll ist. Man tiberzeugt sich leicht, dass die angegebene Lésung 
tatsachlich alle in (2) geforderten Eigenschaften hat. Von einer teilweisen Ab- 
schwachung der an die Funktionen /(P), ®,(P, ¢) gestellten Forderungen sehen 
wir ab, weil in den weitaus meisten praktischen Fallen diese an sich beliebigen 
Funktionen durch solche approximiert werden k6énnen, die den formulierten 
Bedingungen geniigen. 

Jetzt wenden wir uns der Bestimmung der quasistationaren Temperatur- 
verteilung T(P, é) zu. Ein Blick auf (7) zeigt, dass die erste Summe mit wach- 
sender Zeit ¢ beliebig klein wird. Daraus folgt aber sofort die Unabhangigkeit des _ 
quasistationaren Zustandes von der Anfangstemperaturverteilung /(P). Das- 
selbe gilt fiir die erste Summe in der geschweiften Klammer. Dies besagt nach 
der physikalischen Interpretation von (7) [2], dass die Einflussgréssen der 
Anfangstemperaturen ®,(P, 0) der umgebenden Medien ebenfalls keine Bedeu- 
tung fiir T(P, t) haben. Daher muss lediglich noch der Ausdruck 


2a : 
Ze- u,(P) eC 2H-9) dy — HAP, t) (10) 
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untersucht werden. Auf Grund unserer Voraussetzungen und der zusiitzlichen 
Forderung 


jim |®,(P,)} Sm (k=1,2,...,7) (11) 


ergibt sich unmittelbar, dass die H,(P, 2) und somit ebenfalls die Bw (c) fir 
t bzw. t > oo beschrankt bleiben. Dariiber hinaus gilt dann (9) auch fir T > 00.7 
Weil die Steuerfunktionen ®,(P, ¢) ohnehin eine gewisse Grésse (Schmelzpunkt, 
Entziindungstemperatur, absoluter Nullpunkt) nicht iiberschreiten diirfen, 
wenn die Fragestellung physikalisch sinnvoll sein soll, bedeutet die Bedingung 


(11) keine Einschrankung der Allgemeinheit. Fiir grosse Werte der Zeit ¢ wird 
nun 


‘. k 
> co dBM(z) 
a my P 
| i ate) ee alge |e ie as : 
Le PES p=l I 
Daraus folgt aber sofort die Existenz von 
7. oo a Sy aha | 
T(P, 1) = [Jim] 3; ss ~ u,(P) e-# Ht-) dy — H,(P, t) |. (12) 


Hierbei deutet das Symbol [ jim. | an, dass in dem nachfolgenden Ausdruck 


alle Glieder zu vernachlassigen sind, die mit wachsender Zeit ¢ beliebig klein 
werden und nicht etwa der Grenziibergang ¢ > oo schlechthin auszufiihren ist. 
Nunmehr kénnen wir das Ergebnis der eingangs formulierten Aufgabe all- 
gemein angeben. Es lautet: 
Sind die Steuerfunktionen ®,(P, ¢) zweimal stetig nach den Ortskoordinaten 


und einmal stetig nach der Zeit ¢ differenzierbar und bleiben sie dariiber hinaus 


fiir ¢ > co beschrankt, dann existiert ein quasistationdrer Temperaturzustand 


—iT(P, t). Um T(P, t) zu berechnen, bestimmt man zunachst die 7 von dem. 
} Parameter t abhangigen stationaren Temperaturverteilungen H;,(P, t) aus (5). 


Danach lést man die r vollhomogenen ebenfalls von dem Parameter t abhan- 
gigen Warmeleitaufgaben (6). Jetzt bildet man den Ausdruck (10), fithrt die 
Integration aus und vernachlassigt alle Glieder, die mit wachsender Zeit ¢ 


_beliebig klein werden. Das Resultat stellt dann die gesuchte Temperaturver- 
teilung T(P, ¢) dar. 


Abschliessend sei erwahnt, dass aus dem allgemeinen Ergebnis eine ganze 


_ Reihe spezielle Lésungen abgeleitet werden kénnen. Wir beschranken uns hier 
auf Steuerfunktionen der Form 


D,(P, t) = xu(P) Pell) - 


s In diesem Fall wird 
g GE(P, 1, 1-1) = galt) GP, #7) 


oS \ aaa) tC eka 
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gesetzt. Die Gleichungen (3), (5) und (6) gehen dann tiber in 


T(P, t) = T)(P, 2) 


; (13) 
- 3 nt G(P,) + [PE GP, t— 2) dr — pull) He Py 
k=1 6 
AH,(P)=0, L{H,(P)}n,= — Sin tel), | - 
Gis, 2g ial 2 eae 
D{G,(P,)}=0,  G,(P, 0) =H,(P), Me 
LIGUP Mee 0. fe Oq, Wale 
Fir T(P, 2) gilt jetzt somit 
et ol BES Da | [8 LAU 5) By? u,(P) e- @8it-9) de — g(t) H,(P) |. (16) 


Man erkennt, dass in diesem Fall die Randwertprobleme (14) und (15) nicht. 
mehr von dem Parameter t abhangen. 


4. Beispiel 


Als Beispiel wahlen wir das von V. Vopi¢KA auf andere Art behandelte 
Problem. Dieses lautet in unserer Bezeichnungsweise, wenn (r, @, 2) Zylinder- 
koordinaten sind, 


D{T} = 22 _ @ AT 


Pek oT ik igo Th, Qe Oo 
a or ee oe) = 0 


Lr, ©, 2, 9) = 17, @; 2), 


A arr: 
L{T}r, ou bse aya TN Site = ¥1(, 2) e*!, 
nor = 
UT},= (4-347) =neoe™ beso, =~ § 
OL | 
L{T}r,=( 7 Aise a) R= ge . 
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Man sieht sofort, dass die Steuerfunktionen die spezielle Form 
®,(P, t)= xP) v(t) 
haben und alle Voraussetzungen erfiillt sind, um die geschilderte allgemeine 
Methode anwenden zu kénnen. Wir bestimmen deshalb, da wir uns nur fiir den 


quasistationaren Zustand interessieren, nach (14) vorerst die H «(P) = A, (7, p, 2) 
(k =1, 2, 3). Mittels bekannter Lésungsmethoden finden wir 


SF =>" Ose [a cosn y + 5) sinn g] J, (7 k,v) [cosk, z+ y,sink,z], i= [en 


l=1 n=0 
mit 

2m s 
‘aay a cosn : 

a —2 A | / x1(@, 2) | eg | [cosk, z+ y,sink, 2] dg dz 

In ion = 

of | | ik Tq ib LR) + Juli rR)] (2(2-+ vb) Ars + (1— yf) sin2 hs] 

1 


H, =)" Sd” (a®, cosn y + 6°), sinn g] Tl 7) (Q® sinh k,, 2 +.Q® cosh hy, 2] 


mit 
R 2x 
“i - cosn p 
-2 f | al ) ¥ In(m 7) an a dy dp 
é 


0 


= 2 Qy 1275 a) O(2 
a Rn I(2- noe 2g R) + Tel P| [22 28 + 22) 9) 


128 oH x [a® cosn pm + 5°), sinn @] J,(R, 7) (Q® sinhk,, z — 2° cosh k,, z] 


‘m,n 
m=1n=0 
mit 
R 2a 
- CcOSn p 

3) —2 [ [rte P)* I nhm) eas a dy dp 

4in.n 0 0 
x 2 oe 1275, 20 Q (2 3 
b®,) Rn (+ — pepe) Tah P) + LE R)] [2 20) + Bea) 
m 


- Dabei bedeutet der Strich an dem jeweils letzten Summenzeichen hier und im 
folgenden immer, dass alle Koeffizienten mit dem Index Ls 0;1, m beliebig 
- durch zwei zu dividieren sind. Ferner ist unter J, die Ableitung der Bessel- 
~ Funktion /,, nach dem gesamten Argument zu verstehen. Alle in dieser Arbeit 
beziiglich der Bessel-Funktionen verwendeten Relationen findet man in der 
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Literatur [8]. Die k, bzw. k,, sind die Wurzeln der Eigenwertgleichungen 


A A 5 
—ksink s —cosks —khsink s —cosk s - a 
Ay bre Xs 0 J n(k ) es A RR 
d We R : 
UA ehar Seems UG oaee ene: Jn (hk R) “4 
Oy Og 
Schliesslich ist 
A : A : 
a, Pisink;s — cosk;s q, fu sink;s — cosk,s 
Yi 7 7 ? 
— k,cosk,s + sink, s —_k,cosk,s + sink;s 
Ao Xs 
Qo =|(* &,, sinhe, hk, 
mete ae my SD eS = COSI aS 
bere 253) 


Qe — | hm COSh Ry», S + sinhk,, s| 
Fir die G,(P, t) =G,(7, @, 2, t) (k =1, 2, 3) ergibt sich nunmehr nach (15) | 


G,= »; Ds: Doe [Aye cost p + BY? ,sinn | TABS r) 


x [cosk, z + yp, sink, z] et (Hi +m) 


> 


2% s 
eee cosn ; 
ae thr Jilhm®) | [rales 2) {0% 0 [oostys + sink 2] depde 
9 Oy &, 0 —s 
(1) ir n* is Fs 
Bima) Ra thh+Rs)|(1— ==) by R) + Jeb a) [2 (1+ yf) As + (1 — yf) sin2hy 
2 
R22 
ae er? = cosn 
AP, 7 [y, cosk; s — sink, olf ft ) ¥Tn(Rin) ee i dv dp : 
BY) Sa Be | n® 5, iP 7 
mm} Reeth Rb] [(1— 2) TR) + Teh R)) [2 (14+ vf) s+ (1 yf sina 
m ; 
R 2x 
i) ’ 4 Ri [y, cosh; s + sins) | [zs 7) Ta A) tee a dv dp | 
on 0 6 
th f 
B®) Dee n® je Be ] 
Pan) Rave hay (1-2 a) Fa R) + Teh Ry] [21+ wf) s+ (1 — yf) sim 2 
\ m 


: 
: 
| 
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Jetzt bilden wir entsprechend (16) den Ausdruck fiir T(P, t) =T(r, ps2, t). 
Es wird 


3 co 8606©0—‘“# 
T= Sve!) ¥) >; > iF a aN [A/?,,, cosnp +B , sinn g] 


l=1 m=1n=0 Wp +a? (hE + he) 
x Talhy r) [cosk, z+ y, sink, 2] — H,(7, q, z)}. 


Damit ist das Problem in aller Vollstandigkeit gelést. 
Zum Schluss zeigen wir, dass unser Ergebnis mit dem von VopI¢KA 
angegebenen tibereinstimmt. Nach den dortigen Ausfiihrungen ist 


T(r, 9, 2, 2) = alPs G85 b) a 
Wir beschranken uns hier darauf, die Giiltigkeit der Relation 


co ce co 
a ae 1p A™ cosn@a +B". sinn 
PR Ow Ti éo,+@ (A? + R2) =i eae Y l,m,” | 


l=1 m=1n=0 


(" 
x TulRyn r) (cosk, z+ y, sink, 2] — H,(r, 9, 2) 


fiir k =1 nachzuweisen. Unter Verwendung unserer Schreibweise ergibt sich 
fiir w) nach den Uberlegungen von VopIcéKA 


Ug = e'™! S* SUCH) cosn y + Di? sinn g] J,(¢ &7) [cosk, z + y, sink, z] 


1=1 n=0 
mit 
22 s 
cosn 
co | 2k; [ [use 4) id [cosk,z+ y,sink; 2] dp dz 
In das 
ee ee |Z ie Jali er) +J(é ex7)] (2 (1+ vf) Ars + (1 ~ of) sin 2 hs] 
1 
und 


on SEW. 


Weil die linke und rechte Seite der Gleichung (*) denselben Randbedingun- 
gen geniigt, braucht die Aquivalenz beider Seiten lediglich noch fiir das Zylin- 
derinnere nachgewiesen zu werden. Dort gilt aber fiir 6, = €; bzw. 8, = hk; die 
_Entwicklung 
kim Inlhm R) | fis By R) + Jnlt 6. R)| I nlFim 7) 


mat (B+ 3] \(8- e pe) Sil R) +Ji2Rn 2) 


k2 R2 


m 
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Fihren wir nun die Ausdrticke ftir up, Hy, af), 0), Aa Bia encre Df} 
in (*) ein, dann ergibt sich unter Beriicksichtigung dieser Entwicklung unmit- 
telbar die Richtigkeit der Relation fiir k =1. Die Beziehungen fiir Bi 2a 
beweist man ganz entsprechend. 

Wir bemerken noch, dass sich unser Ergebnis sehr leicht in Real- und 
Imaginarteil aufspalten lasst. Die in H,(r, y, z) auftretenden Grossen 


J n(t hi”) 


Fi hy Tid AR) + Juli) 
z 


sind nadmlich fiir » =0,1,2,... immer reell, weil dies fiir J,,(¢k,7) und 
i Jonsi(t Rk, 7) gilt und ausserdem die Beziehung 


Ti by) = > (Trae Ba”) — Inarlé a) 


besteht. Beachtet man dies, so bereitet die Aufspaltung keinerlei Schwierig- 
keiten. 

Die Lésung des von CHow behandelten Problems [3], die hier nicht wei- 
ter erértert werden soll, folgt schliesslich aus (13). Dazu muss lediglich noch 
T,(r, y, 2, t) aus (4) bestimmt werden, worauf wir jedoch verzichten wollen. 
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Summary 


Ina recent paper VopicKa solves a special problem of heat conduction with 
local and time-dependent boundary conditions. As the method does not seem very 
suitable for solving problems of this kind, the solution for the most general case is _ 


given. As an example the problem treated by Vopréka is solved, and the equality 
of both results is shown. 


(Eingegangen: 20. Mai 1958.) 
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Zur Aerodynamik von Hagelkérnern 


Von Roranp List, Weissfluhjoch-Davos!) 


1. Einleitung 


Dem zunehmenden Interesse an Vorgingen der Hagelentstehung ent- 
sprechend, sind in den vergangenen Jahren verschiedene Theorien iiber das 
atmospharische Wachstum von Eiskérnern publiziert worden (siehe zum Bei- 
spiel [1-4]?)). Die dabei gemachten Annahmen iiber Kornformen, Luftwider- 
standskoeffizienten und Fallgeschwindigkeiten basieren alle auf Verhaltnissen 
fiir glatte Kugeln. Daraus ergeben sich jedoch Resultate, die normalerweise 
von den tatsachlichen Umstianden erheblich abweichen. — Um solchen Wachs- 
tumstheorien beziiglich der Aerodynamik von Hagelkérnern eine bessere Grund- 
lage zu geben, seien in der Folge Messungen an natiirlichen Eisschlossen und 
Modellkérpern beschrieben. 


2. Messapparatur und Messmethodik 


Die Messungen der Luftkrafte auf Modellkérper oder Hagelkérner wurden 
im bereits friiher erwahnten Windkanal [5] in einem Kaltelaboratorium des 
Eidgenéssischen Institutes fiir Schnee- und Lawinenforschung vorgenommen. 
Der Aufbau des Versuchskanals ist aus Figur 1 ersichtlich und zeigt, dass die 
angesaugte Luft mit einem Ventilator, angetrieben durch einen Zweiphasen- 
Kurzschlussankermotor von 1,1 PS Leistung, in einen Leitstern und ein Leit- 
gitter getrieben wird. Zwei sich folgende Venturidiisen sorgen fiir die nétige 
Beschleunigung und Geschwindigkeitsverteilung im eigentlichen Messteil. Die 
maximal erreichbare Geschwindigkeit betragt rund 25 m/s. Als weiteres Charak- 
_teristikum des Kanals sei erwahnt, dass bei einer Reynoldsschen Zahl 
- Re=5,1-104 der Widerstandsbeiwert c,, einer glatten Kugel 0,45 + 0,02 betragt. 
_ (Das kritische Gebiet mit dem Abfall des Widerstandskoeffizienten einer Kugel 
und dem zu c,, = 0,3 zugehérigen Re kann infolge ungeniigender Ventilator- 
-leistung nicht ermittelt werden.) Das zu untersuchende Objekt wird auf ein 
Stahlstangelchen mit einem Durchmesser von 3 mm starr befestigt; dieses 
_ wiederum kann in ein Messingrohr eingeschraubt werden. Die Aufhangung des- 
_selben erfolgt gelenkig, so dass eine Verschiebung des auszumessenden Korpers 


1) ERidgendssisches Institut fiir Schnee- und Lawinenforschung. 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 158. 
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auf einem Kreisbogen in Richtung des Luftstromes mdglich ist. Da das Messing- 
rohr in sich verdrehbar ist, konnen im Abstand von 22,5° verschiedene Exposi- 
tionen des Objektes um die gewahlte Drehachse eingestellt werden. Zur Mes- 
sung der Luftkrafte wird eine Federwaage eingesetzt, die unter Beriicksichti- 
gung des Hebelgesetzes den gewtinschten Luftwiderstand ergibt. Die Waage 
wird jeweils so verschoben, dass die Aufhangevorrichtung lotrecht 1 cm vor der 


y 


al 
Ll 


Schematische Darstellung des Versuchskanals: 1 Ventilator; 2 Leitgitter; 3 Leitstern; 4 Unter- 
suchungsobjekt; 5 Haltestangelchen; 6 Messingrohr; 7 Hebelgelenk; 8 Umlenkrolle; 9 Feder- 
waage; 10 Stativ. 


Figur 1 


Kanaloffnung steht. —- Als Korrektur ist jeweils der Luftwiderstand des von der 
Strémung betroffenen Haltestangelchens (ohne den im Objekt eingeschlossenen 
Teil) in Abzug zu bringen. . 

Diese Aufhangung hat zur Folge, dass ein auf das auszumessende Korn 
ausgetibtes Drehmoment um eine Achse parallel dem Gelenk der Aufhangung 
mitgemessen wird. Infolge der Lange des Hebelarmes (85 cm) und den maximal 
beobachteten Drehmomenten von der Gréssenordnung 1 g* cm ist jedoch’ 
dieser Effekt bedeutungslos. 

Die zu untersuchenden Hagelkérner werden auf dem Haltestangelchen auf- 
gefroren, nachdem dessen Spitze auf Rotglut erhitzt in das Objekt hinein- 
gefiihrt wurde. Allfallig entstehende Vertiefungen in der Eisoberflache kénnen 
mit zusdtzlichem Wasser auf die urspriingliche Form gebracht werden. — Holz- 
modelle werden auf die Halterung aufgesteckt. 

Bei der Fixierung der K6rner ist jeweils darauf zu achten, dass die Achse der 
Aufhangung mit einer Achse des das Objekt umhiillenden kleinsten Ellipsoides 
zusammenfallt. Die extremsten Verhaltnisse werden damit am ehesten erfasst. 
Wahrend bei Hagelkérnern normalerweise alle drei «Hauptachsen» auszu- 
messen sind, geniigt es, bei den rotationssymmetrischen Graupelmodellen nur 
eine Achse senkrecht zur Symmetrieachse zu betrachten. 
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Die Bestimmung des Widerstandskoeffizienten Cy erfolgt nach dem Gesetz: 
Weg 2 F. (1) 


Der Luftwiderstand W wird mit der Federwaage gemessen, die Geschwindig- 
keit v., mit einem Prandtl-Rohr. Die zu jedem c,, gehdrende Projektionsflache 
F des Objektes in eine Ebene senkrecht zur Bewegungsrichtung wird ermittelt 
itber eine Photographie eines in der Windkanalachse aufgestellten Photoappa- 
rates. Die Luftdichte 9, ergibt sich aus Lufttemperatur und Luftdruck. 

Da es sich zeigte, dass der Widerstandskoeffizient im Geschwindigkeits- 
bereich von 14 bis 20 m/s nicht variiert, wurde dessen Grésse durch Mittelung 
zweier Messungen an den Grenzen dieses Bereiches berechnet; viel gréssere 
Geschwindigkeiten konnten mit dem Windkanal nicht erreicht werden, bei 
kleineren war die Genauigkeit und Empfindlichkeit des MeBsystems un- 
gentigend. So darf angenommen werden, dass der absolute Fehler der bestimm- 
ten Widerstandskoeffizienten + 5° nicht iiberschreitet. 

Bei den Versuchen musste die Zeitdauer der einzelnen Messreihen méglichst 
beschrankt werden, um ein allzu grosses Verdunsten der Hagelkérner im Wind- 
kanal einzudammen. Zum gleichen Zwecke wurden die Objekte mit Polyvenyl- 
formal, dem Schaferschen Replica-Material [6], eingehiillt (Schichtdicke etwa 
0,02 mm), so dass jeweils nur eine Gewichtsverminderung von ungefahr 4%, 
eintrat. 


3. Entwicklungsstadien eines grossen Hagelkornes 


In der Publikation Kennzeichen atmosphdarischer Eispartikeln, Teil 1 [7], 
wurde bereits auf den Hauptweg der Entstehung der Hagelkérner hingewiesen 
mit den Phasen Eisbildungskern-Schneekristall-Reifgraupel-(Wassertropfen)- 
Frostgraupel-Hagelkorn. Eine prinzipielle aerodynamische Betrachtung der 
einzelnen Stadien ist jedoch nur sinnvoll, wenn nicht zu viele Hauptformen 
parallel nebeneinander bestehen. Der grossen Variation im Aussehen entspre- 
chend, wird also nicht auf Eisbildungskerne und Schneekristalle eingegangen. 
Reifgraupeln wiederum sind erst interessant von der Grésse an, wo der Habitus 
des urspriinglichen Schneeteilchens nicht mehr deutlich sichtbar ist, das heisst 
von einer Grésse von rund 3 mm an. Aus verschiedenen Messungen und Be- 
obachtungen (zum Beispiel [7]) wissen wir, dass die Mehrzahl der grossen Reif- 
graupeln einen kugelsektorartigen Habitus aufweist. So wurden denn zwei un- 
- polierte Holzmodelle nach der oben beschriebenen Methode im Windkanal aus- 
gemessen. Die Resultate sind dargestellt in Figur 2 A und B, wobei Haupt- 
 spannt, Widerstandskoeffizient und Fallgeschwindigkeit als Radiusvektoren 
_ in der Luftstromrichtung vom Nullpunkt weg aufgetragen sind. Bild A gibt die 
 Resultate wieder fiir ein Modellkorn mit einem Kegeléffnungswinkel von 70°, 
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0 05 1 0 Q2 O04 f o4 ; ae 
Vq [m/s] Cy Flom] 
Figur 2 


Hauptspannte (F), Widerstandskoeffizienten (c,,) und relative Fallgeschwindigkeiten (v5) von 
Modellgraupeln in Funktion der Fall-Lage. 


wahrend dieser fiir Bild B 90° betragt. Generell wurden beziiglich der Korn- 
achse symmetrischeWerte gemessen, Abweichungen traten nur auf, wo infolge 
von Kornschwingungen eine genaue Messung nicht méglich war. Der Vergleich 
zeigt, dass der Widerstand in der Hauptachsenrichtung abhangig ist vom «Ab- 
reisswinkel» der Strémung an der ausgepragten Kérperkante. 
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Aus der Gesamtheit der Widerstandskoeffizienten interessiert uns jedoch 
nur der Widerstandskoeffizient fiir stabile Lage bei freiem Fall. L. v. Bucu [8] 
beobachtete bereits 1814, dass Reifgraupeln mit der Spitze nach hinten fallen. 
Darlegungen von ARENBERG [9] scheinen jedoch auch ein Fallen mit entgegen- 
gesetzter Lage zu erméglichen — eine Ansicht, die ein Experiment jedoch nicht 
bestatigte. Bei Graupelfallen im Freien wurden namlich Reifpartikeln in einer 
Glasschale mit einer Wasserschicht von 5 cm Dicke aufgefangen; der gebremste 
Fall verifizierte in jedem Fall die erste Behauptung. 

Fiir den freien Fall von Reifgraupeln kommen also Widerstandskoeffizienten 
von rund 1,0 in Betracht. Dies entspricht auch den Resultaten von NAKAYA 
[10], die umgerechnet analoge Werte ergeben. Wenn wir nun, wie verschiedene 
Autoren — zum Beispiel [3] — in Unkenntnis der Hauptkornformen Kugeln be- 
trachten, so ergibt sich fiir diese bei gleichem Volumen und gleicher Dichte, 
wie sie eine kegelige Reifgraupel vom Typ A aufweise, eine um 87% héhere Fall- 
geschwindigkeit (!). Da die Modelle eine geringere Rauhigkeit aufweisen als 
natiirliche Reifpartikel, so erh6ht sich dieser Wert in Wirklichkeit noch. Die 
Zuverlassigkeit der theoretischen Berechnungen muss also fiiglich bezweifelt 
werden! 

Sobald die auf eine Reifgraupel auftreffenden kleinen Wassertropfen nicht 
mehr relativ schnell gefrieren, hat das Wasser Zeit, in die Kapillarsysteme des 
lockeren Reifgeriistes einzudringen. Es entsteht eine Frostgraupel, deren Cha- 
rakter nach Modell C (Figur 2) gegen die Form D wechselt und zur Kugel wird. 
Figur 3a zeigt, wie bei allen kugeligen Reifgraupeln bei der Umwandlung die 


Reifgraupel Reifgraupel 
t% x 
ge ee ‘ ee ai =. 
Scohmelzen teilweises Schmeizen _Jangsames Gefrieren Schme/zen teilweises Schmelzen  langsames Gefrieren 
Wasserfropfen Wassertropfen 
Gefrieren Gefrieren Gefrisren 
eae 
Kareispartikel Typ IT Frostgraupel Klareispartikel TypI Klareispartikel IypI frostgraupel  Klareispartikel Tyo I 
a b 
Figur 3 


Entwicklungstendenzen von Reifgraupeln: a ausgehend vom kugeligen Korn, b ausgehend vom 
kegeligen Korn. 


Form erhalten bleibt, wahrenddem bei den meisten kegeligen dieTendenz zur 
Kugel vorherrscht (Figur 3). Die stark hervorgehobenen Linien zeigen die 


a 


L 


haufigsten Varianten. 
Eine Deutung der stabilen Fallagen ist ftir die Frostgraupelstadien nicht 
mehr allgemein verbindlich, da auch die Dichteverteilung stark variieren kann. 
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Es ist nur die Aussage méglich, dass bei glatter, das heisst klarer Oberflache der 
Widerstandskoeffizient gegen den der Kugel strebt. — Fiir die Fortsetzung des 
Wachstums lehrt die Erfahrung, dass in der Natur keine glatten kugeligen 
Hagelkérner mit Durchmessern iiber 10 bis 15 mm auftreten. Der Einfluss der 
Oberflachenspannung muss dem Einfluss der Aerodynamik weichen, der das 
fiir klares Eiswachstum notwendige freie Wasser an die Kornperipherie treibt. 
Dadurch entstehen Rotationsellipsoide, deren Widerstandskoeffizienten bereits 
durch PRANDTL [11] gemessen wurden (siehe Figur 11). 


Figur 4 Figur 5 
Korn K 57.A, Diinnschnitt D 1 in Korn 57.A, Diinnschnitt D 1 bei 
durchscheinendem Licht (1:0,67). polarisiertem Licht (1:0,67). 


Ein Wechsel zu schnellem Eiswachstum ergibt eine neue Formvariante: 
dreiachsige, ellipsoidische Typen (siehe auch [12]). Dieser Hauptweg der Ent- 
wicklung grosser Hagelkérner aus rund 50 Sammelproben schwerer Hagelwetter 
kann besonders schén verfolgt werden aus den Figuren 4 und 5, die einen 
0,3 mm dicken Diinnschnitt von Korn 57.A in durchscheinendem und polari- 
siertem Lichte wiedergeben. Die Ursprungspartikel kugeliger Form zeigt im 
polarisierten Licht eine kegelige Struktur — sie ist also aus einer kegeligen Reif- 
graupel entstanden. Dass nach der Umwandlung in eine Frostgraupel diese 
immer noch mit dem Ort der urspriinglichen Spitze nach hinten fiel, geht aus 
dem angesetzten Blasenkranz hervor, der an der Kegelbasis wuchs und nicht 
das ganze Teilchen umschliesst. Die dadurch erfolgte «Auswuchtung» des Pam 
tikels liess die Fallage andern, sehr wahrscheinlich sogar instabil machen Die 
nachste Phase ergibt sich durch die Vergrésserung zum Rotationsellipsaas An 
dieses angelagert ist eine Zone geringerer Dichte, die urspriinglich nur aus Pro- 
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tuberanzen bestand, wie sich aus den dazwischenliegenden Klareiskandlen fol- 
gern lasst. Die Kérperumhiillende kann als beidseitig in der kleinsten Achse 
eingedriicktes Ellipsoid bezeichnet werden. Dieser Wachstumsphase tiberlagerte 
sich jedoch 6rtlich eine zweite, was zur Auffiillung des urspriinglich lockeren 
Eisgeriistes fiihrte. 

In jedem dieser Zwischenstadien kann jedoch die symmetrische Form in 
eine unregelmassige oder kugelartige iibergehen, je nach den Wachstumsbedin- 
gungen. Das heisst also, dass der Habitus der Hagelkérner nicht zufallig ist. 
Die resultierende Kornform ist charakteristisch fiir die Wachstumsbedingungen, 
die die Hagelschlosse im Laufe threr Entstehung durchlaufen hat. Diese Form- 
beziehung ist jedoch auch gréssenabhangig. 

Die Luftkrafte, die auf diese ellipsoidischen und unregelmassigen Haupt- 
formen wirken, sollen nun im nachsten Abschnitt dargestellt werden. 


4. Luftkrafte auf Hagelkérner 


Es wurden folgende vier Haupttypen aus verschiedenen schweren Hagel- 
wettern herausgenommen und ausgemessen: 


Korn K 57.24. Form: beidseitig in Richtung der kleinsten Achse eingedriick- 
tes Ellipsoid (siehe auch [7]). 
Oberflache: relativ glatt. 
Gewicht: 24,3 g*. 
Spez. Gewicht: 0,908 g*/cm?. 
Hauptachsenverhaltnis (ohne Beriicksichtigung der Einschnti- 
rung): 1,0:0,79:0,44. 

Korn K 57.25. Form: einseitig in Richtung der kleinsten Achse eingedriicktes 
Ellipsoid. 
Oberflache: glatte Protuberanzen. 
Gewicht: 32,4 g*. 
Spez. Gewicht: 0,86 g*/cm?. 
Hauptachsenverhaltnis (ohne Beriicksichtigung der Einschnii- 
rung): 1,0:0,78:0,55. 


- Korn K 57.26. Form: kugelartig. 


Oberflache: rauh mit Protuberanzen. 
Gewicht: 18,59 g*. 
Spez. Gewicht: 0,87 g*/cm*. 


- Korn K 57.27. Form: dreiachsiges Ellipsoid. 


VAY RASA SY 
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Oberflache: rauh. 

Gewicht: 23,3 g*. 

Spez. Gewicht: 0,88 g*/cm’?. 
Hauptachsenverhiltnis 1020, 77.0;,64: 
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Die Umrisse der vier Korner, jedes aus den drei Hauptachsenrichtungen be- 
trachtet, sind in Figur 6 dargestellt. 


BDpeP 


5OO@ 
BBD S 


K 57.27 K 57.26 TGS E25 K 57.24 


Figur 6 
Umrisse von Hagelkérnern, beobachtet aus den drei Hauptachsenrichtungen: 
Erste Reihe: Kornachse /, senkrecht zur Bildebene; 
zweite Reihe: Kornachse /, senkrecht zur Bildebene; 
dritte Reihe: Kornachse J, senkrecht zur Bildebene. 
(G25) 


Die Messergebnisse sind den Figuren 7-10 zu entnehmen, wobei wiederum zu 
beachten ist, dass Hauptspannt, Widerstandskoeffizient und zugehérige extra- 
polierte Faligeschwindigkeit je als Radiusvektoren in der Luftstromrichtung 
vom Nullpunkt weg aufgetragen sind. Konsequent wird dabei die langste Korn- 
achse als Drehachse /,, die kiirzeste als Drehachse /, bezeichnet. Zur Kontrolle 
der Messungen kénnen jeweils auch die Werte, die sich bei gleicher Anstrém- 
richtung (in Richtung der Hauptachsen!) einstellen, verglichen werden. Von 
den sechs méglichen Kombinationen sind 3 in den entsprechenden Figuren 7 
bis 10 eingezeichnet. Die Abweichungen bedeuten jedoch nicht unbedingt Mess- 
fehler, sie rithren eher davon her, dass die Hauptachsen nicht ganz senkrecht 
zueinander getroffen wurden. Dass eine ungenaue Festlegung der Nullstellung 
auch gréssere Differenzen ergeben kann, geht speziell aus Figur 8, Drehachse Jz, 
hervor. (Normalerweise wurde die Nullstellung ermittelt durch folgendes Vor- 
gehen: In das Versuchsobjekt wird ein rechtsdrehendes Koordinatensystem 
gelegt, die Austrittspunkte der positiven Achsen werden markiert. Das Auf- 
gefrieren erfolgt jeweils in der Achsenrichtung durch einen so bezeichneten 
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Figur 7 
Korn K 57.24, raumliche Verteilung der Hauptspannte, Widerstandskoeffizienten und Fallge- 
schwindigkeiten. 


Punkt. Die Ausgangsstellung 0° ergibt sich durch Drehen der lingsten «freien» 
Achse in Richtung des Luftstromes). 

Widerstandskoeffizienten und Fallgeschwindigkeiten der verschiedenen 
Korner bei ausgezeichneten Anstrémrichtungen seien in Tabelle 1 zusammen- 
gestellt (siehe auch Figur 6). Ein Vergleich der von vorn und von hinten an- 
gestr6mten Objekte zeigt, dass eine Verdrehung um 180° speziell bei unregel- 
massigen Hagelschlossen sehr unterschiedliche Widerstandskoeffizienten und 
Fallgeschwindigkeiten ergeben kann! — Ein Einfluss der verschiedenartigen 
_ Einschniirungen ist jedoch nicht zu ermitteln. 

Die beobachteten betrachtlichen Differenzen der Resultate je nach der 
_ Anstrémrichtung lassen sofort die Frage nach der stabilen Lage bei freiem Fall 
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Figur 8 


Korn K 57.25, raumliche Verteilung der Hauptspannte, Widerstandskoeffizienten und 
Fallgeschwindigkeiten. 


auftreten. Da die Beantwortung auf Grund der Messung der Luftkrafte nicht ge- 
macht werden kann und Messungen iitber Drehmomente fehlen, wurde der Fall- 
versuch im Wassertank herbeigezogen. Dieser ergab unter Beriicksichtigung 
der Reynoldsschen Ahnlichkeit fiir den freien Fall in Luft, dass die ellipsoidi- 
schen Objekte stets mit der kleinen Achse lotrecht fallen! Ein- oder beidseitige 
Einbuchtungen in Richtung dieser Achse beeinflussen diese Tatsache nicht. 
Ein Einfluss der Achsenverhaltnisse (speziell J,//;) auf den gréssten Wider- 
standskoeffizienten zeichnet sich nicht ab, da die Kérner infolge ihres ver- 
schiedenen Oberflachencharakters nicht ohne weiteres vergleichbar sind. Es 
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Figur 9 
Korn Kk 57.26, raumliche Verteilung der Hauptspannte, Widerstandskoeffizienten und 
Fallgeschwindigkeiten. 
Tabelle 1 
Widerstandskoeffizienten und Fallgeschwindigkeiten von Hagelkérnern 
bei ausgezeichneten Anstromrichtungen 
Korn K 57.24 57.25 57.26 57.27 
Cw O, Cy Ug Gr Ug Cy Vg 
mit Strom 0/36) 4078s) 10,267 9°50,0) 05545 27518) 031 4157 


K 
Reece BiniSencaty ah) 0) S5iabs4 bs5 lO: 2Gen.49,6,1N0,53 ):26)9 0 S6uraae 


Pe, rome 0,33 38,8 | 0,30 42,1 | 0,62 24,8] 0,42 34,4 
Ornacdse "2 gegen Strom | 0,35 37,7 | 0,33 39,4 | 0,64 24,3 | 0,40 35,2 


mit Strom On73 IS /OU0165; 4 23,28 00,700" 214 O62 2357, 


ee epier Serre. eh 065s 21 Li O70" 22 AbOS7 24,2 10,69" 422.5 
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Figur 10 


Korn K 57.27, raumliche Verteilung der Hauptspannte, Widerstandskoeffizienten und 
Fallgeschwindigkeiten. 


ist aber anzunehmen, dass sich dieser Beiwert mit zunehmendem Scheiben- 
charakter bis 1,15 und mehr vergréssert [11]. (Die Grésse des Widerstands- 
beiwertes wird bei unregelmassigen Kérnern auch beeinflusst durch das Ver- 
haltnis Querschnittumfang zu Querschnittflache, wodurch die c,-Werte fir 
Scheiben noch iiberschritten werden kénnen.) . 

Zur Errechnung der Korngeschwindigkeit bei freiem Fall (auf Institutshéhe 
von 2665 m ti. M.) wurde angenommen, dass der Widerstandskoeffizient bis 
zu Reynoldsschen Zahlen von 2-105 konstant bleibt. Diese Extrapolation ist 
jedoch nur nétig fiir Fallagen, die in der Natur nicht stabil sind, die fiir uns also 
nur vergleichenden Charakter haben. Die Verhiltnisse fiir stabile Fallagen 
konnten mit dem Windkanal gerade noch erfasst werden und sind damit ver- 
bindlich. Da bei schweren Hagelfallen mit tiber 3 cm Partikelgrésse rund 
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Tabelle 2 
Vergleich der effektiven mit der theoretischen Fallgeschwindigheit von Hagelkéynern 


| ‘ 
Ab h 
’Hagelkorn °Rugel Beas a 
m/s m/s % 
Hagelkorn K 57.24 ... 19,9 29,6 49 
Hagelkorn K 57.25 ... 22,4 30,4 36 
Hagelkorn) K 57:27... ¢ . 22,5 29,1 29 


80% der Hagelkérner ellipsoidische Formen aufweisen, soll die Tabelle 2 die 
erheblichen Abweichungen der tatsdchlichen Verhiltnisse von den theoretisch 
berechneten aufzeigen. Die Rechnung wurde fiir ein kugeliges «Ideal-korn» 
gleichen Volumens und gleicher Dichte bei einem c, von 0,5 ausgefiihrt. 
Man muss sich bei der Betrachtung dieser Zahlen klar sein, dass es sich nicht 
um extreme Abweichungen handelt. Im Gegenteil, es sind klare Hinweise im 
Aufbau der Hagelkérner vorhanden, dass in friiheren Wachstumsstadien (glei- 
cher Grésse!) die Oberflache rauher war und die Protuberanzen wesentlich deut- 
licher zutage traten (siehe Figur 4). Dass die Rauhigkeit der Oberflache und 
vereinzeltes Auftreten von Protuberanzen Widerstandskoeffizienten und 
Querschnittsflache vergréssern, zeigen speziell auch die Messungen an dem 
ausserlich als kugelig zu taxierenden Korn K 57.26. In jeder Fallage ist der 
Hauptspannt grésser (im Maximum bis 30%) als bei einer Kugel gleichen Vo- 
lumens und Gewichtes, ebenso der Widerstandskoeffizient (bis 40%). 

Damit sind wir in der Lage, Variationsgebiete der Widerstandskoeffizienten 
der Haupttypen der Hagelkérnerim Laufeihres Wachstums anzugeben (Figur11), 
allerdings unter der Beschrankung auf maximale Durchmesser von 6 cm. Die 
obere Begrenzung dieser Gebiete soll fiir rauhe, mit Protuberanzen besetzte 
Schlossen gelten, wahrenddem die untere eher glatten Hagelkérnern entspricht. 

Die Zusammenhange zwischen Geschwindigkeit und mittlerem Haupt- 
durchmesser der Hagelk6rner sind dargestellt in Figur 12. Die wahrscheinlich- 
sten Zonen sind durch entsprechende Schraffierung hervorgehoben, die zu- 
grunde gelegten Werte fiir Dichten sind fiir Reifgraupeln 0,5 g/cm’, fiir Frost- 
graupeln und Hagelkoérner 0,8 g/cm3. 

Neben der Fallgeschwindigkeit ist auch die Wachstumsrate dm/dt von gros- 
ser Bedeutung. Diese ist, unter Annahme gleicher Einfangswahrscheinlichkeit, 
proportional dem Produkt aus Geschwindigkeit und Querschnittsflache senk- 
recht zur Fallrichtung. So ermittelte Werte sind fiir reale Formtypen und 
Widerstandskoeffizienten im Vergleich zu glatten Kugeln gleicher Masse und 
Dichte dargestellt in Tabelle 3. ' 

Ein Vergleich zeigt, dass in diesen Fallen Hagelkorner praktisch iiberein- 
stimmende Wachstumsraten wie ideale Kugeln gleicher Masse und Dichte auf- 
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Figur 11 


Variationsgebiete fiir Widerstandskoeffizienten (c,,) von Graupeln und Hagelkérnern in Funktion 

der Reynoldsschen Zahl Re: Kurve I, Beiwerte fiir glatte Kugeln, gemessen von PRANDTL phils 

Kurve II, Beiwerte fiir glatte Kugeln, gemessen von BiruAm und Rerr [1]; Kurve III, Beiwerte 

fiir abgeflachte, glatte Ellipsoide, gemessen von PRaNnpTL [11]; Kurve IV, Beiwerte fir Ellipsoide 
(Tabelle 1); Kurve V, Beiwerte fiir Scheiben, gemessen von PRANDTr [11]. 
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Fallgeschwindigkeiten von Graupeln und Hagelkérnern (2665 m ti. M.): 1 Kugel, @ = 0,8 g/cm’, 

ea Oioeme Kugel, Q = 0,8 g/cm, c,, = 0,7; 3 Kugel, 9 = 0,8 g/cm’, Cy = 1,0; 4 Kugel, g = 0,5 g/cem?, 

Cy = 1,0; 5 Kugelsektor, Kegeloffnungswinkel 90°, @ = 0,8 g/em?, Cy = 0,8; 6 Kugelsektor, Kegel- 

hay GAP ane f= 0,5 g/cm’, Cy = 0,8; 7 Kugelsektor, Kegeléffnungswinkel 70°, 9 = 0,5 g/em®, 

Cy = 1,0; 8 dreiachsiges Ellipsoid, @ = 0,8 g/cm3, Cy = 0,7; 9 dreiachsiges Ellipsoid, @ = 0,8 g/cm, 
Cy = 1,4. 
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Tabelle 3 


Vergleich der zur Wachstumsrate proportionalen Produkte (v-F) dveier Hagelkirner mit 
idealen Kugeln gleicher Masse und Dichte 


eee ee ee 


(v-F) (v-F) (v-F) 
Hagelkorn Kugel 
m3/s m?/s % 
Hagelkorn K 57.24... 3,30 °10>2 sisilks}en (Obst — 3,6 
Hagelkorn K 57.25 . . . | 4,08-10-2 <All FON +1,2 


Hagelkorn K 57.27... 2;92510-4 a, LOO = —5,8 
|e al eee | 


weisen, sofern die Luft gleich charakterisiert werden kann beziiglich ihres 
Feuchtegehaltes; nur sind wesentlich kleinere Aufwinde nétig, um eine Schlosse 
schwebend zu halten. Umgekehrt bedeutet dies, dass bei einem gegebenen 
vertikal aufsteigenden Luftstrom 1,5- bis 2,5mal schwerere ellipsoidische Kér- 
ner an Ort gehalten werden k6énnen, als es die Rechnung fiir Kugeln ergibt! 
Entsprechende Daten fiir reale «kugelige» Hagelkérner diirften einen Faktor 
von 1,2 bis 2,0 ergeben, je nach Rauhigkeit und Protuberanzen. 

Die natiirlichen Bedingungen aerodynamischer Natur erweisen sich damit zur 
Bildung grosser und schwerer Hagelkorner als wesentlich giinstiger als die Schluss- 
folgerungen aus theoretischen Uberlegungen. 

_ Eine Darstellung aerodynamischer Verhaltnisse bei der Hagelbildung ware 
unvollstandig, wenn nicht auch auf mégliche Rotationen der Eispartikeln hin- 
gewiesen wiirde. In der Tat sind auch untriigliche Anzeichen im Aufbau ein- 
zelner Schlossen zu finden, die auf eine Rotation hinweisen (Figur 19, {12)). 
Da aber der diesbeziigliche Einfluss auf den Widerstandskoeffizienten nicht 
speziell gross sein diirfte, sei auf diese meist um die Fallachse (!) erfolgende 
Wirbelung nicht weiter eingegangen. 

Damit ware also von seiten der Aerodynamik her die Basis gegeben, Be- 
rechnungen tiber das Hagelkornwachstum neu anzusetzen. Es wird aber vor- 
laufig darauf verzichtet, weil das tatsachliche thermische Verhalten der Hagel- 
kérner bei ihrem freien Fall durch die Atmosphare nicht bekannt ist. Zur 
Verfiigung stehende theoretische Uberlegungen basieren ebenfalls wieder auf 
Berechnungen an glatten Kugeln, was noch gréssere Fehlresultate als bei der 
Erfassung aerodynamischer Verhaltnisse ergeben diirfte. Auch ist zu beachten, 
dass die von SCHUMANN [2] gemachten und von Mason [3] und Lupram [4] 
iibernommenen oder weitergefiihrten Ansdtze tiber maximale Anlagerung von 
Eis in Funktion der héchst méglichen Abfiihrung der Gefrierwarme den tat- 
sichlichen Vorgangen nicht entsprechen. Die Natur kennt diese Wachstums- 
beschrankung nicht, da das nicht gefrierende Wasser grosstenteils im Eisgeriist 
eingebaut wird! Ebenso fehlen geniigende Experimente iiber die Uberlagerun- 
_ gen einzelner Wachstumsphasen, die zur Schalenbildung fithren kann. 
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5. Zur Frage des gréssten Hagelkornes 


Brruam und RELF [1] beschaftigten sich 1937 eingehend mit der Frage des 
grdssten spharischen Hagelkornes und glaubten, eine Limitierung des Durch- 
messers statuieren zu kénnen auf Grund des plétzlichen Absinkens des Wider- 
standskoeffizienten fiir Kugeln bei hdheren Reynoldsschen Zahlen. Ihre Uber- 
legungen wurden ins allgemeine Lehrgut der Meteorologie aufgenommen [3], 
ohne dass folgende zwei Punkte, die zu schweren Bedenken Anlass geben, be- 
achtet wurden: 

a) BirHaAm und RE LF fiihrten die Messungen an Kugeln, die an einem Flug- 
zeug angehangt waren, aus. Damit wollten sie jegliche Turbulenzen, die in 
einem Windkanal die Resultate beeintrachtigen, vermeiden. Vergleiche mit 
Messungen von PRANDTL [11] zeigen denn auch, dass dadurch eine Verschie- 
bung des Ablésepunktes oder, als Resultat davon, ein Absinken des Wider- 
standskoeffizienten von Reynoldsschen Zahlen von Re = 2,2-10° nach 
Re = 3,2-10° verschoben werden konnte (Figur 11). Da aber schon durch die 
Anwesenheit vieler Hagelkérner in einem Aufwind einer Gewitterwolke dessen 
Turbulenzgrad erhéht wird, ist es nicht gerade sinngemass, Messungen in spe- 
ziell turbulenzfreier Luft fiir die Hagelbildung herbeizuziehen. 


b) Die Messungen von BiLHAm und RELF wurden an glatten Kugeln aus- — 


gefiihrt. Glatte kugelige Eiskérner sind aber in der Natur héchstens mit Durch- 
messern von 15 mm zu finden. Gréssere Hagelpartikeln weisen bei ihrem Wachs- 
tum wohl 6fters diese Form auf, ihre Oberflache ist jedoch stark aufgerauht bis 
protuberanzenartig. Diese Oberflachenstruktur hat aber zur Folge, dass auch 
bei kleinen Re eine turbulente Grenzschicht vorhanden ist. Damit sind wir zum 
vornherein in einem tiberkritischen Gebiet der Reynoldsschen Zahl — der 
Widerstandskoeffizient wird sich nicht stark verandern, eine aerodynamisch 
bedingte Beschrankung der Hagelkorngréssen wird damit illusorisch. 

Fir die wertvollen Ratschlage zur Interpretation der Ergebnisse bin ich 
vor allem Herrn Prof. Dr. J. ACKERET vom Aerodynamischen Institut der ETH 
zu Dank verpflichtet. Danken méchte ich weiter Herrn Dr. M. DE QUERVAIN 
fiir die Férderung dieser Arbeit und den Mitarbeitern Herrn W. Putre und 
Herrn P. Rass fiir ihre Mithilfe bei den Messungen und Auswertungen. Die 


Untersuchungen wurden im Rahmen des Forschungsprogrammes der Eidge- 


ndssischen Kommission zum Studium der Hagelbildung und der Hagelabwehr 
an der Forschungsstelle Weissfluhjoch ausgefiihrt; zur Finanzierung standen 
teilweise Kredite des Schweizerischen Nationalfonds zur Verfiigung. 
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Summary 


All existing hail theories in mathematical form proceed from the assumption 
that hail stones can be represented as smooth spheres. Observation of the most 
common types of graupels and hail stones shows, however, that a sphenical shape 
is found only in a minority of cases, while the surface is generally more or less 
rough. In order to establish the importance of these obvious deviations from 

_ the theoretical assumptions, hail stones and various models of soft hail (graupel) 
were subjected to wind resistance tests in a wind tunnel. The apparatus and 
method used can be seen in Figure 7; the results are contained in Figures 2, 7, 8, 
9, 10. The values given for the coefficient of wind resistance, for the actual 
section facing into the air stream, and for the expected speed of fall have in each 
case been recorded along the direction of the current outwards from the centre 
of the figure. The measurements were made with respect to a rotation round the 
main axes of the hailstones. 

Experiments in a water tank further showed that all these measured hail stones 
fell in one direction, in which the speed of fall was minimal. 

If we compare this speed of fall with that of ideal, smooth spheres having the 
same mass and density, differences result of between 30 and 50% (Table 2). And 
as the measured hail stones have still a relatively smooth surface, it is reasonable 
to expect that even greater deviations can in fact occur. 

The extent of variation in coefficients of resistance for natural hail stones 
expressed as functions of the Reynolds number has been entered in Figure 77. — 
Various connections between the particle radius (the mean radius of the section 
facing the direction of flow) and the particle’s speed of fall are represented in 
Figure 12 for a number of different hail-stone shapes. In the case of these rela- 

tionships, too, there is considerable divergence from the usual calculations. 

It is also an interesting fact that the products of the speed of fall v and the 
cross-section F are approximately as great for observed hail stones as for spheres 
of equal mass and density (Table 3). Based on the assumption of constant col- 
lection efficiency these products are proportional to the rate of growth. 

Everything goes to show that natural aerodynamic conditions essentially 

favour the formation of larger and heavier hail stones than is allowed for by 
theoretical considerations. — It is also not admissible to postulate at the critical 
Reynolds number a maximum hail-stone diameter on the basis of the decline in 

the coefficient of resistance of a smooth sphere, since equivalent hail stones can 
be treated in the optimum case as rough spheres (with changed character of the 
boundary layer). 


' (Eingegangen: 26. September 1958.) 
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Theory of Oscillation Type Viscometers V: 
Disk Oscillating Between Fixed Plates’) 


By Gorpon Frank NEWELL”), Providence, R. IS USA 


Abstract 


A disk oscillating between two fixed plates is considered as an instrument 
for absolute measurements of viscosity. The existing theories relating the vis- 
cosity to the decrement of oscillation are improved by calculating the effects of 
the edge. This is done by assuming that the separation between the plates is 
small compared with both the radius of the disk and the boundary layer 
thickness. A comparison is made with the experimental data of KEstTin and 
PirarczyK for which the present theory is estimated to be correct to 0:1%. 


1. Introduction 


A disk oscillating in a fluid between two fixed parallel plates has been used 
by several investigators as an instrument for measuring viscosities. The status 
of the experimental and theoretical problems has recently been discussed by 
KeEsTIn and WANG [1]’) and by KEstin and Moszynsk1 [2]. As compared with 
a disk oscillating in an ‘infinite’ fluid described in parts III and IV ({3, 4]4), the 
presence of two stationary plates in the manner shown in Figure 1 tends to 
increase the viscous drag on the disk. This is desirable from an experimental 
point of view in cases where the infinite fluid would give a decrement that is 
too small to measure accurately. | 

The present paper will deal only with the theoretical aspects of the problem 
and will employ most of the same notation and mathematical formulation 
previously described in parts I [5], IT [6], III [3] and IV [4]. The only difference 
between the present problem and those discussed in the previous parts is a 


*) This research was supported by the United States Air Force through the Air Force Office 
of Scientific Research of the Air Research and Development Command, under contract No. 
AF 18 (600) 1548. Reproduction in whole or in part is permitted for any purpose of the United 


States Government. Part of the work was done while the author was being supported by a grant 
from the Alfred P. Sloan Foundation. 


2) Brown University. 
3) Numbers in brackets refer to References, page 173. 
4) Equation numbers and other references to earlier parts of this series will be denoted by the 


prefix I, 1, If or IV. For the most part, however, the present paper is dependent only on part I 
which also contains a table of notation. | 
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difference in geometry. As compared with the infinite fluid case of parts III 
and IV, the only new feature is the additional boundary condition that the 
fluid velocity must vanish on the plates instead of ‘at infinity’. 

The main deficiency of present theories is that they do not properly take 
into account the influence of the edge. Whereas KEsTIN and WANG have made 
a semi-empirical study of this aspect and based their considerations upon the 
experimental data of KesTIN and PILARczyxk [7], the present analysis will be 
directed toward obtaining the edge correction from first principles. This will 
permit one to use the viscometer for absolute measurements of viscosity 
instead of just relative measurements. . 

Briefly the problem is to solve equation (I, 13), namely 

0? w 3 dw 0?w 


a 


eee LOR Rec Oe ty OG (1) 


subject to the boundary conditions 


mE, %, Ss) = 1 on the disk A, 


w(é, 7, s) =0 on the plates 
and evaluate the integral, equation (I, 16), 


Aw 
D(s) =o 0° 5 [> doc. (3) 
Hl on 


The angular displacement «(z) of the disk consists of a transient plus a decaying 


oscillation of the type, equation (I, 18) 


a(t) = A e~4* cos(w t+ y) | (4) 
in which 


S=+(—A)o (5) 


are the roots of the equation (I, 40) 


(S + Ay)? +1+ D(S) =0. (6) 


The reader is referred to part I for a detailed discussion of notation and of 


“the origin of these equations but particularly relevant here is the fact that all 


lengths are measured in units of the boundary layer thickness 6 = (¥/@p) 


1/2. 
’ 


_ y is the (unknown) kinematic viscosity and «, the natural frequency of oscilla- 


A 


/ 
- 


tion. In the same units, we let & and 7p denote the radius and half thickness of 


the disk respectively as in parts III and IV. 


_ ZAMP X/11 


wou 4 
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To our previous notation, we must add two new lengths, namely the two 
distances between the plates and disk which we denote by f, and Bz (see 
Figure 1). We will assume, however, that B, and f, are nearly equal, 


beh at a 
ALL, 
s | w-0 z 


Geometry of a disk oscillating between stationary plates. 


In addition we assume that 
20+ Bit Bo<éSo (8) 


(the separation between plates is small compared with the radius of the disk) 
and 


27 + Bit Bol (9) 


(the separation between plates is small compared with the boundary layer 
thickness). We will not find it necessary, however, to make any further restric- 
tions on the absolute magnitude of &). 

Equations (7) and (8) are restrictions on geometry alone whereas equation (9) 
is in effect a restriction on 0. If, in contrast, with equation (9), we should have 
By ~ B, > 1, the plates would be located so far from the disk that they would 
have practically no influence on the fluid motion and the whole problem would 
reduce essentially to that of the infinite fluid considered in parts III and IV. 
If the plates are to serve some useful purpose that could not be achieved with 
the simpler physical arrangement without the plates, then certainly the plates 
should be placed at least within a distance 6 from the disk and preferably so 
that B, ~ B, < 1. If m) and f, are of comparable size then this condition is 
equivalent to equation (9). 

We have chosen to measure lengths in units of 6 in order to maintain a 
notation consistent with that of the previous parts, even though there is no 
particular advantage to this choice in the present problem. The boundary 
conditions, equation (2), along with equation (9), will force w to decrease from 
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1 to 0 in a distance small compared with 6. The physical significance of 6 as a 
‘boundary layer thickness’ is therefore essentially lost and we would expect 
the fluid motion (i. e. w) to be nearly independent of this, at present, rather 
artificial unit of length. 

The mathematical interpretation of this can be seen by noting what would 
happen to equation (1) if we should measure distances in units of (2+ Pi + Be) 


instead of 6, i.e. let & = (2 9+ Bi + By)-2& and v’ = (2 
; No + Bi + Bs) 7- 
Equation (1) would then become : ’ at 


02w 3 02 w 02w 


Og’? ar BOF + one 8 (2 79 + By + Be)? w. (10) 


Since we are interested in the solution of this for s ~ + 3, the right-hand side 
of equation (10) is small and, in some sense at least, w should be given approxi- 
mately by the solution of equation (10) or (1) with s = 0. We expect then that 
a power series expansion of w in powers of s should converge rapidly for | s | 
of order 1. 

Thus far we have not taken advantage of equation (8) which we now use 
to remove the last unpleasant feature of equation (1), the term 3 &-! dw/0é. 
We know that over most of the disk, not too near the edge, both dw/0& and 
0?w/0&? will be small. The neglect of both of these terms in equation (1) is, in 
fact, the basis of most previous treatments of this problem. ‘Not too near the 
edge’ could only mean within a distance of order (2 4) + 6, + 62) and here we 
expect 0?w/0n’?, 0?w/dé’? and Ow/0E’ in equation (10) are of order 1 but 
3(é')-! dw/0E’ is of order (€’)-1 ~ (2 ny + By + Bo) /So < 1. 

The smallness of this term corresponds physically to the fact that to an 
observer situated this near the edge, the edge appears to be nearly straight 
and the oscillations appear as linear rather than rotational oscillations. We 
anticipate, therefore, that an asymptotic expansion of w in powers of &)* will 
converge rapidly. 

The key to the solution of this problem is the above observation that to a 
first approximation w satisfies LAPLACE’s equation. This approximation to w 
is valid everywhere including the region near the edge. It will contribute to 
D(s) a first correction for the edge and will in this sense represent a second 
approximation to D(s) as compared with earlier theories. 


2. First Approximation 


Our goal is to obtain eventually an expansion for D(s) in powers of s and 
€,1. The contribution to D(s) from the edge will be smaller than that from the 
flat surfaces of the disk by a factor of order (2 4) + 8; + B2)/So oF smaller so 


we first present a brief review of the theory neglecting the edge effect in order 
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to obtain the terms in D(s) of lowest order in € * and certain qualitative 


properties of the disk. 
Neglecting the edge effect is equivalent to neglecting the radial dependence 


of w. To obtain all terms of lowest order in 51, we consider the approximation 


sinh (7) + B, — ||) s*? for é <2oy ee 


sinh pf, st? 
3 ey it i 1/2 
ea = mee ie IE for.£-<aeon oy] see Ay 
0 fore Se 


which is a solution of equation (1) except for the discontinuity at € = €. The 
lowest order in both €51 and s is obtained from the simpler approximation, 


ea i for§< £9, n> 0, | 
il 
Wy = “ea ln for <f,4 <0, * 
2 
0 . for E> &. | 


The expression for D(s) resulting from equation (11) is 
1 a — 
Dy(s) = &4 @ 6° I> 882 = [coth(A,)/s) + coth(f2//s)] , (13) 


a formula which has been derived and used previously [1]. We shall expand 
equation (13) in powers of s and write the result in the form 


Dy(s) = Dols) {1 + (3) 8 Br Bs — (ou) 8? Bt + 62) B 
+ (sag) G+ PB 4 
peri Dols) = Miro) s, (14a) 
p= ae. (14b) 
M =2x(€ 6)* I-1(8 6)-* we? (14c) 


D,(s) is the lowest approximation to D(s) and could be obtained directly 
from equation (12). It has been separated into the product of three factors M, 
(»@) and s the first of which depends only upon the suspension system but not 
the fluid (since the quantities 6 6 and & 6 are lengths measured in the laboratory 
system of units), while v 9 depends upon the fluid but not the suspension system. 
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Equation (6) has two complex conjugate roots which determine A and w, 
equation (5). It also has infinitely many real negative roots which determine 
the transient part of «(t). One can show, however, that the transient decays ina 
time short (of order 6-') compared with the natural period of oscillation and 
the transient can be eliminated from the observation of the main oscillation 
by waiting a short time for it to decay. 

Equation (6), being a complex equation, is equivalent to two real equations 
and, therefore can, in principle, be used to compute two unknowns. It is a 
practical experimental necessity, however, that A be small (less than 1/10 at 
least) and by solving equation (6)-with A < 1, one finds that 


l1-—-w <A <i. 


This means that the change in w due to the fluid is very small and in present 
applications, it has not been possible to measure the frequency change 1 — w 
with sufficient accuracy to be of much use. 

The imaginary part of equation (6) gives the real equation 


2@ (4 — Ay) = ImD (iw — Ao) , (15) 
which with D(s) ~ D,(s) yields as a first approximation 
2(4—Ay)) ~Mre. (15a) 


One may use this equation to compute v from measured values of M, @ and 
(A — A,). The higher approximations to equation (15) will, however, also 
contain w and one has the option of using either an observed value of w or a 
value computed from the real part of equation (6). Although the latter would, 
in most cases, be the more accurate value of w (otherwise one would be able 
to use the measured value of w — 1 to compute something else, perhaps @, 
more accurately than it is already known), the difference is usually insignificant 
as compared with other errors and one should use whichever value of @ is more 
convenient to determine. 

We now turn our attention to the problem of finding a more accurate 
expression for D(s) which will include the contribution from the edge and 
give a more satisfactory formula than equation (15a). 


3. Edge Effects | 


To find the next approximation to D(s), we say that w is approximately 
equal to the solution of LAPLACE’s equation for the geometry of Figure 1. To 
find this w, we simplify the problem further by taking 8, = 6, = f. The method 
to be used here could also be used for B, + f, but the edge contribution to 
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D(s) will be relatively small (of order f/&) and the error caused by a small 
difference in p,, By is expected to be smaller yet by a factor of order (8, — f)?/B?. 


Yl KN/, WO , eo 
Y} ~L+i(Tg*B) l bot! No B) Fy ‘ 


Figure 2 


The geometry in the complex z-plane showing boundary conditions on a polygon. 


We now modify Figure 1 slightly so as to correspond to Figure 2. We choose 
a coordinate system with its origin on the walls of the disk and consider 


w= Eg—é& (16) 
and 7 as coordinates in a complex plane 
2 re (17) 


Boundaries at x = €& and x = — L have been added to show boundary con- 
ditions on a polygon in the z-plane. 

The proposed scheme of expanding w in powers of &)1 and s will, at best, 
determine w to some finite power of >’. We anticipate that the influence of the 
edge decays exponentially and for € = 0 (x = &,) is proportional to the negative 
exponential of a number of order &/8 > 1. Such quantities are smaller than 
any finite power of £/& and can therefore be neglected. The problem, in question, 
therefore is to find w for the limiting case in which the boundaries of Figure 2 
at x = €) and x = — L all recede to infinity. 

We now determine w by conformal mapping or more specifically by means 
of a Schwarz-Christoffel transformation. One can easily determine a transfor- 
mation which maps the polygon of Figure 2 into a half plane with the boundary 
of the polygon mapping onto the real line. It can be chosen so that the boundary 
w = 0 maps into the negative half of the real line and the boundary w= 1 
maps into the positive half of the real line. It is convenient to make still a 
second mapping, a logarithmic transformation, which transforms this upper 
half plane into a horizontal strip. The pair of transformations then maps the 
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polygon of Figure 2 into the strip of Figure 3 with the numbered points of 
Figure 2 mapping into the corresponding numbered points of Figure 3. If we 


iv 
(ML ] 


| 1/2 w-0 


Figure 3 
Geometry in W-plane resulting from a conformal mapping of Figure 2. Numbered points are 
images of corresponding points of Figure 2. 


denote the coordinates of Figure 3 by 


W=U+iV (18) 
and define 
och fe. 
y = cos (1 + B ) (19) 


then the mapping in question is given by the equation 


~ 2th oe | sini W . _,/tanhw 
aoe {= Sin ( ar ) + coshy sin ( Gane ) \. (20) 


One must be careful, however, to use only the analytic continuation of this 
formula along paths lying entirely within the horizontal strip of Figure 3 of 
that branch of the multi-valued mapping for which z = 0 when W = 0. 

The purpose of the conformal mapping is to make use of the fact that 
LAPLACE’s equation is invariant to such a transformation, i. e., w(2(W)) con- 
sidered as a function of W also satisfies LAPLACE’s equation. The solution in 
the W-space is trivial, however; 


24W 2 
w(2(W)) = Re(1+ —T—) =1- (21) 
We need not find w as an explicit function of z because equations (20) and 
(21) already give a convenient parametric representation of the solution. For 
any value of W, equation (20) defines a point z(W) and equation (21) gives the 
value of w at that point. The point-by-point evaluation of w is not, however, 
_ of primary concern. We are mainly interested in evaluating D(s), which can be 
_ done directly in terms of w as a function of W. 
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Before we do this, it is interesting to investigate a few qualitative properties - 
of dw/dn in the z-space since this measures, in effect, the drag per unit area of 
the fluid on the disk and it appears in the integrand of equation (3). 

The present approximate solution for w gives 


ow 2 0V 2 dW cosh U ¥ 
On & OF in dz (cosh2y — cosh?U)!? 


along the end of the disk, | U | < y and 
Ow 2 NOV eee, cosh U (22) 


On 2 On a dz  £ (cosh?U — cosh®y)1? 


along the top or bottom of the disk; | U| > y. This along with equation (20) 
defines a parametric representation for 0w/dn as a function of z. Near the 
corner x = 0, 7 = 4) (U = y), 0w/0n becomes infinite. As a function of U, it 
has a — 1/2 power singularity but one can show from equation (20) that as a 
function of x for 7 = yy) and 0 < x <€ Np, 


Ow 1 —1/3 
ea E 37 p* x tanh | 
thus 0w/07 has a — 1/3 power singularity as a function of x. 
Of more importance is the behavior for U > y or equivalently, x > p. 
From equation (22) we see that 


Ow iy) 1 _ _ cosh?y 
an B " 2 cosh? U = | 
and from equation (20) that 
U~ an x + [coshy log (coshy + 1) — (coshy — 1) logsinhy] + ---. (23) 


Together these two equations show that 0w/07 is equal to B-! plus a term that 
decreases as exp (— x x/B). The term $-! corresponds to that which obtains 
from Wy, equation (12), while the exponential term confirms our earlier 
prediction that the edge effect would decay exponentially. 

The surface integral in equation (3) can be separated into two terms, one 
from the walls of the disk and the other from the top and bottom. By using 
the cylindrical symmetry and the reflection symmetry with respect to the 
plane 7 = 0, equation (3) can be reduced to the single integrals 


5 i &o No 
Die ees [aes +8 [dy at (24) 


0 
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We are interested here only in the excess drag caused by the edge. It is con- 
venient therefore to subtract the part that arises from the constant value 
0w/0n = B-* since its contribution toD(s) has already been included in the 
analysis of the previous section without some of the simplifying assumptions 
being used here. If we denote by D,(s) the contribution to D(s) from the edge, 
write the integrals in terms of x and expand & = &} (1 — x/&,)3 in powers of x, 
we obtain from equation (24) 


D,(s) -- ao | [as sree. e a + fe rae 
ee (25) 
SR, VD l ° 
+] dx = Ere Mes x al nee (5 | 


The third integral is of order 6/& and we shall discard it for now but 
consider it later in a higher approximation. The first two integrals form a 
single line-integral along the curve w = 1 of Figure 2 and they can be translated 
into a corresponding line-integral along V = 0 of Figure 3. By using the fact 
that W(z) is an analytic mapping, one obtains the form 


4 B D4(s) 
Eo 


U(&) F ae | 
D,(s) = / ja é 
0 


Vivo JB 


in which U(é,) is the value of U for which x(U) = &). By using equation (21) 
we obtain 


Da(s) = 4  &* Dols) [= Ul) — 8 & 


This, in turn, can be evaluated to within an exponentially small error from 
equation (23). We finally obtain, 


D,(s) = 8 B Dj(s) 271 £5 ' [coshy log (cosh y + 1) — (coshy — 1) logsinhy]. (26) 


D,(s) is to be added to D,(s) or D,(s) to give the term of order ree the 
expansion of D(s)/Do(s) in powers of s and 1. The quantity in brackets is a 
function only of y, equation (19) which, in turn, depends only upon the geo- 


metry of the suspension system. For B, = f2 = f, equation (26) gives exactly 


the value for the term of order £51! s° in the expansion mentioned above. The 


fact that we have used a solution of LAPLACE’s equation instead of equation (1) 


excludes the possibility of also finding the smaller terms of order £57, 5 


ako 


 etc., directly from the above results. 


ee AY 
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4. Higher Approximations 


Although it is very difficult to obtain a more accurate approximation to w 
from which to calculate still higher-order terms in the expansion of D(s), one 
can make more efficient use of the solution at hand by employing the varia- 
tional form for D(s) introduced in part III, section 4. If w is only an approximate 
solution of equation (1), then the formula 


ed's [f.. Ow 0?w 3 Ow Ow 
Dis) = 25 [eae e—fface olor tate 3 sw] (27) 
A V 


(in which V is the volume of the fluid and dr the element of volume) is an 
extremal and therefore insensitive to small errors in w. One should be able 
correctly to obtain from this at least the two terms of order aoe and of order 
s &,1 for 6, = B2 using only the approximate w of the last section. 

The first term of equation (27) is equivalent to equation (3) or (24) but we 
previously neglected the third integral in equation (25) and we should now 
include it with the terms of order &>*. Since our w is a solution of LAPLACE’s 
equation, the second derivative terms in equation (27) cancel but the first 
derivative term will also contribute something of order >? while the last term — 
of equation (27) will contribute something of order s. 

The integral containing 0w/0& can be evaluated to the desired accuracy by 
the following sequence of steps 


[face w3e Be ~- 308 [fan ag 
i 


N+ B 


No 
=62§ jf wl + [ay let ~ 20 & (3 mo + B) - 
0 No 


If we denote by D,(s) the contributions of order >”, then the above term 
combines with the corresponding term from equation (25) to give 


Ds(s) = 2 Dy(s) (£-)° 


(28) | 


(+3) +6/4e§ [ae 


n=No “ a ; 


We were not able to evaluate the integral in equation (28) explicitly 
although the integrand for integration with respect to U instead of x can be 
written explicitly and the upper limit of integration can be taken as oo. It 
is known, however, that the integral is a function of Mo/B only and that it is 
always negative, in contrast with the first term of equation (28). For any given 
geometry 7/f is fixed and the value of the integral is simply a constant number 


a 


Spaeth} 
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that can be calculated by numerical integration. In a specific example (with 
HolB ~ 1/2) a numerical integration showed that the contribution from this 
integral was only about one tenth aslarge as the contribution from the first 
term of equation (26). Although the integral was considered to be negligible in 
this instance, it is not obvious that this should always be the case. 


The only remaining significant contribution from equation (27) is the last 
term 


@ o> st [— [ax &2 w? , (29) 


This term includes all terms proportional to s? including those of relative magni- 
tude & and €,1. If we write w? = w? + (w? — w2), then the contribution from 
wy gives the part of order £} and the contribution from w? — w? gives the part 
of order €,'. The former part can be evaluated explicitly but it gives nothing 
new because this has already been evaluated by other means and appears as 
the second term in equation (14). 

The integral with w* — ws could not be evaluated explicitly but again this 
can be reduced to an integral that depends only upon 7/8. For any given 
geometry it could be evaluated numerically, if necessary. For the intended 
applications, however, this term is quite small. Compared with D,(s), it is of 
order (8/&)) 6? with both £/& < 1 and # < 1. Furthermore, it is proportional 
to s? and its contribution to Im D(i w — Aw) in equation (15) will be of order 
(B/&)) B? A smaller than Im Dy with A < 1 also. 

We will neglect this last contribution and take as our final approximation 
to D(s), the sum of D,(s), equation (14); D,(s), equation (26); and D,(s), 
equation (28), 


Dis) ~Mvos{C+ > 5B,Bo— op 8 B+ BB+ cay B+ ADB +} 
) 


(30 
with 


=1+ 8271 £ [cosh y log (coshy + 1) — (coshy — 1) logsinhy] 
0 


~ 1 
+2(£) iva sofarg (p (G+ 50) af 


Note that this differs from D,(s) only in that the first term of equation (14), 
the number 1, is replaced by C. C is a function of B/& and y only, however, and 
_ for any fixed geometry is just a new constant. 


(30a) 
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We could obtain as many terms as we like in the expansion of D,(s) in 
powers of s although we already have included more than would usually be 
necessary. Disregarding this as a source of error we see that equation (30) 
admits fractional errors in Im D (i w — Aw) and in calculated values of » which 
are of order 


(ZY. £ (4). £ pA, and fp. (31) 


The first error arises from terminating the series for C, the second from taking 
, = By in Dy, the third from neglecting the term from equation (29) and the 
fourth from the edge correction proportional to s°. 

The actual calculation of » from known values of all the other physical 
quantities can be done very easily and certainly with no greater difficulty with 
the above value of C than for C = 1. D(s), equation (30), is substituted into 
equation (15) and the latter is solved for y by a method of successive approxi- 
mations. Since all but the first term of equation (30) are relatively small, they 
can be approximated by using a preliminary estimate of » to determine 6 and 
thus B, and f,. The only other place that y appears in D(s) is as a numerical 
factor which one now calculates. If necessary, one can recalculate 6, and fy, 
to obtain new estimates of the small terms in equation (30) and repeat the 
calculation of ». 

In the experimental arrangement used by KeESTIN and PILARCZYK, 
Piléy = 00282, Ba/£, = 0-0209, B/E, = 0-0276 and 4,/f = 0-500. The resume 
values of D,/D , D;/D,) and C are respectively 0:0958, 0-0037 and 1-0995. The 
first two errors in equation (31) contribute errors in C of order 10-4 and the 
experimental measurements of the geometry also contribute to C an error of 
order 10-4. We estimate that C = 1-:0995 + 0-0005. 

To indicate how well the present theory would give absolute values of 
viscosity, we can make a calibration curve similar to that made by KEsTIN and 
WANG [1] (their value of C has a slightly different meaning from the C used 
here, however). Using accurately known values of viscosity along with the data 
of Kestin and PILarczyx [7], one can compute what value of C would give 
exact agreement between the values of viscosity computed from the data of 
Kestin and PILARczykK and the known values. Figure 4 shows some values of 
C calculated in this way. 

Certainly the theoretical value of 1-0995 is a decided improvement over the 
value C = 1 but the agreement between theory and experiment is not as good 
as expected. On the basis of (31), the theory should be correct to O-1°, 10a 
as large as about 0-4, Furthermore, the dominant error from a constant value 
should be nearly proportional to £4. The experimental measurements were 
reproducible to about 0-1% but the errors in some of the instrument constants 
might have been somewhat larger than this. Such errors would have little 
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effect on the comparative viscosity measurements reported by KeEstiIn and 
WANG but would have a significant effect upon absolute measurements. At the 
present time, the only plausible explanation that can be offered for the remain- 
ing discrepancy is that they arise from errors in the measurement of both Ay 
and the plate separations. Unfortunately there is no way of checking this from 
the available data®). ; 
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Comparison of the theoretical formula with data of Kerstin and PirarczyK. Previous theories 
neglecting the edge effect give a value C = 1. @ nitrogen; © hydrogen; x helium; A argon. 
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Zusammenfassung 


Es wird untersucht, ob eine zwischen zwei feststehenden Platten schwingende 
Scheibe als Gerat fiir die Absolutmessung der Zahigkeit von Fliissigkeiten 
benutzt werden kann. Die dafiir bestehenden Theorien, bei denen die Viskositat 
aus dem logarithmischen Dekrement der Schwingung berechnet wird, sind hier 
durch die Beriicksichtigung der Randeffekte erweitert worden, unter der Voraus- 
setzung, dass der Abstand zwischen den beiden feststehenden Platten klein ist 
im Verhaltnis zum Radius der Scheibe und zur Grenzschichtdicke. Die Theorie 
stimmt mit den von KesTIn und PILARczyK experimentell gefundenen Werten 
auf 0,1% genau iiberein. 
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The Seeding of Cumulus Clouds 
by Ground-Based Silver Iodide Generators’) 


By VINCENT J. SCHAEFER, Schenectady, N. Y.?), 
and JouHN H. Dietericn, Missoula, Mont.%), USA. 


The control of thunderstorms is a challenging objective in experimental 
meteorology. For the past five years it has been the long-range goal of Project 
Skyfire [1]*). 

The lightning storm season in western Montana, northern Idaho, eastern 
Oregon and Washington normally extends from June to September with the 
greatest storm intensity occurring in July and August. Over the past twenty 
year period approximately 70% of all fires occurring in the forests of the north- 
western United States were caused by lightning. In 1948 a start was made 
toward a study of these storms by SCHAEFER and GisBoRNE [2]. In 1952 seven 
time lapse cameras were installed on mountain summits of the Northwest to 
obtain visual records of the types of storms which occur during the summer 
nie ae ane cheat Ts National Meeting, American Meteorological Society, New York 

*) The Munitalp Foundation, Inc. of Hartford, Conn., and Nairobi, Kenya. 


i The US. Forest Service, Intermountain Forest and Range Experiment Station. 
) Numbers in brackets refer to References, page 187. 
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period [3]. Fire lookout towers such as shown in Figure 1 were used for this pur- 
pose. The following year lookout observers manning a network of 25 mountain 
tops in that area attended the first ofa series of cloud survey schools sponsored 
by the US. Forest Service and the Munitalp Foundation [4]. The daily reports 
of cloud and storm patterns radioed to headquarters of Project Skyfire laid the 
foundation for a better understanding of the nature of the lightning storm prob- 
lem in that area. 

Although this five-year study established that certain areas were likely to 
generate storms at specific locations, it also served to emphasize the great 
variability of storminess and the complexity of the storm types which occurred 
across the area. 

The studies made during the five-year Skyfire cloud and storm survey 
clearly show the great difficulties which would accompany an attempt to use 
a purely statistical approach toward establishing the success or failure of cloud 
seeding techniques for moderating lightning storms. To have any chance for 
valid results, a statistical program must be based on the premise that a suffi- 
cient number of ‘pairs’ occur which can be treated and untreated for com- 
parative testing. 

Highly variable weather conditions and the occurrence of several types of 
thunderstorms makes it difficult to obtain a number of ‘pairs’ which are similar 
in basic characteristics. Therefore, in the 1957 Project Skyfire program, it was 
decided to carry out experimental cloud seeding on each day meeting opera- 
tional specifications and to study these results in relation to unseeded con- 
ditions in surrounding areas. This approach would broaden the experience in 
experimental cloud seeding and would provide basic information for the design 
of future research. 

The test area consisted of a pie-shaped region covering about a thousand 
square miles extending about 50 miles northeast of the Montana-Idaho border 
with its apex in the vicinity of Lolo Pass in the Bitterroot Mountains as depicted 
in Figure 2. This area in turn is embedded in the much larger Skyfire network 
illustrated in Figure 3. 

Initial studies of this area were made in August 1956 [5]. Although the moun- 
tains and valleys in this area comprise a highly complex terrain, a definite behav- 
ior pattern became apparent after a month of field observations and activities. 

Based on these observations and experiences, a program was laid out for 
the 1957 season centered about an assembly of 30 silver iodide ground gener- 
ators designed by WELLS and Fuquay [6]. These were placed at 1/2 mile inter- 
vals in the form of a cross. By operating those located along the southwest- 
northeast line, a narrow silver iodide plume of high concentration could be 
generated. The generators on the northwest-southeast legs on the other hand 
~ produced a much wider plume of lower concentration. With such an arrange- 
_ ment a great variety of tests could be made. In addition to producing wide and 
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narrow plumes of relatively high concentration, the effect produced by all 
thirty generators could be studied as well as that caused by a single one. 

By late July it had been firmly established that, although most of the gener- 
ators were located near the top of a mountain ridge, little silver iodide would 
be carried aloft until after 11 o’clock in the morning. Previous to this time of 
the day, the effluent would be carried by drainage winds into the deep valleys 
extending to the southwest and northeast of the generator lines. This phe- 
nomenon might be exploited when fully understood. By operating the gener- 
ators continuously during the night and early morning hours, large masses of 
valley air might be loaded with nuclei so that, as thermals develop, the cloud 
area infected by the generators would extend over a much larger zone; including 
areas ten or more miles in the direction which would ordinarily be considered 
upwind of the generator site. 

The objective of the 1957 season was directed toward determining the degree 
of lightning supp-essi-n (if any) which might be achieved through cloud seeding 
cperations. Previous studies during the summer months in the northwest had 
established that the concentration of natural ice nuclei at temperatures warmer 
than — 20°C was either very low or undetectable. These determinations were 
made, using both surface and airborne cold chambers and by watching super- 
cooled clouds in the area. Except in extremely rare instances, evidence of ice — 
crystals in clouds never was observed until the clouds subcooled to at least — 
20°C. Consequently, the primary objective of the 1957 field program was the 
determination of whether or not clouds could be consistently nucleated on a 
massive scale within the test area as contrasted to the adjacent regions. The 
two small aircraft available to the project were to be used for aerial temperature 
and moisture sounding of the atmosphere, plume tracing, determining the 
temperature of cloud tops and bases, the concentration of natural ice nuclei, 
the behavior of clouds in regions beyond the test area and the degree of thermal 
activity over and downwind of the generator sites throughout the day. 


Project Skyfire Field Activities on an Operational Day 


In order to acquire sufficient quantitative data for analytical use in evalu- 
ating the results of the seeding of cumulus clouds from ground-based generators, | 
a considerable amount of organizational planning and activity was necessary. 
A typical day’s operation involved the field activity of twelve to fifteen persons. 
About half of the field staff had a scientific background, the rest being students 
near or at post-graduate level. 


Sounding Flight 


The first activity of the day involved an atmospheric sounding flight using 
a single engined monoplane, the Cessna 180. A recording aerograph was flown 
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to an altitude of 16000 ft shortly after sunrises The data from this flight was 
immediately transcribed and in addition a report was made of the condition 


of the sky in the general area of the test region. This information was then 
turned over to the forecaster. 


Weather Forecast 


At about 08-00 M.S.T. a special forecast for the day was completed. This 
was based on standard weather reports including facsimile weather maps and 
special local weather data such as obtained by the sounding flight of the previous 
hour and a special pilot balloon observation made in the center of the test area 
at about the time the sounding flight was being made. This was obtained by 
the field crew charged with operating the silver iodide ground generators. This 
information was supplemented by wet and dry bulb observations made at a 
group of mountain-top fire lookout stations scattered throughout the area and 
sent in by radio and telephone as the day’s forecast was being prepared. 

Based on this information, an hour-by-hour forecast was prepared. This in- 
formation was presented in schematic form and showed the probable cloud 
bases and heights likely to occur during the day, the location of the freezing 
level, the — 4°C level (where the first effects of silver iodide would occur) and 
the — 39°C level (where spontaneous nucleation occurs). The presentation also 
showed the temperature with height profile, and the wind velocities and direc- 
tions to be expected. It also included a brief forecast for the area and the loca- 
tion of the major axis of the nearest jet stream. 

The forecaster also prepared a summary of the cloud development which 
occurred on the previous day, the precipitation and the lightning storm activity 
as reported by the nearly 100 fire weather stations scattered throughout the 
mountain area and the 22 Project Skyfire special fire lookout stations. 


Briefing for the Day’s Operation 
At 08.15, the project leaders gathered at the airport for a briefing. If the 


- forecast was favorable for an operational day, the project leaders checked over 


their duties, supplies, the manpower available and, after agreeing on the pro- 
cedures to follow, departed for their respective observation posts. Since cumu- 
lus clouds in the test area rarely formed before 11.30 and the trip to the various 
locations required about an hour of driving, it was quite feasible to have the 
posts manned and all observational equipment operating before cumulus activi- 
ty was underway. At the end of the briefing, word was also sent by telephone 


to the generator operations crew so they could arrange their activity to fit into 


the day’s schedule. 


> ZAMP X/12 


ey ey 
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Manning of Observation Sites 


About 11.30, all posts were in operational order. As each station was set up, 
the word was transmitted by radio to the control station so the individual 
charged with coordinating activities knew that proper observational data 
would be obtained. During most of the 1957 field season, the control site was 
located at West Fork Butte L.O. where both radio and telephone ground line 
were available. This site is shown in Figure 4. With the latest field information 
available on the current weather and seeing the condition of the sky, the con- 
troller then telephoned the flight crew, who were in standby status, and in- 
structed them on take-off schedule and operational procedure. At this time, 
the flight observer was given a rough schedule for his operation. This would be 
supplemented by radio instructions during the flight in the area. 


Exploratory Flight Operations 


During each operational day, plans called for an exploratory flight during 
the active part of the day. The Cessna 180 would be flown to the test area, 
following instructions telephoned from the control site. Ordinarily, the flight 
was called for when the first small cumulus clouds appeared in the area. Upon 
take-off, the plane would be headed for the region above the generator lines, © 
climbing as it traveled. In this manner, when that area was reached, the plane 
would be close to cloud base altitude. The temperature and height at cloud base 
would be measured and this information relayed to the control station where 
it was checked against the forecast. Information on up-drafts was also given 
by the pilot and observer. The plane crew was then released for a period, to 
attempt to survey the location of the plume of silver iodide particles rising from 
the ground generators. Measurements were made with a light, portable, cold 
chamber, samples being taken both in the vicinity of the generators and in 
areas beyond possible influence of the generators, the latter observations being 
made to establish the background level of natural ice-forming nuclei in the air 
which might also modify the clouds, if present in sufficient concentration. The 
silver iodide plume boundary was generally established by sampling the con- 
vective air immediately downwind of the generator line and in similar locations 
and altitudes farther downwind. After the nuclei observations were completed, 
the plane was then flown into the test area and observations made of cloud-top 
elevations and temperatures. In addition, the observer was often requested to 
determine the nature of the clouds, using optical effects or taking samples of 
glaciation, virga or precipitation. They might also be requested to survey and 
report on the behavior of clouds in areas adjacent to the test region and beyond 
the possible influence of the silver iodide emitted by the ground generators. 
Occasionally the plane was utilized to do some localized seeding, using either 
an airborne silver iodide generator or small quantities of dry ice. 
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The use of the Cessna 180 proved to be eminently satisfactory. Not a single 
case of difficulty was encountered due to malfunction of the aircraft. Whenever 
a flight was needed, it was available;and by flying the plane properly, it was 
found to be eminently satisfactory for scouting, sounding, seeding and general 
observational use below 20000 ft. 


The Use of Time Lapse Cameras 


The time lapse camera provides one of the best possible records of the reac- 
tions which occur during an operational day. A number of such cameras were 
used, including a pair of full sky cameras for making triangulation measure- 
ments of cloud position, an extreme wide angle camera (84°) for covering most of 
the test area and the adjacent regions, and several ordinary time lapse cameras 
with 9-5 mm, 17 mm and 25 mm focal lengths. These were used to obtain de- 
tailed observations of reactions occurring in the test area. Inanumber ofinstances 
it was possible to orient the camera so that it included both the plume area 
and the region adjacent to it. In this manner both seeded and unseeded clouds 
appeared in the same sequences. Two-second interval periods were used for the 
reaction studies, 60 to 120 s intervals for the full sky and 84° horizon came as. 


The Use of Still Cameras 


While rate of growth and cloud location data could be obtained from the 
time lapse cameras, the 16 mm pictures are too small for illustrative use and 
ready reference. A group of still cameras were used for this purpose, including 
35 mm color cameras, specially constructed panoramic cameras providing 
about 90° of the horizon on a film 2!/,"" x 6’, several 4” x 5’ press cameras 
and a 3!/,’’  44/,’’ Land Polaroid Camera. This latter camera was used exten- 
sively to provide a current record of the typical clouds and effects observed 
during the active part of the day. The photographs obtained with this camera 
were used to illustrate the log obtained at the control site and proved to be 
extremely useful in post-briefing discussions of the operation and for classifying 
the general nature of the operational day. 


The Use of a Mobile Radar Unit 


About 91/, miles NNW of the control site at West Fork Butte, a mobile 
radar unit was located on the peak of a grass-covered knob near the top of 
a long narrow ridge close to the northern edge of the test area. This was an 


A APS-2 10 cm aircraft radar, modified so it could be controlled from within a 


small 14 ft house trailer. A retractable plastic dome made it weather proof. 


A time lapse camera photographed a slave scope and related clock, elevation 


eA 1? 


indicator and title board. This radar was activated by portable generators and 


was operated whenever there were suitable clouds within sixty miles of the area. 
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The Use of Cloud Theodolites 


During the period of 11.30-14.30, which was generally the most active part 
of the day, cloud growth studies were often made, using cloud theodolites [7]. 
In addition to rate of growth studies, the exact location and elevation of cloud 
base and top were determined, using azimuths obtained from control and ra- 
dar sites. An accurate field map at each site permitted the observers to triangu- 
late the cloud and thus have current information on the cloud behavior, proba- 
ble temperatures and the likelihood of it being affected by the plume from the 
silver iodide ground generators. | 


The Release of Free Lift and Zero Lift Balloons 


During the course of an operational day, and especially during the period 
of initial cloud activity, a number of small 30 g balloons were launched from 
the several field stations. These were used to determine air motions both in 
terms of local trajectories and the presence or lack of convective movements. 

The balanced or zero lift ballon was particularly useful at the generator site, 
especially when the lines were first operated. A 30 g balloon was inflated with 
helium or hydrogen and then weighted with a small tree branch. Twigs or leaves — 
were then removed until the balloon showed just a trace of buoyancy. It was 
then released and followed by theodolite or binoculars until lost to view. Ordi- 
narily, the balloon could be followed for seven to ten minutes, if convection was 
underway. Earlier in the day, if drainage winds were dominant, the balloon 
would either remain near the release point or would quickly be lost from view 
as it moved down into the valley. 

If, at either the control or radar site, the local winds or cloud motions 
showed a considerable variation from the forecast, a 30 or 100 g free balloon 
would be released and a pilot-balloon run made to determine any changes 
which had occurred in the wind field. 

Occasionally zero lift balloons were released from the radar site and followed 
by a 1 m range finder. This method was found to be very useful and extremely 
easy to accomplish, especially when a thin strip of 1/2’ wide 0-001” thick 
aluminium foil was used as the weight. The sparkling foil was an ideal object — 
for obtaining an accurate fix on the balloon location. 


The Use of a Tethered Balloon or Kite for Local Wind Observations 


Since the wind field in a mountainous area follows an extremely com- 
plex pattern and modifications occur quite frequently, experiments were con- 


ducted to determine if a method could be devised for making continuous 
observations. 
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A 500 g balloon was tethered by nylon cord*with smaller balloons on 10 m 
strings placed at 100 m intervals along the cord. This method provided very 
good indications of the air motions-wnder conditions of light winds. However, 
as the wind increased with altitude, the balloon tended to lay over more and 
more until its usefulness was nullified. 

A few experiments with a box kite indicated that such a device would solve 
this problem. The stronger the wind, the higher such a kite climbs, so that it 
is quite possible the air motions at least to cloud base could be continuously 
observed. Small balloons or streamers of different color placed at 200 m spac- 
ings would provide indications of the local wind field. Not enough time was 
available to develop this technique, but it would provide extremely useful data, 
if successful. 


Observations of Local W eather 


At each location hourly weather observations were made. These included 
wet and dry bulb temperatures, wind velocity and direction, condition of the 
sky, amount and distribution of cloudiness, location and number of lightning 
strokes, if any, storm tracks and precipitation. 

If any special changes occurred between the regular observation periods, 
these were included in the records. Time of appearance and location of the first 
cumulus clouds, the first glaciation, virga, precipitation and lightning were also 
noted. 


Observations of Atmospheric Electricity 


During the 1957 operation, the main emphasis of the field operation was 
directed toward the development of adequate observational techniques, to 
obtain sufficiently detailed records for reconstructing changes which occurred 
in the clouds over the test area and the immediately adjacent regions, to permit 
an objective evaluation of cloud modification effects. Although the subject of 
atmospheric electricity was given a high priority listing, it was not feasible to 
get the necessary field installations in time for the field operations. Some corona 
current measurements were obtained from time to time [8]. The results merely 
emphasized the importance of such measurements and the interesting effects 
which could be obtained if suitable instrumentation were available. This activi- 
ty resulted in the decision to pay particular attention to this subject during 
the 1958 field season. 


Data Collection — 


For each operational day a mass of data was accumulated. To keep this in- 
formation in order and to make sure that adequate records were being obtained, 
a folder was prepared for each active day’s operation. A check list fastened to 
the outer cover indicated the completeness of the file. A certain degree of lag 
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in completing the file was inevitable, especially with respect to the photo- 
graphic records. However — by using nonoperational days to work up some of 
the current data, the individuals charged with conducting the analysis quickly 
discovered the importance of closing the gaps and keeping their own contribu- 
tions up to date. This helped a great deal in having good data for the arduous 
job of evaluation which was carried out for several months after the field opera- 
tions ended. 

Based on these studies, a number of forms were developed to greatly sim- 
plify data recording. These forms not only improved the reporting results but 
served to insure that the data obtained included the information necessary for 
satisfactory evaluation operations. 


Conclusions 


The results of the 1957 field season were highly encouraging [9]. From the 
field experiences emerged a much better understanding of atmospheric behav- 
ior in the northern Rockies. As this knowledge increased we found it possible to 
anticipate certain reactions and to recognize unique effects which might possibly 
become as important to the field experimentalist in meteorology as the rain 
gauge and stream flow measurements are to the flood forecaster and hydrologist. 

In previous years attempts have been made to tag and trace air parcels 
with aerosols made of fluorescent powder [10]. Silver iodide is a far more sen- 
sitive and effective substance [11]. It is entirely feasible to produce at least 
1016 detectable particles from 1 g of material. As experience is gained of the 
unique effects caused by the introduction of silver iodide into subcooled clouds, 
it will probably become unnecessary to establish its presence by any other 
means than the unusual visible changes which occur in them. 

The 1957 results of Project Skyfire showed conclusively that it is quite 
feasible to introduce enough silver iodide into cumulus clouds by suitably lo- 
cated silver iodide ground based generators to produce unnatural effects where- 
ever the clouds were colder than — 4°C. The complete change from subcooled 
cloud droplets to masses of ice crystals occurs rapidly and produces spectacular 
changes. These are three basic changes commonly seen in clouds being modified 
by silver iodide. The first effect is the glaciation of the edges as shown in 
Figure 5. This is observed at the dissipating edges of the clouds. If ice crystals 
have not formed in these portions, the fringes disappear by evaporation within 
60 s, often in a third of that time. If ice crystals have formed, the fringes persist 
for many minutes as fuzzy diffuse areas of ice crystals forming false cirrus 
clouds. Optical effects due to the ice crystals are rarely seen. 

The second effect which follows the initial glaciation is illustrated in Fig- 
ure 6. If the cloud base is near the 0°C level 13,000 to-15,000 ft. M.S.L. (of 
frequent occurrence in the northern Rockies), the shift to ice crystals is often 
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noticed to begin at the bottom of the cloud and progress upwards until the 
entire cloud has turned to a mass of ice crystals. Natural clouds tend to show 
the change to ice crystals at the top, the base rarely, if ever, showing this change 
in phase. The illustration is an excellent example of natural and unnatural ef- 
fects, since the silver iodide plume causing the change of the middle cloud was 
from a narrow stream of silver iodide diverted by a local wind shift from a 
plume whose major axis was at right angles to this local wind development. 

The third effect is illustrated in Figure 7 which is essentially an extended 
phase of the effect shown in Figure 6. In this case a sustained thermal forms a 
cloud over an active cloud breeding spot. In addition to carrying moisture up- 
ward and into the cloud base, silver iodide is also in the rising air. As the cloud 
droplets subcool to about — 5°C, ice crystals form and grow at the expense of 
the cloud droplets. Depending on the concentration of silver iodide carried into 
the clouds, the ice crystals which form on the iodide particles either immedi- 
ately use up all of the available moisture, or, if the concentration is low, the 
crystals which form sweep up some of the droplets to form soft hail or snow 
pellets (Graupel). Figure 8 illustrates an excellent example of changes in cloud 
structure due to seeding from ground-based generators. 

By shifting clouds to ice crystals continuously at an early stage of their 
development, as rarely happens under natural conditions, it is hoped that the 
charge separation process which leads to lightning generation may be disrupted 
or at least modified to the extent the cloud-to-ground strikes are reduced or 
prevented. Whether this will be the direction that will be successfull in achiev- 
ing lightning storm modification can only be answered by additional studies. 
Besides studying changes in the nature of the cloud particles, more study of the 
electrical effects occurring in natural and seeded clouds is needed. The 1958 
Skyfire program is directed toward obtaining data on the electrical properties 
of both seeded and unseeded clouds. 
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Zusammenfassung 


Wahrend den Monaten Juli und August 1957 sind 30 AgI-Bodengeneratoren, 
die bei einem gegenseitigen Abstand von einer halben Meile in zwei gekreuzten 
Reihen auf einem hohen Kamm der Bitterrooth Mountains langs der Grenze Mon- 
tana-Idaho angeordnet waren, eingesetzt worden, um das Ausmass der Verdnde- 
rungen zu erforschen, welches erreicht werden kann, wenn unterkiihlte Wolken mit 
Silberjodid-Rauch beschickt werden. 

Auffallende Veranderungen sind wahrgenommen worden, wenn der Generatoren- 
rauch Kumuluswolken verschiedener Grésse erreichte. Die Wirkungen sind von der 
Leitstation des Projekts Skyfire auf West Fork Butte aus beobachtet und im Licht- 
bild von den Autoren, die wahrend des ganzen Sommerprogrammes des Skyfire- 
Projektes den Posten betreuten, festgehalten worden. Die Lage der AgI-Fahnen ist 
einerseits durch Verwendung von Pilot- und Schwebeballonen, die von den Gene- 
ratoren aus gestartet wurden, und andererseits durch den Einsatz von tragbaren 
Kuhlkammern, die im Auto oder Flugzeug transportiert wurden, ermittelt worden. 
Die Ergebnisse der AgI-Impfungen sind derart ungewohnlich, dass erwartet werden 
kann, weitere Studien wiirden die Méglichkeit erschliessen, durch die Beobachtung 
der tatsachlich sich einstellenden Veranderungen an Wolken, die Rauchfahne 6rt- 


lich festzulegen, insbesondere wenn die Wolken klein und nicht kalter als — 10°C 
sind. 


(Received: May 5, 1958.) 
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Figure 1 Figure 4 
Forest fire lookout cooperating with Project West Fork Butte L. O. Site of 1957 
Skyfire. Gisborne Mt. L. O., Idaho. 


control 
station of Project Skyfire. 


Wieak 


a 


TEWART 
Woe ”) a MINERAL PK 
A fre Ns 7K y = 
A 
| pae 
bs) é 
rs / 
« 
x 
, \ 
< \ 
2. > SS 
TS WASTE MIN \== 
4 al \ 
an 1 | eee 
, ae ; a, AWS >t 
= Beye Le “POINT = Xt RVER 
/ ae 4260.0. >» 
| Se z= 
U WEST FORK BUTTE Pes 
ay ee te 
GENERATORS YB \ 
= SS aa Me ‘ a 
y. ’ 
ROCKY PT { fe 2 
Q 
iy 
A Round ror Ht & 
wht: E 
aie ( ee “sr atany &)° STEVENSVILLE 
| 
VICTOR, 
| 2 
\ ) 
@®  arReorT 
4S. RANGER STATION 
A FS. Looxour 
e TOWNS o 2 


Figure 2 
Location of Project Skyfire test area in Western Montana. 
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Figure 3 


Location of test area with respect to Project Skyfire cloud observation lookout-stations. 


Figure 5 


Initial modification of supercooled cloud by silver iodide — the ice halo or glaciation halo effect. 
Note unaffected cloud beyond edge of plume. 
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Figure 6 
Effect produced in a cloud by silver iodide. A narrow finger of the plume reached the center cloud 
without modifying the clouds on either side. August 22, 1957. 


Figure 7 


; jodi i - yene Pe ota {te 
Major effect caused in cloud due to silver iodide seeding by ground generator. August 22, 1957. 
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Figure 8 


ical example of effective cumulus cloud seeding observed during a Project Skyfire operation. 
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Eine wahrscheinlichkeitstheoretische Begriindung 
der Integrationsformeln von Newton-Cotes 


Von WERNER UHLMANN, Hamburg, Deutschland!) 


Ein bestimmtes Integral einer reellen Funktion soll mit Hilfe eines Mittel- 
wertverfahrens mit aquidistanten Stiitzstellen approximiert werden. Ohne Be- 
schrankung der Allgemeinheit kénnen wir annehmen, dass das Integrations- 
intervall symmetrisch zum Nullpunkt liegt. Es soll also 


+nh 
pe +H 


f= / x()dt ndaherungsweise gleich 7 = DS a, x(t h) 
ee = —n 


sein. Es kommt dann darauf an, die Koeffizienten a; méglichst giinstig zu 
wahlen. Dafiir gibt es zwei Méglichkeiten, die naturgemass zu denselben Werten 
a; fuhren (Formeln von NEwron-Cortes): 1. Man fordert, dass fiir Polynome 


bis zu méglichst hohem Grade J = J ist. Neben den guten numerischen Erfah- 
rungen besteht eine gewisse theoretische Rechtfertigung im Approximations- 
satz von WEIERSTRASS. 2. Man verlangt, dass die Taylor-Entwicklungen nach h 


von J und 7 bis zu Gliedern méglichst hoher Ordnung iibereinstimmen. Dieses 
Vorgehen hat den Vorteil, dass es zugleich Fehlerabschatzungen liefert; dafiir 
sind aber auf der anderen Seite Differenzierbarkeitsvoraussetzungen zu 
machen, die fiir das blosse Integrieren unnd6tig sind. 

Es ist daher von Interesse, das Problem ohne solche Differenzierbarkeits- 
voraussetzungen zu behandeln. Eine Méglichkeit dazu besteht in der Heran- 
ziehung der Theorie der stochastischen Prozesse. Auf dieser Grundlage hat 
Cu. BLAnc [1]?) teilweise zusammen mit W. LINIGER [2] den Fehler fiir ver- 
schiedene gegebene Naherungsverfahren untersucht. In dieser Arbeit*) soll nun 
umgekehrt festgestellt werden, welches das giinstigste Verfahren im Sinne 
dieser stochastischen Fehlerauswertung ist und welche Eigenschaften es hat. 


1) Institut fiir Versicherungsmathematik und mathematische Statistik der Universitat 


Hamburg. : 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, S. 206. ; 
3) Der Vereinsbank in Hamburg ist fiir die Bereitstellung von Mitteln zur Durchftihrung dieses 


Forschungsvorhabens zu danken. 
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1. Problemstellung 


Es sei &(¢) ein stochastischer Prozess (vgl. [3] oder [4)). Jede Realisierung 
dieses Prozesses ist eine reellwertige Funktion x(¢) der reellen Variablen ¢ im 
gewohnlichen Sinne. Wir setzen voraus, dass die folgenden Erwartungswerte 
existieren, und zwar sei fiir alle ¢ 


E(é]=0 und E[é()])=1. (1) 


= 


Dann ist 7(¢, t’) = E[E(t) &(t’)] die Korrelationsfunktion, und wir setzen voraus, 
dass sie nur von ¢ — ¢’ abhangt: 


E[é(t) &(t’)] = RU — #) . (2) 


Wir haben also einen stationdren stochastischen Prozess im weiteren Sinne 
(von 2. Ordnung) vorliegen. Diese Voraussetzungen kénnen wohl fiir unsere 
Zwecke als sinnvoll und gerechtfertigt angesehen werden. 


Wir setzen 
Se [Rm dt. (3) 
Dann ist ersichtlich : 
R(0)=1, RA) = RG), SO]0, Se) —=sre (4) 


Die Frage, welche Funktionen als Korrelationsfunktionen auftreten kénnen, 
beantwortet der Satz von KHINTCHINE [5]. Zwei Méglichkeiten sind zum 
Beispiel: R(t) = sinw ¢t/w t, w > O (vgl. [2]) und R(t) = e-*?”, o > 0. 

Unter diesen Voraussetzungen ist (4 > 0,  natiirliche Zahl = 1) 


+nh 
+n 


n= | EO) dt — YY a, (i h) (5) 


—nh Ce 


eine zufallige Variable mit E [7| = 0, die fiir jede Realisierung x(t) den Fehler 
der naherungsweisen Integration darstellt. Es sei bemerkt, dass hier die Anzahl 
der Stiitzstellen gleich (2 7 + 1) ist. Man kann aber sofort auch den Fall einer 


geraden Anzahl von Stiitzstellen erfassen, indem man von vornherein n unge- 
rade und a,; = 0 wahlt fiir 


i oO ce Ts cate ores 


2 
74 


nur ist die Schrittweite dann 2 h. 
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Als Mass fiir die Giite der Naherung wollen wir den Erwartungswert E/[72] 
benutzen. Bei vorgegebener Korrelationsfunktion ist 7 und damit auch E[7?] 
eine Funktion der Schrittweite / und-der Koeffizienten aj. 
Man rechnet leicht nach (vgl. [1]), dass 


2nh 
7 bn 


E[r?] = 2 | S(t) dt — aD, a, {S((n + i)h) + S((n — 7) h)} | 


q Hite (8) 
+ SY” a, a; R((i — k) h) | 


kis —n 


ist. Beschrankt man sich auf den Fall symmetrischer Integrationsformeln (das 
heisst a; = a_,; fir? = 1, 2,..., ), so lasst sich (6) umformen zu 


2nh 
Ely] = 2 | S(t) dt — 4 ay S(n h) — 4 Sa, {S((n + 1) h) + S((n — 1) h)} 
6 v=1 
+ a +4 Saya, Rh) +2 7 aa, {R((i + 8) h) + Ri — A) AD}. 


| (7) 
k,t=1 

Unser Problem ist es, die Koeffizienten a; so zu bestimmen, dass E[7?] minimal 
ausfallt. Jedes (2 + 1)-Tupel (a_,,..., @,), fiir das E[n?] = Min ist, bezeich- 
nen wir als Minimallésung. Bei gegebener Korrelationsfunktion hangt jede 
Minimallésung ersichtlich nur noch von h ab. 

Es wird sich zeigen, dass jede Minimallésung naherungsweise gleich den 
Koeffizienten von NEwTon-CoTEs ist, und zwar unabhangig von der — immer 
etwas schwierig zu rechtfertigenden — Wahl der Korrelationsfunktion. Uber- 
dies wird bewiesen werden, dass E{7?] fiir die Minimallésung von derselben 
Gréssenordnung ist wie fiir die Formeln von Newron-Cotes. Zugleich erhalt 
man Schatzwerte fiir den Fehler, die den bekannten Restgliedern bei NEwWToN- 
Cores entsprechen, die aber statt der héheren Ableitungen des Integranden 
hdhere Differenzen enthalten. 

Zuvor wollen wir einen Zusammenhang zwischen der Korrelationsfunktion 
R(t) und den héheren Differenzen von é(t) herstellen. Wie tiblich setzen wir 

we 
A*E(t) = > (*) (-1y’ &(t + (u—v)h) fir =0,1,2,.... 


v=vU 


Dann ist offenbar 
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Wir setzen nun voraus, dass R(t) (2 + 2)-mal stetig differenzierbar ist und 
sich also darstellen lasst in der Form 


lh 
R(t) = YY be B° + O***) , (8) 
o=0 


wobei b, = 1 wegen (4) und O das Landau-Symbol ist. 
Wir wollen zeigen, dass fiir 4 = 0, 1, 2,... gilt: 


E((A"E(t))"] = (1)! (2p)! by, WPM + OPE?) . 9) 


Beweis: Fiir = 0 ist (9) offensichtlich richtig. Fir « 2 1 benutzen wir, 
dass 


“ (—1)* w! firt=yp, 
Focyr-| 
v=0 O fir tT<p 


ist, was sofort aus dem Mittelwertsatz fiir hdhere Differentialquotienten ange- 
wendet auf die Funktion y = x7 folgt. Nun ist unter Verwendung von (8): 


E{(A*é(t))"] = 7 | Ag) (aN aad ae a 0 | by #2 + O(W2#*?) 


v,o—0 


+ O(R#*?) = (—1)" (2)! by, WM + O(RPH*?) . 


2. Existenz und Charakterisierung der Minimallésung 


Nach (6) ist der Erwartungswert E[7?] als Funktion der (2 + 1) Koeffi- 
zienten a; 1m ganzen (2 + 1)-dimensionalen Raum definiert. Da 7? = 0 ist, 


. . . . pall 
ist es auch E{7?]. Folglich ist die quadratische Form >” a, a, R((i — ) h) 
k,t=—n 


positiv semi-definit und die symmetrische Matrix (Ri — k) h)), deren Eigen- 
werte also = 0 sind, lasst sich nicht-singular auf Diagonalgestalt transfor- 
mieren, Durch diese Transformation geht (6) itber in eine Form, die — von einer 
additiven Konstanten abgesehen — nur noch quadratische Glieder enthalt und 
der man also unmittelbar ansehen kann, dass E{y?] sein absolutes Minimum 
annimmt, und zwar auf einer linearen Mannigfaltigkeit, deren Dimension 
= (2+ 1) — [Rang von (R((i — k) h))| ist. Weiter folgt, dass diese Mannig- 
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faltigkeit genau die Lésung des linearen Gleichungssystems 


¥ a, R((i — 7) h) = S((n + 7) h) + S((n — 7) h) | 


i=—n 


(10) 
(=H, ...,—1,0,1,...,.7) | 


ist, das man aus (6) durch partielle Differentiation von E[7?] nach den a; erhiilt. 


Damit sind die Minimallisungen in bekannter Weise als Lésungen von (10) 
charakterisiert. Sie sind nicht notwendig eindeutig bestimmt, aber jede Lisung 
von (10) ergibt denselben Wert von E{n?). 


Fiir die Praxis interessieren nur symmetrische Lisungen (das _heisst 
a; = a_,), und in der Tat rechnet man sofort nach, dass zu jeder Lésung von 


(10) — sagen wir a‘_,,..., a, — eine symmetrische Lésung existiert, namlich 

a; = a_; = (a; + a’_;)/2. Fir symmetrische Lésungen reduziert sich (10) auf 

ay RG) + 3° a,{R(E+7) k) + Ri) A} eee, 
t=1 [ 


= S((n +7) h)4+ S((n—9)h) G =0,1,2,..., 0), 


und man hat dann also nur noch (m + 1) lineare Gleichungen zu lésen. 

Fiir den Fall einer geraden Anzahl von Stiitzstellen (1 ungerade, @,; = 0) 
sind nur die (nm + 1) Koeffizienten @_,,... , @1, &, a, ..., @, zu bestimmen. 
Ersichtlich bleiben die obigen Betrachtungen sinngemiss richtig, nur entfallen 
natiirlich in (10) und (11) alle Gleichungen fiir gerade 7. 


3. Die Integrationsformeln von Newton-Cotes 


Bevor wir die Minimallésungen naher untersuchen, wollen wir einige be- 
kannte Eigenschaften der Formeln von NEwron-Cotes zusammenstellen (vgl. 
etwa [6], teilweise auch [8)). 


1. Ungerade Anzahl von Stiitzstellen 


x(t) sei im Intervall (— 7 h, + hj eine (2 + 4)-mal stetig differenzierbare 
Funktion. Das Integral 


+h 
[xo dt 
—nh 


soll angenahert werden durch 


Ny ha,x(ih) mit a=a,. 


t=—Nn 
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Aus der Taylor-Entwicklung des Fehlers 


+nh fu 


y(h) = | xi di — ST hee; x(6 b) 
—nh eas 


nach # erkennt man unmittelbar, dass die «, eindeutig bestimmt sind als Lésung 
des linearen Gleichungssystems 


+n 


j20 = 2 2o+1 aD Ss 
oT, a; = ore i pee uae orp ON ony Te (12) 
und dass 
y(h) gta | nee Y anttg,| f2+3 x2n+2) (0) + O(n2n+8) (13) 
: (2n-- 2)! (20 3 ~ i 


ist. Unter Benutzung eines Interpolationspolynoms fiir x(t) ergibt sich anderer- 


seits 


Z 
y(h) = ee ih 7 AAN+3 y2n+2(q) , 


wobei t passend aus [— h, +7 h] und 


ar 
rr=| [[@—7) at+ 0 (14) 
qj 7 
ist. Also ist auch 
y2nts n ‘ ; 
inoue ae Oz (15) 


2. Gerade Anzahl von Stiitzstellen 


Hier sei x(¢) (n + 3)-mal stetig differenzierbar im Intervall [—nh, +n h), 
wobei 7 ungerade ist. Das Integral 


soll angenahert werden durch 


vi h a; x(t h) mit a; = Kx und Oa =O ir ee i aero hae. 


1=—%H 


Carre . : 
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Die verbleibenden «; sind eindeutig bestimmt als Lésungen von 


ME Sen os 1 mere +1 ae Were 
pes = ee a Mies 0. 1,2 ..., ——, (16) 


t=] 4 


wobei wir die Schreibweise 2 {7 fiir 7 gerade bzw. 2 + fiir 7 ungerade benutzen. 
Fir den Fehler 


tnh 
~ bm 


y(h) = | x(t) dt eae h x; x(i h) 


J 
nh % 


oti 


ergibt sich 
2 y 2 y(n 2 yp pn+2 4(n 
y(h) = (iy ABT i yesd (7) — Tait Ae +? xm+0(0) + O(h"+4), (17) 
wobei t passend aus [—” hh, +” h] und 
so wa = ain 5 ie : “ 
re a aa aaa i US tear Re (18) 
2H 0 24) 


ist. 


4. Taylor-Entwicklung der Minimallésungen 


Wir denken uns jetzt die Korrelationsfunktion R(t) vorgegeben. Dann 
hangen die durch (10) bestimmten Minimallésungen nur noch von / ab. Wir 
wollen beweisen, dass unter gewissen Voraussetzungen die Taylor-Entwicklung 
jeder Minimallésung mit den Koeffizienten von NEwTon-CorTes beginnt, dass 
also die Koeffizienten von NEwTon-CoTes naherungsweise Minimallésungen 
sind. 

Zunichst wollen wir den Fall einer ungeraden Anzahl von Stiitzstellen 
behandeln. 

Wir setzen voraus, dass R(t) (4 + 5)-mal stetig differenzierbar ist. Dann 
lasst sich R(t) schreiben in der Form 

2n+1 
RQ) = DF ba, At + O(han +4) (19) 

1-0 

mit by) = 1. Denkt man sich das Gleichungssystem (10) bzw. (11) nach dem 
GauBschen Eliminationsverfahren aufgelést und die fre: wahlbaren a,(h) hin- 
reichend oft differenzierbar gewahlt, so erkennt man, dass alle a,(/) ebenfalls 
(4 + 5)-mal stetig differenzierbar sind und dass die Taylor-Entwicklungen 
der a,(h) nur ungerade Potenzen von / enthalten. Wir konnen also schreiben : 
a,(h) = 2} et, pit A Anise S (20) 


eS 
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Weiter setzen wir voraus, dass fiir 7 = 1,3,5,...,2”+1 


det [(H 1+ 0 (21) 


sb Ae by. 


ist; dabei durchlauft der Zeilenindex yu die Werte tT, tT + 2, xs, 20 =] bond der, 
Spaltenindex y die Werte 1, 3, 5,. 

Diese Voraussetzung (21) ist zum Beispiel erfiillt fiir die im 1. Abschnitt 
erwahnten Korrelationsfunktionen, wie die folgenden beiden Hilfssdtze zeigen. 


Hilfssatz: Fiir R(t) = sinw t/w ¢ mit w > 0 ist fiir alle t (t ungerade) 


Dia det |e 5: 


ee ev 


a0. 


Beweis: Fir 


: 20 
RQ) = =5* = D1) (2 a Ter 


wot =O 
ist 
oh Sie a ee 
park lena (u — 7)! (x — v)! (4 Sys Ty 


Da es uns nur auf das Verschwinden bzw. Nichtverschwinden von D ankommt, 


bezeichnen wir mit c,, ¢y,... passende reelle Zahlen +0. Wir kénnen dann — 


schreiben ; 
D =e, det (_, = at 


Wir setzen (t+ 1)/2=m, w=1t+21—2,»=1—2k+ 2 und erhalten 


i 
De, det ant op Saez 7) ; 
wobei 7, k = 1, 2,3,..., m. Um zu beweisen, dass D + 0 ist, zeigen wir gleich © 


allgemeiner : Fiir jede natiirliche Zahl m = 1 und jede natiirliche Zahl M = list 


M 
1 
vet (TT ea Ee rae (22) 


wobei der Index m andeuten soll, dass 7 und k die Zahlen 1, 2, ... , m durch- 


laufen. Wir beweisen (22) durch vollstandige Induktion nach m. Fiir m = 1 ist 
(22) offensichtlich richtig, und fiir den Schluss von (m — 1) auf m benutzen 
wir, dass bekanntlich fiir jede beliebige Determinante 


det (a n= ( Tn ) dat (ef tien. Gere ee (23) 


Im, k am, k+1 


ist, wenn nur @,,, + O ist fiir k=1,2,..) ,m (m> 1). Damit rechnet man 
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sofort nach, dass 


1 M 1 (M-+1 1 
clet ee $2 = ; 
Ussca¥ezs) é det ( JT aT) 
d “ m—1 


ist, und die rechte Determinante ist + 0 nach Induktionsannahme. 


Hilfssatz: Fiir R(t) = e-" mit o > 0 ist fiir alle t (t ungerade) 


D, = det [(4=”) b, J £0. 


L\—t/ “uy 


Beweis: Fiir 


cl) eo cane ee meee ay Pe 


a_i 
ist 
(forty dos = (a =v} (—1)7¢ (—1)-"? 


eT ‘ 


w* (= M/E ua) 2 ot fe 12°F” 


Mit den oben eingefiihrten Bezeichnungen ist also 
e 


~ (u — »)! (2i+2k—4)! 
D, = ¢ det (Arr) = ea det ( Crk al jae 


Wir fiihren den Beweis durch vollstandige Induktion nach m. Der Fall m = 1 
ist trivial, und fiir den Schluss von (m— 1) auf m benutzen wir (23). Man 
rechnet leicht nach, dass A 


: xi ie 
Dy = oy det (2i +2 4)! (2i+2k elite 


(G+ &— 2)! ),, = eadet ( Gi +b 2)! 


ist, und die rechte Seite ist + 0 nach Induktionsannahme. 
Wir behaupten nun: 


Satz: Unter den einleitend angefiihrten Voraussetzungen (19), (20) wnd (21) 
sind die a,;,, (i =—n, ..., +n) die Koeffizienten von NEWTON-COTES, das hersst, 
ste gentigen den Gleichungen (12). 


Wir beweisen gleich etwas mehr: 


Satz: Unter den Voraussetzungen (19), (20) und (21) erfiillen die a; de 


folgenden Gleichungen fiir t = 1,3,5,...,2n+1(n 21): 


+n 


{7-1 = > n , (24) 
‘ : T 
i=—n 
+n 
yi? o,,=0 fir y= 3,5,....7, (25) 
t=—n 


- ZAMP X/13 


 —-S. *” 
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+H 


+n 
(i 1) Gol = 2, ORE) = 
(a a Dies >? ay ng pi: f= 4 Oe v t %; , 
fe 


i=—-n 


fiy 6 =1,3,5\4., 6 Wide = 1 eae i 


Beweis: Es sei a;(h) eine symmetrische Minimallésung der Form (20). Die 
a,(h) geniigen also dem Gleichungssystem (10), und zwar identisch in h. Wir 
differenzieren nun (10) u-mal nach h und erhalten fiir 7 = —n,..., +n” 


SY YO ay (pyr RO (Gh) 

= (w+ HRY ((n + j) h) + (mt RY ((n— jh) 
und insbesondere fiir h = 0 
EY YO a = jy? ROM) = [in +)!" + (=H R“MD). 27) 


Wegen a)(0) = 0 fiir » gerade und R)(0) = 0 fiir » ungerade verschwinden 
die rechte und die linke Seite von (27) fiir uw gerade. 
Fir ungerade uw und 7 = —n,..., +m folgt aus (27) 


t=—n vy=0 


2h 
und das ist gleichbedeutend mit (beachte «; , = «_; ,) 


Wir setzen nun o = « — o und erhalten damit 


S ¥ aC es g aK; » Drgay ‘<i [ (2 al 7) ais (m > 1)"] ae , (28) 
u-1 ; | " U-oO +n I . 
je lane en oy ju-P-e af 2 ene es | == OY (29) 


Sr nlecd i. Doe i -$ tS rar “| : 


o=1 1=—-Nn = 1=—-N 

246 | 24 : 

fiir « ungerade und j = —n,..., +n. (30) 9 
Nun ist die linke Seite von n (30) ein Polynom in 7 vom Grad S yp — 1 mit 

(2+ 1) verschiedenen Nullstellen (j = —n,..., + n). Daher verschwinden 


fiir w < 2m + 2 alle Koeffizienten, und es preter eich also 


+1=—-n 1=—nNn 


o[Er-se] Senn Ben lao onl 


Fe em 35,00 fina Co 


. 


Es 
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Unter Benutzung von (30) und (31) beweisen wir nun (24), (25) und (26) 
durch vollstaéndige Induktion nach r. 

Verankerung: Fir t = 1 folgt (24) sofort aus (31) fiir o = jo = 1. (25) ent- 
fallt fiir t = 1, und (26) ist wegen (31) schon bewiesen fiir f#=1,3,...,2n4+1. 
Wegen (24) verschwindet in (30) das Glied mit o = 1. Fiir b=2n-4 3 ist 
daher (30) ein Polynom in j vom Grad <2n mit (2+ 1) verschiedenen 
Nullstellen. Es verschwinden also alle Koeffizienten, das heisst, (26) gilt auch 
fiir p= 2n+ 3. 

Der Satz sei nun richtig fiir alle ungeraden Zahlen < t (1 <7 <2n+ 1), 

Schluss von (t — 2) auf t: Wir betrachten Gleichung (26) fiir o = t. Wegen 
] =2n + 1 ist sie insbesondere erfillt fir «4 =+,1+ 2,...,r 1+. Dies 
kann man auffassen als homogenes lineares Gleichungssystem 


t+ 1 . : 
(i 2— Gleichungen) fiir die a « Unbekannten » 


+n 9 +n +n 
ptr panetee aah aos “TT 
as, 24 F |. a paar: het 8 Site de CN 
i=—n : 


i=—n i=—n 


Da die zugehérige Determinante det {(/—”) b,_,] als + 0 vorausgesetzt wurde 
(21), gibt es nur die triviale Lésung, das heisst, (24) und (25) gelten auch fiir t. 
Daraus folgt weiter, dass in (30) alle Glieder verschwinden fiir o = 1, 3,...,T. 
Also ist fiir ~ = 2” + 2+ 7 die linke Seite von (30) ein Polynom in 7 vom 
Grad <2n mit (2” +1) verschiedenen Nullstellen. Folglich verschwinden 
wieder alle Koeffizienten, und damit gilt also (26) auch fir w= 2n+2-+1, 
was zu beweisen war. 

Die Gleichungen (24) und (25) schliessen zunachst nicht aus, dass die Koef- 
fizienten von NEwTon-CorTes nicht nur naherungsweise, sondern auch exakt 
eine Minimallésung bilden. Dass dies nicht der Fall ist, zeigt sich schon fur 
m = 1. Man rechnet namlich leicht durch Auflésung der Gleichungen (11) nach, 
dass die Minimalldsung (falls 6 b, + 03 ist, ist sie fiir = 1 eindeutig bestimmt) 
die Form hat 


1 2S(h)| 
R(h) S(2H), 1 2 fb, £30 5, 
Bald i ea(8) = | ea ae alt sts re ee 
R(h) 1-+ R(2h) 
Sr eevee ce ee Se a 
BM aca ERE Lo Se es GSD eh Ne 
‘Rg 14.82 h) 


Fiir R(t) = sinw t/w ¢ ebenso wie fiir R(/) = e"/2 sind die Koeffizienten von 
_h® sicher + 0. 
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Wir wollen nun noch in aller Kiirze auf den Fall einer geraden Anzahl von 
Stiitzstelien eingehen. Es sei also m ungerade und 


e y 
dy; (hy-S Opti iesl0y ae eos +455. 


Wir setzen voraus, dass R(t) (2 + 3)-mal stetig differenzierbar ist, und be- 
nutzen die Bezeichnungen (19) und (20) unter sinngemasser Abanderung der 
Summationsgrenzen. Die Voraussetzung (21) sei erfullt fur 7 = 1,375; 5,72 
Dann lisst sich wieder zeigen, dass die Ol die Koeffizienten von NEWTON- 
Cores sind, und zwar indem man allgemeiner den folgenden Satz beweist: 


Satz: Unter den oben angegebenen Voraussetzungen erfiillen die Ob; y 2 Oe 


(See cs atl [ila =a dh Ds Dys ens TEN ungerade = 1) die folgenden Glei- 
chungen: 


i! 
Pipi Nap ee ee 
ay. Oe Natal A (32) 
24% 
Dit Se Oat Si 3 eae (33) 
t=1 
24% 
(ga) “ual Pas a; de n° + » (eae) be, P ie ok a; | oy 0 | 
4=1 v=3 1=1 
24% 2Q4y 243 (34) 
flr 0 = 13, 5)n.0 he Und! 1S 5, ee eee | 


Da der Beweis vollig analog zum vorhergehenden Satz erfolgt, konnen wir hier — 


auf seine Wiedergabe verzichten. 


5. Erwartungswert des Fehlerquadrats fiir die Formeln 
von Newton-Cotes 


Wir wollen hier den Erwartungswert E{7?] fiir die Formeln yon NEwTon- | 


Cores berechnen, einerseits um die Analogie zu den Fehlergliedern (13) bzw. 
(17) aufzuzeigen, und andererseits um im nachsten Abschnitt einen Vergleich 
mit E[7*] fiir die Minimallésungen durchfiihren zu kénnen. Ausserdem ist es 
fiir die Anwendung der Integrationsformel von Interesse, den tatsichlich zu 
erwartenden Fehler mit Hilfe des Erwartungswertes E [7] zu schatzen. 


Satz: Wenn die Korrelationsfunktion R(t) (4 n + 8)-mal stetig differenzterbar 
und die Anzahl der Stitzstellen (_nh,..., —h,0,h, 2h,... , nh) ungerade ist 


, 


: 


| 
| 
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so gilt fiir den Fehler n {siehe (5)\ bei Benutzung der Koeffizienten «,; von NEWTON- 
COTES 


mF | ? 4a. aye ments — 
EWl= lac at ters - Br sn a(t 
x (4+ 4)! Dgn gg hAMt8 + O(hAn+8) 
| ; : (35) 
= 2 ER: pan 2 ; 
~ \(2n+ 2)r | fore me." bes | 
x BE [(4aaee E(t) )?] A O(ht+8) : 
Dabei ist 
2n+3 
RH = >» by, h?* + O(hin+8) ; (36) 


+=0 


Beweis: Wir setzen a;(h) = a_,(h) = ha; (x; unabhangig von h) in Formel (6) 
ein und erhalten mor 0" (37) 


_ 2 She {(n + 1) R((m + 2) h) + (nm — 1%) R((n — i) h)} 
Led hajayR(G— 2) pied riers (i —k) R((i—&) h), 
ijk=—n i,k=—n 
E'(y?] |,-0 = 9 - (38) 


Allgemein ergibt sich durch o-maliges Differenzieren (9 = 2,3, ...,4"+ 7): 
E®[y?)] = 2 (2 n)? R°-®) (2 wh) 


= 29 Sa {(n + ie? REM((n +i) h) + wi)" RE-*((n — 1) b)} 


+=—n 


7, ay? ha, {(n + i)? R°-((n + 4) h) + (m — i)? RE-Y((n — 1) A)} 


+o(e—1) Dy %i a, i = k)e-? Re-?)((i — k) h) 
t,k=—n 
+20 Sha, a (i — ke? RE-?((i — R) h) 
i,k=—n 


= Ss Rh? x; oy (i — k)? R®((i — k) h) 


1,k=—n 
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und 


E®[y?] |,-9 = 2 R& (0) lo n)? — ay a, [(n + 4)8-? + (nw — 4/27] | 


| oer 


+ (8) Savas — Het. | 


ik=—n 
Daraus folgt sofort: 
E®(y?)| |,» =0 fiir o ungerade. (40) 
Mit Hilfe von 


und einiger elementarer Umformungen kénnen wir aus (39) fiir gerade 0 folgern: 


e—2 


Oe aco = Re™ (0) a (3 : fl 


t=—n 


+n 
+1) Say — 20 | 


Blo (41) 
(oe —o—1) Ds, Rs ee oo : | 


k=—n 


Wir nutzen nun aus, dass die «,; die Koeffizienten von NEwTron-CoTEs sind und 
also den Gleichungen (12) geniigen. Dann folgt aus (41) sofort: 

E®) (7?) yoo =O Ur pie 4 5 6 ee (42) 
und 


2 
4n-4 +n 


E@n+9 [72] |, 9 = ROn+9(0) é . a 3 E n + 3) ps G2N+2 gy, 2 y2n+s f. (43) 


1=—N 


Mit Hilfe von (37), (38), (40), (42) und (43) ldsst sich nun die Taylor-Entwick- 
lung von E{n?] aufstellen, und man erhalt den ersten Teil von Formel (35). 
Der zweite Teil folgt sofort aus Formel (9) fiir 4 = 2 + 2, was zu beweisen 
war. 

. Am Beispiel = 1 sei gezeigt, dass man in Formel (35) im allgemeinen 
nicht auf O(A4"* 8) verzichten kann. Man rechnet namlich leicht fiir » = 1 nach: 


224 9 
E[r?] = = by hi + se biy 2 + O(h14) 


il 
= ope # E((A* &(¢))?] — a by) M2 + O(AA4) . 
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Ganz entsprechend wie oben beweist man: 
Satz: Wenn die Korrelationsfunktion R(t) (2n + 6)-mal stetig differenzter- 
bar und die Anzahl der Stiitzstellen (—nh,...,—3h,—h,h,3h,...,” h) gerade 
ist, so gilt fiir den Fehler n bet Benutzung der K oeffizienten a; von NEWTON-COTES 


t=1 


- 2 n+2 ae é. 2 
7 tae | writ pre ies as (2+ 2)! bay ight t4 4 BS ies 


Jf 2 ynrte2 an = |" 
meat |wee Oy Pd] | Blan GO) + Onte 


Man beachte dabei, dass fiir die Integration die Schrittweite = 2h, fiir A" +1 &(t) 
aber nur = / ist. Wenn man will, kann man das leicht abandern, indem man 
etwa in (9) h durch 2h ersetzt und Formel (44) entsprechend abandert. 


6. Erwartungswert des Fehlerquadrats fiir die Minimall6sungen 


Wir wollen nun beweisen, dass der Erwartungswert E[7?] fiir die Minimal- 
l6sungen von derselben Gréssenordnung ist wie fiir die Koeffizienten von 
NewrTon-CotEs [vgl. (35) und (44)]. Damit kann von unserem Standpunkt aus 
fiir praktische Anwendungen die Benutzung der Formeln von NEwron-CoTEs 
als gerechtfertigt angesehen werden. Denn bei gleicher Grdssenordnung des 
Fehlers haben die Minimallésungen den Nachteil, komplizierter zu sein, und 
ausserdem hiangen sie von der nur schwer zu treffenden Wahl der Korrelations- 
funktion ab (vgl. dazu auch fiir 7 = 1 Abschnitt 4). 


Satz: Die Anzahl der Stiitzstellen (—nh,..., —h,0,h, 2h, ..., wh) set un- 
gerade. Die Korrelationsfunktion R(t) sei (4 n + 7)-mal stetig differenzierbar, und 
ease fart = 1, 3;5, << -2n+ 3 


det [(42”) b, J+ 0. (45) 


pat 


Dann ist fiir jede (sicher existierende) symmetrische, (4 n + 7)-mal stetig differen- 
zierbare Minimallésung der Erwartungswert E(n?] ein O(h4"* 8), aber kein O(h*" +"). 


Beweis: Wir benutzen sinngemass die Bezeichnungen (19), (20) und (21). Es 
geniigt offenbar zu zeigen, dass fiir die Ableitungen von E [7?] nach hf gilt: 


EM?) | 9 =O fin r= 0,1,2,..., 44 9 (46) 


und 


EGn+9 (7?) |, 9 + 0. (47) 


(44) 
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Setzt man in (7) eine symmetrische Minimalldsung a,(/) — sie geniigt (11) — 
ein, so erhalt man 


2nh 


Ely?) waft ( de — 2S" a S((n +2) h). (48) 
Daraus folgt 
Ely, =e (49) 


E®[r?] =4n S(2nh) —2 Xa S((m + i) h) + a,(h) (n + 1) R((n +2) h)}, 


i=—n 


E®[y?] est =O (50) 


nud alisemem tur ye = 12a 


EM Dy] = 2 (2 nj? RED (2 nh) 
+m perl 


— 2. dui On (M33) al (A) (m + aye tt? SH tI" (nm + 4) h), 


t=—n v=0 


EU+N (yp?) [sts a? (2 n)i" +1 R—(0) 


+n le (51) 
ae Died pl valde cana gee mee 
i=—n v=0 
Fiir gerade ys ergibt sich 
EMV yy ae =0, (52) 


und fiir ungerade mw erhalt man 


EMM |p 9 = 2 (nye? (uw — IID 


ul 


+H 


aes ey ae . es 
2 5 MEDyeiy Dy ebayer Tal ty 
ote i=—n 
=4 3) pa cig) 
o=1 
240 
um pt+l—y ‘ 
1)! . +H 
—2Z (uw + 4h u+1—p e a+ 1—y— 
a & remy wena ey fen 2S 
24y 2 4p te ee 
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nw 
2) 
OI 


Wir setzen w + 1 — 9 = o und erhalten durch ¢lementare Umformung 


( Lyt 
| saa ol co = — 2 pa Fata + a = hog 


a) > y 
u-1 s Y -o—1 Aten ott ¥ S ’ ;o—v | 
i O.% b Ky o ” i > ce 4) ot v v ay u 


e=2 | 
Pee atid a 
( 


Nach (26) folgt 
E+ fy) |, 9 =O fir w=1,3,5,...,4n43. (54) 


Durch (49), (50), (52) und (54) ist (46) bewiesen. 
Die Behauptung (47) beweisen wir indirekt, machen also die Annahme, dass 


EGY ttl Ly = 0 (55) 
ist. Wir setzen 
+H 
1) 1 
eS ee F ul eo se Ae +h =o) 2), tras 
rng ees t=—n 


und erhalten damit aus (55) und (53) 
4n+6-—cC Can+5,0 


Wegen (24) und (25) ist c,,,5,=—0 fir o=1,3,5,...,2%+1, und also ist 


4n-+5 n+ 6— 
o ni” 6—o . 
4An 6—q 4n+5,6 
a=2n+3 r 


20 


Nach (30) ist — wieder unter Benutzung von (24) und (25) — ausserdem 


GN (56) 


4n+5 
eR ON! Pi gl es aa he eh S alice lg 
ao=2n+3 
240 
Also ist 
4n+5 
4n+5—o 7 
P(t) PD ste ° (445,60 
o=2n+3 
246 


ein Polynom in ¢ vom Grad <2n+2 mit den Nullstellen |¢| = 1, 2, 3,. 
und (wegen o ungerade) der doppeiten Nullstelle ¢ = 0. P(t) lasst sich ree 


Ss = 
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schreiben in der Form 


Nach (56) ist nun 


und wegen (14) folgt daraus c = 0, das heisst P(t) = 0, und also C4n45,6 = 9 
auch firo =2n+3, 2n+5,...,4n+5. Zusammen mit (26) ergibt siem 
damit insbesondere 


sn o +n ae 
(29 ha Sera — 2a] Bey a | Sore, [eo | 
te (57) 


i=—n 3 1=—n 
fir g=2n+3 und w= 20-3, 207 Op eee 


Nach Voraussetzung (45) folgt daraus, dass alle eckigen Klammern in (57) ver- 
schwinden und insbesondere also, dass 


+n 
2 
/2 2 as 2 Mae: 

1 Diy, ary a. 2n+3 nN fF == (() (58) 


ti=—n 


ist. Nun sind nach Abschnitt 4 die «, , gleich den Koeffizienten «; von NEWTON- | 
Cores, und also steht (58) im Widerspruch zu (15). Wir haben also (55) zum ~ 


Widerspruch gefiihrt, was zu beweisen war. 
Ganz analog beweist man: 


Satz: Die Anzahl der Stiitzstellen (—nh,..., —h,h, 3h,..., wh) set gerade | 


(also n ungerade). Die Korrelationsfunktion R(t) sei (2+ 8)-mal stetig diffe- 
tenzierbar, und. es set fuy t= 1,3,5, 1., War 2 


det [(#=”) b 


H—t 


al #0. 
Dann ist fiir jede (sicher existierende) symmetrische, (2 + 8)-mal stetig diffe- 


renzierbare Minimallosung der Erwartungswert E{n®] ein O(h2"+4), aber kein 
Open >). 
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Summary 


A definite integral is approximated by a linear combination of values of the 
integrand at equidistant points. As a measure for the goodness of fit we use the 
mean square deviation; for this we take a wide sense stationary stochastic process 
as basis. We prove that there exist always a best integration formula in this sense. 
The Taylor expansion of its coefficients begins with the coefficients of NEwron- 
Cotes. It is shown that the error of the best integration formula is of the same 
order of magnitude as in the well-known formula of NEwton-Cortes. Therefore we 
get a justification for the formulae of Newron-Coress also if the integrand is not 
differentiable. Besides we obtain estimates of the error, which use higher differences 
instead of higher differential quotients of the integrand. 


(Eingegangen: 11. August 1958.) 
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Zeros and Poles 
of Output Voltage of 3-Terminal Potentiometer Networks 


By Moues Aziz AppeEt-Messin, Alexandria, Egypt’) 


1. Introduction 


In the past few years resistive networks containing fixed resistors and linear 
potentiometers?) have been developed for the construction of analogue computing 
instruments. GREENWoopD, Horpam, and Macrae [1]*) have used potentiometer 
networks for the generation of some prescribed functions. The paper of SPEISER [2] 
contains methods of performing the elementary arithmetic operations and of 


1) Faculty of Engineering, Alexandria University. é 

2) Hereafter such networks will be referred to as potentiometer networks. A linear potentio- 
meter is that in which the relationship between the shaft angle and the variation of resistance 
between the contact point and one end of the winding is linear. 

3) Numbers in brackets refer to References, page 215. 
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curve fitting for functions of one and of two variables. The work of ABDEL- 
Messin [3] comprises all possible networks containing one and two potentio- 
meters for one and for two variables, and gives a design for a network to generate 
the ballistic functions in air. 

Consider a 3-terminal network (Figure 1) containing fixed resistors and a finite 
number of linear potentiometers mounted on the same shaft. Let a represent the 
angular displacement of the shaft from one end of the potentiometer winding, 
expressed as a fraction of its fullrange;0 S¥ Sl. 


= 


Figure 1 


3-terminal potentiometer network. 


If the network is supplied by a constant voltage £,, then the voltage at any 
point of the network, output voltage, is a function of 7; E,(%). SPEISER [2] has 
stated that £,(%) is a rational function of ¥, which has no poles in the interval 
0 = % = land no zeros within 0 =< 47 =< 1; 

In this paper we give more precise properties of the zeros and poles of the 
function £,(%), stated in the following theorem, and prove these properties. 


2. Theorem 


The output voltage E,(%) of a 3-terminal potentiometer network is a rational — 


function of x with the following properties: 


(i) The zeros of E,(%) are either real outside the interval 0 < x < 1, or occur 
in conjugate complex pairs. 
(ii) The poles of E,(7) are simple real negative or > 1. 


3. Proof of Property (i) 4) 


In a potentiometer network the resistance elements appear in the form 7%, ¥ 7%. _ 


and (1 — *) v3, where 0S ¥< 1. 


If three such resistances are connected in parallel as in Figure 2, then their 
total conductance is 


a) 4 1 : 1 
Vo he (ils Sophie 
Putting 
1 
tai Pa sieeve 


‘) The problem is treated in a method analogous to that used by FairKow and Grrsr [4] for 
the transfer function of RC networks. 


| 
: 


——— 
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then 


where 


1 
Bee — fe mee ah ae = Die , €=— (a,b,¢c20). 

15 
The network of Figure 1 is considered on nodal system. The ground terminal 
is taken as node 0, the other input and output terminals as nodes 1 and 2 respec- 


node i node J 


Figure 2 


Parallel connection of resistance elements. 


tively. The remaining nodes are identified so that each branch is an 7, an * 7, an 
(1 — *) rv; in parallel as in Figure 2. The conductance g,, (i +7) of the branch 
between nodes 7 and 7 is according to (1) of the form 


ot eas (C¢5 + ep ¥ + Or, ¥*)  (AeGs O35, C55 2 O) - | 


Denote 


n (2) 
Dy bis = Bi (eee aera ' | 


7=0 
¢42% 


where + 1 is the total number of nodes. (Of course g;; = g;; and g;; has the 
same form as g; ;.) 

The only driving current source is J, impressed on node 1. Let E; denote the 
voltage at node i measured with reference to the ground node 0, 

Applying KrrcnHorr’s current law for the equilibrium of the nodal system, 
we get (see, for example, Bove [5], p. 12): 


81 Ej = ig Ey — — Bint, = 1h, 
— Bor Ey + be2 Ey — fan En = 9, 
ri Eni Ei — €ne Ls — + nn En = 9 
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where D is the system determinate, 4,, and Aj, are cofactors®) of D. Hence 


A 
aE ee ; 
Lig 
where 
821 83 San | 
ae — 831 833 — &3n H(y) 
12) oe | yt 
Ae an ae settee unt: con cos | 
esarian 5n3 Enn 


Here H(y) is a polynomial in y. 

On account of (2) it can be proved that the coefficients of H(y) are all real 
positive ®) (cf. [4], pp. 116, 124). 

It follows that the zeros of A,,, and hence of F,, are either real non-positive y’s 
or occur in conjugate complex pairs’). 

Refering to our variable 7 = 1/(1 + y), we find that the zeros of E,(%) are 
either real outside the interval 0 < x» < 1 or occur in conjugate complex pairs. 
This establishes property (1). 


4. Proof of Property (ii) 


We first examine the zeros and poles of the so-called driving-point and transfer 
resistances. 


4.1 The Driving- Point Resistance 


Consider a 2-terminal potentiometer network as shown in Figure 3. 


Figure 3 
2-terminal potentiometer network. 


The driving-point resistance is defined as the ratio 


E 
Rives 
1 qT, 


5) . . . . 
I, is a rational function of x. To discuss the zeros and poles of F,, we make use of 


energy considerations. If we consider the network upon loop system, then the 


°) The (tj-cofactor) = (—1)'+4(ij-mi ij-mi i i 
¢ j-minor), where the (ij-minor th i 
canceling the row i and the column j in D. Mme 


®) Note that 4,, and A,, may have common factors. When dividing A, [or H(y)] by such a_ 


factor, some coefficients of the numerator of E, may become negative 


{ 


7 Si argent 
) Zeros at the origin or at infinity may occur according to the degree of y in the denominator 


of Aj, and A,,. 
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energy loss or ‘dissipation function’ F is given by (see, for example, CAUER [6] 
p. 79, or GUILLEMIN [8], p. 5) 


P 
eb 

BF =>» rT 1;, (3) 
1,j=1 


where 


p number of loops in the system, 
I; current circulating in loop 7, 

ig total resistance in loop j, 

;j(¢@+7) resistance common to loops i and j. 
Of course 


(If the loop currents are all assumed positive in the same path direction, clockwise 
for example, then all 7; ;-terms are positive and all the r; ;-terms (i +) are negative.) 
On the other hand the energy dissipated in the network is known to be equal to 


Ry Ij, (4) 
so we get 
P 
Ra Tbs piety (5) 
f= 


In a loop system, the loop resistance is a series combination of our resistance 
elements 7,, x v, and (1 — %) 73, so that we can write 


%,5= 4, +%b,,;+(1—-*)¢,; (a4,6,c20 and 0S ¥S 1). (6) 
Substituting in (3) we get 
p p p 
Fa aj 1,1p + 4D) dig Tel + lo) Dy; 1, I, - (7) 
i,j=l i,j=1 ij=l 


The quadratic forms on the right side of (7) are positive definite *). Write 


p P Pp 
Dali =A, Voi l=B, Sel =C. (8) 


t,j=1 4,4=1 1,7=1 
Equating expression (7) to zero, we have 


C+A 
ype awed 9 
Parte. Ss (9) 
which is not an explicit relation of « since A, B and C are implicit functions in x, 
but the positive definiteness of A, B and C nevertheless shows that all zeros of F 


implied by (9) must be real negative or > 1°). Since from (5) 
Ry= (F eer 


it follows that the zeros of R,, are also real negative or > 1. 
8) This can be seen from their physical meaning since they represent energy (ci. [51] p: 128). 
Mathematically this can be proved by showing that the corresponding determinants and their 
principal minors are positive. For a criterion for positive definiteness see GUILLEMIN [9], p. 150. 
z ®) This is analogous to the treatment of LC, RL, and RC networks by GuiL_emin [8], pp. 51, 


58 and 64. 
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The energy (4) can be expressed by its equivalent form 


baie 
Rag 
which is again equal to F, hence 
1 
=2 (1B) hes 
Ruy ( )By=1 


This establishes that the poles of R,, are also real negative or > 1. 
To show that the zeros and poles of R,, are moreover simple, consider the 
voltage equilibrium equations of the system (see, for example [7], p. 227): 
ty Ty + tia Lg + ot pl = Eas 
tog Ly te Yagi lait +o tangle OG 


a, 6 %6 00 6) © 0 6 eee ‘e 6) @ © ee 18 iia: oe, Le) Ole, eam 


(10) 


Differentiating these equations with respect to « and indicating such differenti- 
ation by means of primes, we get 


V5 Tie te Te FE i 
Ye Tp tog TE tt tgp ty = Sy, (11) | 
Holy Xp y LEE oe ee 
where ! 
Sy = Pp ly Vg ds ae a Ope (12)q 
Solving (11) for J; we obtain | 
p 

—1 

i, = = Dee a. (13) | 


where D is the system determinant and 4,, are cofactors of D. 
But from (10) 


r= 21 Au 
a: D a 


Substituting for 4;, from (14) in (13) and noting that 


we get 


On the other hand 


Taree Bp (4) 


2 
Ri, 


. 
3 
| 
| 
| 

Hoes ae 

Ry 

from which by differentiation with respect to ¥ we have 

dx 
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Equating (15) and (16) we get 
Ey 


E,\? (dRu a 
(Bi) (= Sis. 
Substituting for S; from (12) we obtain 


Ey (aku) _ yo, 


afd 


Differentiating (6) with respect to x and substituting for y,; in (17) we get 


where B, C are represented by (8). 


Since the zeros of FR, are implied by (9), we see that for the internal zeros!) of 
Ry, (where C — B + 0) 


It follows that the zeros of R,, are simple (since at a zero of higher multiplicity 
the function must have zero slope). 

By the duality principle ([7], pp. 246-250) it follows that for the conductance 
function G,, = 1/R,, the zeros are also simple. This establishes that the zeros and 
poles of F,, are simple. 

The results of this part can thus be summarized in the property: the zeros and 
poles of the driving-point resistance R,, are simple, real negative or >1. 


4.2 The Transfer Resistance 


The aim of the following treatment is to show that the poles of the transfer 
resistance R,, (defined in 4.3) are contained in the poles of the driving-point 
resistance F4,. 


Figure 4 


4-terminal potentiometer network. 


If to any two nodes of the network of Figure 3 two terminals 2-2’ are fastened 
we get a 4-terminal network as shown in Figure 4. 

Consider this network upon a loop basis, where the voltages E, and E, are 
assumed as forming the closing loops to those meshes of the network which would 


10) Without considering the zeros at infinity. 
11) For negative zeros (C -B<0), dR,/dx 1s positive; whereas for zeros >1, the slope of the 


function R,, is negative. 


-ZAMP X/14 


po \ ie 
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be left open at the terminals 7—7’ and 2-2’ respectively. Designating these meshes 
as loops 7 and 2 and applying KircHHoFrr’s voltage law, we get the following 


system of equations (cf. [6], p. 3, and [7], pp. 134-136): 


ty Ty tg Ig to tha lp =, 
Vm fact fog tae oe Map I, =£,, 
¥o 1, + %92 Tg +o + %ap Lp = 9, 
Voi dy eV pata t “Fag lp = Oe 


where p, I;, 7;; denote the same quantities as in 4.1. 
Using CRAMER’Ss rule to solve this system for J, and J, we get 


I, = Gy Ey + Gy Eg, Tg = Gy Ey + Goo Ep, 
in which 


Aj, 
jb 2 


Gy =— >, Gpe=Ga= 


(D is again the system determinate and 4;, are cofactors of D; A;; = A;,). 


Solving (18) for E, and E, we have 
Ey=Ryl, + Ryl,, FE, = Ra li + Reels, 


where 
G G G 
R=, Rp = Ra ==... Rio = pe Gh GG 
1 |G| 12 21 |G| 22 |G| | 11 “22 12 
Substituting (19) in (21) we obtain 
Rite orege ee ee at 


Ay ° Ayo = Ais ‘ Ay ¢ Ago - Ais 


from which we see that a pole contained in R,, is also contained in R,. 


(20) 


(21) 


Since from 4.1 the poles of R,, are simple, real negative or >1, it follows that 


the poles of R,, are also simple, real negative or >1. 


4.3 Conclusion 


If in (20) we regard the end 2-2’ as open, then J, = 0, and we have 


EB 
Rife tay ee, ee 
- - Ey 


FR, is, therefore, defined as ‘open-circuit transfer’ resistance. From (22) 


Ror 


11 


E,=E} 


The poles of E, are thus the zeros of FR, and the poles of Ry. 


——— 


(22) 


After 4.1 and 4.2 we conclude that the poles of E, are simple, teal, negative 


Or #15 
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Since our 3-terminal network (Figure 1) maybe considered as a 4-terminal 
network with one input and one output terminal joined together, the above 
property of the poles of E, also applies to the 3-terminal networks, 

This establishes property (ii). 
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Zusammenfassung 


Betrachtet wird eine 3-Pol-Schaltung, die feste Widerstande und endlich viele 
lineare Potentiometer auf einer gemeinsamen Achse enthalt. Wenn diese Schal- 
tung von einer konstanten Spannung £, gespeist ist, dann ist die Ausgangsspan- 
nung £, an einem beliebigen Punkt der Schaltung eine rationale Funktion von 4, 
wobei x den Drehwinkel der Achse bedeutet (0 S$ 7 <1). 

In dieser Arbeit wird folgender Satz bewiesen: «Die Nullstellen von F, sind 
entweder reell ausserhalb des Intervalls 0 < * <1, oder konjugiert-komplex. Die 
Pole von E, sind einfach, reell, negativ oder >1.» 


(Received: August 11, 1958.) 
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Tagung fiir Elektronenmikroskopie in Freiburg i. Br. 
vom 18. bis 21. Oktober 1959 


Die diesjahrige Tagung der Deutschen Gesellschaft fiir Elektronenmikroskopie 


e. V. findet vom 18. bis 21. Oktober 1959 in Freiburg i. Br. statt. Anfragen sind zu 
“richten an Dr. H. KEHLeR, Farbwerke Hoechst AG, Frankfurt a. M.-Hoechst. 


Tea. 


H. KEHLER 


216 Buchbesprechungen — Book Reviews — Notices bibliographiques ZAMP 


Advances in Space Science 


Academic Press, Inc., Publishers, have completed plans for a new serial publica- 
tion entitled Advances in Space Science, for which FREDERICK I. Orpway has 
assumed editorial responsibility. Mr. ORDWway is associated with General Astronau- 
tics Corporation, Washington, D.C., and will be aided by an Editorial Advisory 
Board of international composition. 

The following have agreed to serve on this Board: WERNHER VON Braun (US), 
FREDERICK C. Durant (US), EuGEN SANGER (Germany), LesLiE R. SHEPHERD 
(Great Britain), GEorGE P. Sutton (US), and ETIENNE Vassy (France). 

Advances in Space Science will be devoted to critical reviews in the whole field 
of astronautics, with particular emphasis on the disciplines of mathematics, astron- 
omy, geophysics, geology, geography, and biology. 

L. HENLEIN — 
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Handbuch der Physik - Encyclopedia of Physics. Herausgegeben von © 
S. FLicce. Band 34: Korpuskeln und Strahlung in Materie II (Springer-Verlag, 
Berlin 1958). 316 S., 213 Fig.; DM 78.-. 

Es kann kaum die Aufgabe des Rezensenten sein, den Inhalt dieses reich- 
haltigen Bandes eingehend zu wiirdigen. 

Eréffnet wird die Folge der 6 Artikel durch den Beitrag von R. KOLLATH © 
iiber den Durchgang langsamer Elektronen und Ionen durch Gase. Es handelt 
sich dabei zum grossen Teil um klassische experimentelle Ergebnisse, an deren 
Gewinnung KoLiatH und RAMSAUER sehr massgebend beteiligt waren. 

Der zweite Artikel von 85 Seiten stammt von R. D. BirKHorFF und behandelt 
den Durchgang schneller Elektronen durch Materie. Dieses Gebiet ist vor allem — 
auch fiir den Strahlungsbiologen von Bedeutung. Neben der Diskussion der — 
Elementarprozesse (Streuung von Elektronen oder Positronen an Elektronen, 
Ubereinstimmung mit den Resultaten der Quantenelektrodynamik) und der 
klassischen Theorie der Bremsung von Elektronen wird auch der Stoss mit dem 
Plasma der Leitungselektronen behandelt. Nach einem Abschnitt iiber «stragg- 
ling» folgt die Diskussion der Streuung am Kern und abschliessend der Durch- . 
gang von Elektronen durch dicke Schichten (mit Reichweite-Energie-Relation). — 

Darauf folgt der Artikel von L. Simons (27 Seiten) iiber Positronium. Die 
Theorie des Positroniums findet sich im Artikel von G. KALLEN, Bd. 5, Teil 1. 
Hier werden die experimentellen Ergebnisse und Methoden dargestellt. Bekannt-_ 
lich ist die Ubereinstimmung der experimentellen Resultate mit der Theorie eine 
hervorragende Bestatigung der Quanten-Elektrodynamik. 

_ Der vierte Artikel von E, Merzpacuer und H. W. Lewis (27 Seiten) behandelt 
die theoretisch weit weniger durchsichtige Erzeugung von Réntgenstrahlen durch 
schwere geladene Teilchen. Er zerfallt naturgemass in einen Abschnitt, der die 


Ionisation der inneren Schalen behandelt in ei i i 
: , und in einen 
R6ntgenlicht diskutiert. meena Si 


ee 


eA” 
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Der Artikel von W. WHALING (22 Seiten) gibt hauptsachlich in Form von 
Kurven und Tabellen das experimentelle Material iiber die Bremsung von Pro- 
tonen und «-Teilchen fiir Energien unter 10 MeV. 

Der Schlussartikel von R. D. Evans (80 Seiten) gibt eine ausfiihrliche Dar- 
stellung des Compton-Effekts. Dieser hat, wie bekannt, bei der Entstehung der 
neuen Quantentheorie eine hervorragende Rolle gespielt. Es ist deshalb sehr zu 
begriissen, dass der vorliegende Artikel auch der historischen Seite des Gegen- 
standes gerecht wird. Eingehend werden daher die Fragen der Gleichzeitigkeit 
der Emission des Elektrons und die Folgerungen aus den Erhaltungssatzen 
behandelt, bevor die ausfiihrliche Diskussion der Klein-Nishina-Formel statt- 
findet. Es folgt dann ein Abschnitt iiber Compton-Schwachung mit sorgfaltigem 
Vergleich mit der Erfahrung und einer iiber den Compton-Effekt an gebundenen 
Elektronen, Den Abschluss bildet die in neuerer Zeit besonders wichtige Streuung 
an orientierten Elektronen. 

Dieser Handbuchband wird zweifellos vielen Fachleuten, zur zuverlassigen 
Orientierung dienen. Der Kreis der Interessenten erstreckt sich dabei iiber ein 
breites Spektrum, angefangen beim Theoretiker, der sich etwa iiber die experi- 
mentelle Verifikation der Quantenelektrodynamik orientieren will, iiber fast 
jegliche Art von Experimentalphysiker bis zum Strahlungsbiologen. 

Die Ausstattung des Bandes ist ausgezeichnet. R. Jost 


Introduction to Multivariate Statistical Analysis. Von T. W. ANDERSON. 
(John Wiley and Sons, New York 1948). 374 S, 16 Fig.; $12.50. 

Der bekannte amerikanische Statistiker gibt in diesem Lehrbuch eine aus- 
gezeichnete Darstellung der wichtigsten Methoden und Ergebnisse aus der Sta- 
tistik mehrerer Variablen. Ohne Kenntnis der entsprechenden Resultate fiir ein- 
dimensionale Statistiken ist dieses Buch nicht verstandlich. Der Verfasser benutzt 
den Matrizenkalkiil, um formal méglichst Anschluss an eindimensionale Darstel- 
lungen zu gewinnen. Neben einer ausfiihrlichen Besprechung der Korrelation, 
Likelihood-Methode im Falle mehrerer Variablen, Konfidenz-Bereiche, T?-Test, 
Unabhangigkeit, werden auch die zeitabhangigen Serien und stochastischen 
Gleichungen kurz besprochen. W. SAXER 


An Introduction to Probability Theory and Its Applications. Von. 
W. FELLER. Band 1, 2. Aufl. (John Wiley and Sons,. New York 1957). 461 S., 
12 Fig.; $10.75. 

Sechs Jahre nach Erscheinen der 1. Auflage dieses ausgezeichneten Lehr- 
buches fiir Wahrscheinlichkeitsrechnung ist bereits die 2. Auflage erschienen. 
Nach teilweise Umarbeitung einzelner Kapitel gibt der Verfasser noch wesentlich 
mehr. Anwendungen als in der 1. Auflage, vor allem aber solche aus der Physik. 
Wir erwdhnen insbesondere Kettenreaktionen, rekurrente Ereignisse, stocha- 
_ stische Prozesse, Irrfahrten (Brownsche Bewegung), ergodische Satze. 

Das Buch kann sowohl Mathematikern als auch andern Wissenschaftern zum 
Studium warm empfohlen werden. W. SAXER 


An Introduction to Combinatorial Analysis. Von J. Ricorpan (John 
Wiley & Sons, New York 1958). 244 S:.: § 8.50, ates 
Der Verfasser, Mitarbeiter der Bell Telephone Laboratories, gibt in diesem ori- 
-ginellen Buch eine gute Ubersicht iiber die Theorie der Kombinationen und ihrer 
Anwendungen. Neben den iiblichen bekannten, elementaren algebraischen Resul- 
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taten wird vor allem der Zusammenhang mit den erzeugenden Funktionen, der 
Theorie der Gruppen und Netzwerke ausfiihrlich dargestellt. Ganz besonders wird 
die bekannte Untersuchung von G. Porya iiber kombinatorische Anzahlbestim- 
mungen fiir Gruppen, Graphen und chemische Verbindungen einlasslich besprochen 
und weitergefitihrt. ee 

Neben den Mathematikern diirften sich vor allem Statistiker und Elektro- 
ingenieure fiir das Buch interessieren. Sehr viele Aufgaben erhéhen seinen Wert; 
es kann zum Studium warm empfohlen werden. W. SAXER 


Handbuch der Physik — Encyclopedia of Physics. Herausgegeben von 


S. Friicce, Band 28: Spektroskopie II (Springer-Verlag, Berlin 1957). VI+ 446 S., 
222 Fig.; DM 98.-. 


Der 28. Band des neuen Handbuches der Physik enthalt Beitrage tber Mikro-— 


wellenspektroskopie (W. Gorpy, 78 S.), kontinuierliche Spektren (W. FINKEN- 
BURG und TH. PETERS, 117 S.), Kristallgitterspektren (E. Fick und G. Joos, 89 S.), 
Zeeman-Effekt (VAN DEN Boscu, 32 S.) und natiirliche optische Aktivitat (J. P. 
MATHIEU, 106 S.). 

Gorpy folgt dem Schema friiherer Darstellungen: experimentelle Technik, 
Molekiil-, Atom- und paramagnetische Elektronenresonanz-Spektren im Mikro- 
wellengebiet sind im Sinne einer Ubersicht beschrieben. Daher fehlt die Wieder- 
gabe sowohl allgemein wichtiger Ergebnisse (wie zum Beispiel Tabellen der Spek- 
tren asymmetrischer Kreisel ohne und mit Quadrupolkernen, Kreisel mit interner 


Rotation usw.) wie auch vieler spezieller Untersuchungen, die fiir ein Handbuch ~ 


dieses Ausmasses angemessen waren. Die Literatur ist bruchstiickweise zitiert. 
Der Beitrag tiber kontinuierliche Spektren zeichnet sich durch knappe Fassung 
und Ubersichtlichkeit aus. Die im ersten Abschnitt zusammengestellte Theorie wird 


systematisch zur Beurteilung kontinuierlicher Emissions- und Absorptions-Elek- — 


tronenspektren von Atomen und Ionen, zum Beispiel in Bogenentladungen und 
Schockwellen, sowie freier Elektronen herangezogen. Die umfangreichen empiri- 
schen Kenntnisse tiber Molekiilkontinua, hauptsdchlich zweiatomige, sind zum Teil 
nach den Potentialkurven, zum Teil nach Elementen und nach Erzeugungs- 


methoden geordnet. Das Literaturverzeichnis zu diesem Beitrag scheint bemerkens- — 


wert vollstandig. 
Das Kapitel iiber Kristallspektren ist beschrankt auf die Spektren einheitlicher 
und gemischter Ionenkristalle, homéopolarer Kristalle und Molekiilkristalle. Der 


klassischen Theorie der Termaufspaltung im Kristallfeld von BETHE (mit tabella-_ 


risch wiedergegebenen Resultaten) geht eine kurze gruppentheoretische Einleitung 
voraus. Neben der Kramers-Entartung findet die Methode der Elektroneneigen- 
funktion des Kristalls (HELLWEGE) breite Behandlung, wobei einige Resultate fiir 
elektrische und magnetische Dipoliibergange tabelliert sind. Ferner ist der Uber- 
gang schwaches > mittleres > starkes Kristallfeld (BETHE folgend) im Zusammen- 
hang mit den Racah-Koeffizienten, den Stevenschen Matrixelementen des Kristall- 
felds und dem Zeeman-Effekt naher diskutiert. Die neueren Untersuchungen tiber 
die Linienbreite im Kristall werden nur gestreift. Letzteres gilt auch fiir die Terme 
unpolarer Molekiilkristalle, zum Beispiel die neueren Arbeiten iiber «site symmetry 
groups» und «libration» im «orientated gas»-Modell. Die Diskussion experimenteller 
Daten ist im wesentlichen beschrankt auf Ionenkristalle. 

Der schone Beitrag iiber den Zeeman-Effekt ist hauptsachlich eine Zusammen- 
stellung der theoretischen Ergebnisse, umfassend schwaches Feld und die drei 
Kopplungsfalle, starkes Feld, Regeln fiir den g-Faktor, Intensitat der Zeeman- 
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Linien, Hyperfeinstruktur und die Analyse beobachteter Zeeman-Spektren von 
Atomen. In diesem Sinne ist der Artikel als vollstandig zu bezeichnen. 

Im Sinne einer allgemeinen Ubersicht hat Matureu die optische Aktivitat be- 
handelt: Ausschnitte iiber phanomenologische Beschreibung, experimentelle Me- 
thodik (unvollstandig), Stereochemie (an Beispielen), klassische und quanten- 
mechanische Theorien (ausfiihrlich) leiten zur Diskussion empirischer Gesetz- 
massigkeiten des natiirlichen Drehvermégens isotroper und anisotroper Kérper 
uber. Wahrend man hier wenig bekannte Arbeiten iiber die Wirkung dusserer 
Fehler auf das Drehvermégen (KALLEMANN) findet, fehlt zum Beispiel die Bespre- 
chung moderner wichtiger Theorien der optischen Aktivitat schraubenlinien- 
formiger Molekeln. 

Der 28. Band des Handbuches halt sich in bezug auf Individualitat und Stil 
der einzelnen Beitrage im Rahmen anderer Bande des Werkes. 

Hs. H. GUNTHARD 


Physik und Technik der Atomreaktoren. Von F. Cap (Springer-Verlag, 
Wien 1957). 487 S., 100 Abb.; sFr. 65.50. 

Das Buch Physik und Technik der Atomreaktoren richtet sich an Studierende 
héherer Semester und an Reaktorfachleute. Es will nicht nur als Lehrbuch, sondern 
auch als Nachschlagewerk dienen. Das Werk ist in neun Kapitel eingeteilt, in wel- 
chen die kernphysikalischen Grundlagen, die Bremsung der Neutronen, die Diffu- 
sion der Neutronen, die Theorie des homogenen Reaktors, die Theorie des heterogenen 
Reaktors, der Bau von Reaktoren, der Betrieb von Reaktoren, die Beschreibung 
von verschiedenen Reaktortypen und die Verwendung von Reaktoren behandelt 
werden. Die theoretischen Kapitel des Buches sind klar geschrieben und bereiten 

einem beim Durcharbeiten Freude. Die beigefiigten Ubungsbeispiele erganzen in 
schéner Weise den gebotenen Stoff. Leider sind die Resultate nicht angegeben. Bei 
einzelnen Aufgaben ware ein kurzer Hinweis auf den Loésungsgang sehr wertvoll. 
Obwohl das Buch als Reaktorlehrbuch in deutscher Sprache empfohlen werden 
kann, mdchten wir doch nicht auf die amerikanischen Lehrbiicher wie etwa: 
GLASSTONE und Epitunp, The Elements of Nuclear Reactor Theory; MuRRAY, 
Nuclear Reactor Physics oder das neu erschienene Werk von WEINBERG und WIG- 
NER, The Physical Theory of Nucleay Chain Reactors verzichten. W. WINKLER 


Differential Equations: Geometric Theory. Yon Soromon LEFSCHETZ 
(Interscience Publ., Inc., New York 1957), X+364 S., 73 Fig.; $ 9.50. 

Das in der wohlbekannten Reihe «Pure and Applied Mathematics» als Band VI 
erschienene Buch von S. LerscHetz stellt eine wesentliche Erweiterung des frii- 
heren Werkes Lectures on Differential Equations (1946) dar. Die Darstellung macht 
weitgehenden Gebrauch von der Vektor- und Matrizenschreibweise sowie von 
topologischen Begriffen, wobei das Wichtigste dazu in den einleitenden Abschnitten 
sowie in einem Anhang geboten wird. Einen breiten Raum nimmt die Behandlung 
der nichtlinearen Systeme von der Form dv/dt = P x + q(x; ¢) ein, wo P eine kon- 
stante oder von ¢ abhangige Matrix ist. Eingehend werden die so wichtigen Fragen 

der Stabilitat der Lésungen und der Periodizitat behandelt. Speziell werden hierauf 
die zweidimensionalen Systeme ausfiihrlich im Hinblick auf ihre isolierte Singulari- 
taten und ihr Verhalten im Unendlichen untersucht. Die beiden letzten Kapitel 
‘schliesslich bringen noch eine Reihe moderner Ergebnisse zu den Differentialglei- 
‘chungen 2. Ordnung und den nichtlinearen Schwingungen. — Die Darstellung der 
“Materie ist klar und iibersichtlich, und besondern Dank wird man dem Verfasser 
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_ Beriicksichtigung des erhéhten Risikos bei der Bemessung der Praémie oder bei 


 zweite Buch nicht nur dem Nachwuchs, sondern auch den in der Praxis stehenden 
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wissen fiir die eingehende Beriicksichtigung der wichtigen russischen Arbeiten aus 
der jiingsten Zeit, die ja haufig nur sehr schwer zuganglich sind. Das Buch wendet | 
sich nicht allein an den Mathematiker, sondern vor allem auch an den Physiker 
und Ingenieur, fiir den die behandelten Probleme von grosser Bedeutung sind. 


Diesem empfehlenswerten Band ist eine weite Verbreitung sicher. 
E. RotH-DESMEULES 


Versicherungsmathematik, 2. Teil. Von WALTER SAXER, mit einem An- 
hang von HEINRICH JECKLIN (Springer-Verlag, Berlin 1958). 283 Seiten; DM 45.-. 

Den ersten Teil seiner Versicherungsmathematik, der bereits vor drei Jahren 
erschienen ist, hat der Verfasser auf die Bediirfnisse des Versicherungsmathema- 
tikers in seiner praktischen Tatigkeit ausgerichtet, indem er mit grosser Sorgfalt 
die fiir diesen Aufgabenbereich ben6tigten theoretischen Grundlagen in metho- — 
disch durchdachter und straff begriindeter Form bereitstellte. Yon Anfang an ~ 
plante er, sein Werk zu vervollstandigen durch einen zweiten Teil, der tiber den — 
iiblichen Rahmen der Lehrbiicher fiir Versicherungsmathematik hinausreicht, — 
mit Problemstellungen und Lésungen, die bisher nur in Spezialabhandlungen ein- 
schlagiger Zeitschriften behandelt worden sind. Dieser von den Fachvertretern — 
mit Spannung erwartete zweite Teil des Saxerschen Lehrbuches ist soeben er- 
schienen. : . 

Das Buch geht aus von einer «Theorie der Versicherungsfunktionen». Die von 
STIELTJES und SCHARF ausgearbeiteten Integral-Begriffe ermdglichen eine ein- 


heitliche Darstellung der meist nur stufenweise beobachteten Entwicklungen und ~ 


der mit Vorteil herangezogenen mathematischen Funktionen. Einer konzentriert — 
durchgefiihrten Darstellung der «Mathematik der Lebensversicherung» folgt eine — 
noch in keinem Lehrbuch dargelegte «Mathematik allgemeiner Risikoversiche- — 
rungen», die sich nicht auf das Gebiet der Lebensversicherung beschrankt. Als — 
mathematische Hilfsmittel zur Erfassung stochastischer Vorgange werden die 
wichtigsten Verteilungsfunktionen beschrieben. Zur Einfiihrung in die weitschich- 
tigen Probleme, mit denen sich die Risikotheorie befasst, werden die Bestimmung 
der Pramie bei verschiedener Abgrenzung der Leistungspflicht sowie die Auswir- 
kung von Schwankungsreserven und die Ermittlung der Risikogewinne behandelt. 

Ein weiteres, in sich abgeschlossenes Kapitel ist der «Erneuerungstheorie» 
offener Gesamtheiten gewidmet, deren Elemente nach bestimmten Gesetzen aus-_ 
scheiden und durch neue ersetzt werden. Wahrend die Gewinnung von Sterbe- 
tafeln aus statistischen Erhebungen schon im ersten Teil des Lehrbuches behan- 
delt wird, befasst sich der zweite Teil mit der « Ausgleichung von Sterbetafeln», 
indem er praktisch bewahrte Verfahren und Teste zu deren Beurteilung schildert. 

Zum zweiten Teil des Lehrbuches von Saxer fiigt JeEcKLIN einen Anhang 
hinzu, der sich der technischen Behandlung erhéhter Risiken in der Lebensver- 
sicherung widmet und die Gewinnung: der statistischen Grundlagen sowie die 


— 


der Beschrankung der Versicherungsleistung darstellt. 

Die einzelnen Kapitel des zweiten Teiles der Versicherungsmathematik von 
SAXER lassen sich auch unabhangig voneinander lesen. Insgesamt bringen sie eine 
Vertiefung der im ersten Buch dargelegten Kenntnisse. Von besonderer Bedeutung 
scheint uns die Ausweitung des Blickes auf die Fiille der Probleme, vor die sich 
heutzutage der Versicherungsmathematiker gestellt sieht. Daher diirfte dieses 


sn 


bewahrten Versicherungsmathematikern reiche Anregung und wertvolle Stiitze 
bieten. ; : H. Wyss 
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SELECTA HERMANN WEYL 


Herausgegeben. zu seinem siebzigsten Geburtstag 
von der Eidgenéssischen Technischen Hochschule in Ziirich 
und vom Institute for Advanced Study in Princeton 


(2956) 592 Seiten, Ganzleinen Fr./DM 48.90 


B. L. vAN DER WAERDEN in «Zeitschrift fiir angewandte Mathematik und Physik»: «In der immer an- 
schwellenden Flut der Fachliteratur ist es sehr wichtig, die wertvollsten Arbeiten der wahrhaft grossen 
Forscher in Buchform vereinigt zu haben. Zu diesen Grossen gehérte WEYL.» 


H. KNeEser in iJariethericht der Deutschen Mathematikervereinigung»: Ein schénes Bildnis Weyls « er: 

dffnet den Band; am Schluss zihlt ein Verzeichnis seiner Veréffentlichungen 16 teilweise mehrfach aut- 

gelegte und iibersetzte Biicher und Einzelschriften sowie 163 Aufsitze in Zeitschriften auf. Was da- 

zwischen steht, hat die Mannigfaltigkeit und Geschlossenheit eines Kunstwerks, auch in dieser Hinsicht 
_ ein Bild der einzigartigen Persénlichkeit Hermann Weyl» 


‘When the greatest mathematician of his generation died on the 8th December, 1955, this magnificent 

volume in honour of his 7oth birthday, produced jointly by the Zurich’ Federal Institute of Technology 

and the Princeton Institute of Advanced Study, was passing through the press. 

The paper and printing are of a fine quality all too rarely seen today and worthy of the material they 
display».’ R. L. GoopstEIn. 
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Von HERMANN WEYL 
- Professor an der Eidg. Technischen Hochschule in Ziirich und am Institute for Advanced Study in Princeton. 


Aus dem Englischen tibersetzt von Lutu BecHTOLsHEIM, Professorin an der University of Redlands (USA) 
(1955) 160 Seiten mit 72 Figuren. In Ganzleinen Fr. 17.70 (DM 17.70) 


Inwat: 1. Bilaterale Symmetrie —2. Translatorische, rotative und verwandte Symmetrien — 3. Or- 
mamentale Symmetrie — 4. Die allgemeine mathematische Idee der Symmetrie — 5. Zwei Anhange iber 
endliche Rotationsgruppen in drei Dimensionen 


‘ 


O.-H, Ketzer in «ZAMM»: «Das Buch wendet sich weniger an den mathematischen Spezialisten; der 

Mathematiker wird aber durch dieses Buch mit grossem Gewinn die Verwirklichung der mathematischen 

Dinge erleben kénnen, wie umgekehrt der Nichtfachmann den mathematischen Kern in den Erschei- 

- nungen zu Gesicht bekommt. Dem Buch ist eine weite Verbreitung in den Kreisen der gebildeten Fach- 
leute und der gebildeten Nichtfachleute zu wiinschen.» 


« Mathematisch-Physikalische Semesterberichte»: «Das Buch ist fiir jeden Mathematiklehrer an hoheren 
Schulen, der sich um die Durchdringung seines Unterrichts mit abbildungsgeometrischen bie tes a 
bemiiht, eine Fundgrube fiir wertvolle Anregungen. » 


B. Huppert in ¢ Jahresbericht der Deutschen Mathematiker-Vereinigung»: «Das Biichlein zeigt nicht 

nur den grossen Mathematiker Hermann Weyl in seiner vollendeten Kunst der Darstellung, es gibt mit 

der Fiille der Beispiele aus Biologie und bildender Kunst auch einen Eindruck von der Vielseitigkeit 
Hermann Weyls.» 
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